# Specific Domains

**URL:** https://discourse.julialang.org/c/domain/10.md?page=262

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

**Page:** 263

---

## [Broadcasting an expression-evaluator function](https://discourse.julialang.org/t/broadcasting-an-expression-evaluator-function/85127)

<div class="topic-metadata">

**Author:** [@RXGottlieb](https://discourse.julialang.org/u/RXGottlieb)\
**Replies:** 2\
**Last updated:** [August 2, 2022, 2:21pm UTC](https://discourse.julialang.org/t/broadcasting-an-expression-evaluator-function/85127 "2022-08-02T14:21:46Z")

</div>

Hello, I have a symbolic expression that I’d like to evaluate a large number of times on a GPU. A few days ago I made a post about creating an evaluator function (see: the previous post), which led me to Symbolics.build…

---

## [Linearizing ModelingToolkit ODESystem](https://discourse.julialang.org/t/linearizing-modelingtoolkit-odesystem/61387)

<div class="topic-metadata">

**Author:** [@baggepinnen](https://discourse.julialang.org/u/baggepinnen)\
**Replies:** 3\
**Last updated:** [August 2, 2022, 11:58am UTC](https://discourse.julialang.org/t/linearizing-modelingtoolkit-odesystem/61387 "2022-08-02T11:58:48Z")

</div>

I am having my bi-monthly go at ModelingToolkit and have reached a point where I can create an ODESystem from ControlSystems.jl StateSpace and connect a few of those togeher and simulate, nice :slight\_smile: As a next s…

---

## [Agent Based Modeling resources other than Agents.jl documentation](https://discourse.julialang.org/t/agent-based-modeling-resources-other-than-agents-jl-documentation/85035)

<div class="topic-metadata">

**Author:** [@franfram](https://discourse.julialang.org/u/franfram)\
**Replies:** 2\
**Last updated:** [August 2, 2022, 8:39am UTC](https://discourse.julialang.org/t/agent-based-modeling-resources-other-than-agents-jl-documentation/85035 "2022-08-02T08:39:32Z")

</div>

Hello everyone. I’m starting to use Agents.jl for a PhD project on animal movement and epidemiology. I was wondering if anyone knew of additional resources to learn other than those of Agents.jl. George Datseris youtu…

---

## [Nested data structure with value by reference](https://discourse.julialang.org/t/nested-data-structure-with-value-by-reference/85147)

<div class="topic-metadata">

**Author:** [@Thomberger](https://discourse.julialang.org/u/Thomberger)\
**Replies:** 0\
**Last updated:** [August 2, 2022, 8:24am UTC](https://discourse.julialang.org/t/nested-data-structure-with-value-by-reference/85147 "2022-08-02T08:24:54Z")

</div>

I would like to have nested structure with minimum memory allocation. This would need passing data by reference but i’m not able to enforce that explicitly. I have a small number (i.e. 5) of different geometries. Each g…

---

## [Is there any way to optimize array additions and multiplications with transposes?](https://discourse.julialang.org/t/is-there-any-way-to-optimize-array-additions-and-multiplications-with-transposes/80268)

<div class="topic-metadata">

**Author:** [@Geositta](https://discourse.julialang.org/u/Geositta)\
**Replies:** 27\
**Last updated:** [August 2, 2022, 8:06am UTC](https://discourse.julialang.org/t/is-there-any-way-to-optimize-array-additions-and-multiplications-with-transposes/80268 "2022-08-02T08:06:32Z")

</div>

The question is related to python - Optimizing array additions and multiplications with transposes - Stack Overflow Given A is an array, I want to evaluate sum += p(A) \* factor(p) where p is a permutation/transpose, fa…

---

## [Are there any interesting visualization examples in Julia？](https://discourse.julialang.org/t/are-there-any-interesting-visualization-examples-in-julia/62661)

<div class="topic-metadata">

**Author:** [@Wind](https://discourse.julialang.org/u/Wind)\
**Replies:** 16\
**Last updated:** [August 2, 2022, 5:37am UTC](https://discourse.julialang.org/t/are-there-any-interesting-visualization-examples-in-julia/62661 "2022-08-02T05:37:48Z")

</div>

I found a web page with a lot of dazzling and interesting javascript visualization examples. It’s true that programming language is the best drawing tool. Since Julia is very fast and easy to use, I assume it can do a l…

---

## [Should model matrix for nested factors be full-rank?](https://discourse.julialang.org/t/should-model-matrix-for-nested-factors-be-full-rank/84733)

<div class="topic-metadata">

**Author:** [@PharmCat](https://discourse.julialang.org/u/PharmCat)\
**Replies:** 15\
**Last updated:** [August 1, 2022, 8:45pm UTC](https://discourse.julialang.org/t/should-model-matrix-for-nested-factors-be-full-rank/84733 "2022-08-01T20:45:00Z")

</div>

In some cases when levels of one factor all “inside” levels of another we have not full-rank model matrix from StatsModels. For example - calssical bioequivlence crossover study (2X2) with factors: Var ~ Formulation + …

---

## [Makie inconsistency between GLMakie output and saved mp4 output](https://discourse.julialang.org/t/makie-inconsistency-between-glmakie-output-and-saved-mp4-output/85112)

<div class="topic-metadata">

**Author:** [@gottacatchenall](https://discourse.julialang.org/u/gottacatchenall)\
**Replies:** 0\
**Last updated:** [August 1, 2022, 6:37pm UTC](https://discourse.julialang.org/t/makie-inconsistency-between-glmakie-output-and-saved-mp4-output/85112 "2022-08-01T18:37:41Z")

</div>

Hi, I’m running the following code to create an animation with two axes within a figure. using Makie, GLMakie using Distributions fig = Figure(resolution=(1600,800)) ax1 = Axis(fig\[1,1\], title="Sampling Locations"…

---

## [Symbolics.jl | Lie series expansion for PDEs](https://discourse.julialang.org/t/symbolics-jl-lie-series-expansion-for-pdes/84885)

<div class="topic-metadata">

**Author:** [@Audrius-St](https://discourse.julialang.org/u/Audrius-St)\
**Replies:** 2\
**Last updated:** [August 1, 2022, 7:43pm UTC](https://discourse.julialang.org/t/symbolics-jl-lie-series-expansion-for-pdes/84885 "2022-08-01T19:43:31Z")

</div>

Hello, I’ve used Sympy for several years now and am interested in giving Symbolics.jl a go with an eye to switching over completely. I’ve previously used Sympy to solve systems of ODES using the symbolic Lie series met…

---

## [Euler Equation Iteration to solve an Optimization problem](https://discourse.julialang.org/t/euler-equation-iteration-to-solve-an-optimization-problem/84135)

<div class="topic-metadata">

**Author:** [@Albert\_Zevelev](https://discourse.julialang.org/u/Albert_Zevelev)\
**Replies:** 21\
**Last updated:** [August 1, 2022, 7:20pm UTC](https://discourse.julialang.org/t/euler-equation-iteration-to-solve-an-optimization-problem/84135 "2022-08-01T19:20:11Z")

</div>

I would like to solve a simple/routine economics problem with Euler Equation Iteration. The following code works well: using Interpolations, Plots, LinearAlgebra # PolyBasis(K) = \[ones(size(K)) K K.^2\] PolyGetCoef…

---

## [Plots.jl contour line colors](https://discourse.julialang.org/t/plots-jl-contour-line-colors/85084)

<div class="topic-metadata">

**Author:** [@gideonsimpson](https://discourse.julialang.org/u/gideonsimpson)\
**Replies:** 3\
**Last updated:** [August 1, 2022, 6:40pm UTC](https://discourse.julialang.org/t/plots-jl-contour-line-colors/85084 "2022-08-01T18:40:09Z")

</div>

When using contour (not contourf) with Plots.jl, is there a way to force all the contour lines to be a single fixed color?

---

## [Modelling uncertainty propagation in forecasts 🔮](https://discourse.julialang.org/t/modelling-uncertainty-propagation-in-forecasts/85089)

<div class="topic-metadata">

**Author:** [@agucova](https://discourse.julialang.org/u/agucova)\
**Replies:** 1\
**Last updated:** [August 1, 2022, 2:58pm UTC](https://discourse.julialang.org/t/modelling-uncertainty-propagation-in-forecasts/85089 "2022-08-01T14:58:21Z")

</div>

Hi everyone, I’ve been playing with Causal, a business modeling platform which does automatic uncertainty calculations in forecasts. This works by allowing one to set model parameters in terms of ranges, which are real…

---

## [Kriging on surface of sphere](https://discourse.julialang.org/t/kriging-on-surface-of-sphere/85058)

<div class="topic-metadata">

**Author:** [@yakir12](https://discourse.julialang.org/u/yakir12)\
**Replies:** 5\
**Last updated:** [August 1, 2022, 1:47pm UTC](https://discourse.julialang.org/t/kriging-on-surface-of-sphere/85058 "2022-08-01T13:47:39Z")

</div>

I have a few random points on a sphere’s surface and I’d like to interpolate and extrapolate those using kriging. But I’m getting artefacts in the result that look like an issue with the fact that the coordinates are sph…

---

## [Vector to lists of indices grouped by key](https://discourse.julialang.org/t/vector-to-lists-of-indices-grouped-by-key/85050)

<div class="topic-metadata">

**Author:** [@FireCrumb](https://discourse.julialang.org/u/FireCrumb)\
**Replies:** 9\
**Last updated:** [August 1, 2022, 11:00am UTC](https://discourse.julialang.org/t/vector-to-lists-of-indices-grouped-by-key/85050 "2022-08-01T11:00:46Z")

</div>

Hi, I have a vector of numbers: (in reality much larger) arr = Vector{Int64}(\[6, 9, 9, 4, 1, 1, 2, 7, 8, 3\]) I want to be able to generate the list of indices for each value: \> p = pairs(arr) pairs(::Vector{Int64})(.…

---

## [Makie.jl figure with broken axis](https://discourse.julialang.org/t/makie-jl-figure-with-broken-axis/85062)

<div class="topic-metadata">

**Author:** [@Kevin\_Kleinbeck](https://discourse.julialang.org/u/Kevin_Kleinbeck)\
**Replies:** 3\
**Last updated:** [August 1, 2022, 7:09am UTC](https://discourse.julialang.org/t/makie-jl-figure-with-broken-axis/85062 "2022-08-01T07:09:05Z")

</div>

I intend to create a figure with a broken Axis, which can be useful when the data range is relatively large but relevant details occur only close to the maximal / minimal values. For example, let’s say we want to create…

---

## [GPU performance and switching from tabular to recurrent data format for Flux.jl](https://discourse.julialang.org/t/gpu-performance-and-switching-from-tabular-to-recurrent-data-format-for-flux-jl/84969)

<div class="topic-metadata">

**Author:** [@JLDC](https://discourse.julialang.org/u/JLDC)\
**Replies:** 5\
**Last updated:** [August 1, 2022, 6:52am UTC](https://discourse.julialang.org/t/gpu-performance-and-switching-from-tabular-to-recurrent-data-format-for-flux-jl/84969 "2022-08-01T06:52:48Z")

</div>

When training a recurrent network using Flux.jl, a dataset with K features, S samples, and a sequence length of L should take the form of a vector of length L where each element is a K × S matrix. Let’s say that I have …

---

## [Ensemble simulations with random initial conditions that sum to a fixed value](https://discourse.julialang.org/t/ensemble-simulations-with-random-initial-conditions-that-sum-to-a-fixed-value/84951)

<div class="topic-metadata">

**Author:** [@Dragolonth](https://discourse.julialang.org/u/Dragolonth)\
**Replies:** 6\
**Last updated:** [August 1, 2022, 6:51am UTC](https://discourse.julialang.org/t/ensemble-simulations-with-random-initial-conditions-that-sum-to-a-fixed-value/84951 "2022-08-01T06:51:59Z")

</div>

Hi all. For my ensemble problem, I would like to randomize the initial conditions for each trajectory. However, I would like the sum of certain species to remain a fixed number. For example, if u0 is my set of initial co…

---

## [BifurcationKit MethodError](https://discourse.julialang.org/t/bifurcationkit-methoderror/84904)

<div class="topic-metadata">

**Author:** [@Dragolonth](https://discourse.julialang.org/u/Dragolonth)\
**Replies:** 3\
**Last updated:** [August 1, 2022, 6:14am UTC](https://discourse.julialang.org/t/bifurcationkit-methoderror/84904 "2022-08-01T06:14:37Z")

</div>

Hi all. Currently, I have a reaction network defined in Catalyst, and I’ve been following the tutorial at Bifurcation Diagrams · Catalyst.jl. However, I’ve been running into an error and am not entirely sure how to resol…

---

## [Mixed-Integer Multi-objective Constrained Optimization Problem Solving](https://discourse.julialang.org/t/mixed-integer-multi-objective-constrained-optimization-problem-solving/83030)

<div class="topic-metadata">

**Author:** [@Magy](https://discourse.julialang.org/u/Magy)\
**Replies:** 4\
**Last updated:** [July 31, 2022, 11:43pm UTC](https://discourse.julialang.org/t/mixed-integer-multi-objective-constrained-optimization-problem-solving/83030 "2022-07-31T23:43:51Z")

</div>

Dear All, I wonder if there is a direct package that can be used to solve a mixed-integer multi-objective constrained problem. I tried to use “Metaheuristics.jl” but failed to represent my model’s mixed-integer variable…

---

## [Allowing Agents.jl OpenStreetMapSpace to read GeoJSON](https://discourse.julialang.org/t/allowing-agents-jl-openstreetmapspace-to-read-geojson/85055)

<div class="topic-metadata">

**Author:** [@imantha](https://discourse.julialang.org/u/imantha)\
**Replies:** 1\
**Last updated:** [July 31, 2022, 7:15pm UTC](https://discourse.julialang.org/t/allowing-agents-jl-openstreetmapspace-to-read-geojson/85055 "2022-07-31T19:15:40Z")

</div>

Hi I was wondering whether it would be too difficult add GeoJSON reading capabilities to agents.jl’s OpenStreetMapSpace. I am asking this as open source tools such as QGIS allows easy conversion between shp and GeoJSON …

---

## [Solve many first-order ODEs](https://discourse.julialang.org/t/solve-many-first-order-odes/85060)

<div class="topic-metadata">

**Author:** [@qwerty](https://discourse.julialang.org/u/qwerty)\
**Replies:** 6\
**Last updated:** [July 31, 2022, 6:25pm UTC](https://discourse.julialang.org/t/solve-many-first-order-odes/85060 "2022-07-31T18:25:21Z")

</div>

What is the best way to solve so many decoupled ODEs? Possibly in parallel. Should I use Modeling Toolkit or DifferentialEquations? The equations are of the type: u’ = -a^2k^2u where a is the same for all, then I have…

---

## [Any conjugate gradient method allowing a custom inner product?](https://discourse.julialang.org/t/any-conjugate-gradient-method-allowing-a-custom-inner-product/84982)

<div class="topic-metadata">

**Author:** [@ranocha](https://discourse.julialang.org/u/ranocha)\
**Replies:** 5\
**Last updated:** [July 30, 2022, 5:57am UTC](https://discourse.julialang.org/t/any-conjugate-gradient-method-allowing-a-custom-inner-product/84982 "2022-07-30T05:57:01Z")

</div>

Is there any conjugate gradient (CG) implementation in Julia allowing one to specify the inner product (with a fixed symmetric and positive definite matrix M, say)? In my use case, I know my system matrix A is symmetric…

---

## [Control of Time-Step size in OrdinaryDiffEq](https://discourse.julialang.org/t/control-of-time-step-size-in-ordinarydiffeq/84949)

<div class="topic-metadata">

**Author:** [@Whyborn](https://discourse.julialang.org/u/Whyborn)\
**Replies:** 4\
**Last updated:** [July 29, 2022, 11:39pm UTC](https://discourse.julialang.org/t/control-of-time-step-size-in-ordinarydiffeq/84949 "2022-07-29T23:39:52Z")

</div>

Hi, I am trying to use OrdinaryDiffEq for the time integration of a hyperbolic PDE system in a finite-volume context. I want to control the size of the time-step through a CFL condition, rather than by some error norm. …

---

## [Nonlinear curve fitting in Julia](https://discourse.julialang.org/t/nonlinear-curve-fitting-in-julia/84999)

<div class="topic-metadata">

**Author:** [@Jian\_ZUO](https://discourse.julialang.org/u/Jian_ZUO)\
**Replies:** 4\
**Last updated:** [July 29, 2022, 9:20pm UTC](https://discourse.julialang.org/t/nonlinear-curve-fitting-in-julia/84999 "2022-07-29T21:20:33Z")

</div>

Hi all, My question is to do a nonlinear curve fitting task (bounded) for estimating 10 parameters from several give data points, presented as follows: using Optim # my data points I = \[100.254, 94.261, 88.224, 82.193,…

---

## [JuMP 2-d array Decision Variable computationally expensive](https://discourse.julialang.org/t/jump-2-d-array-decision-variable-computationally-expensive/84987)

<div class="topic-metadata">

**Author:** [@sophiepavia](https://discourse.julialang.org/u/sophiepavia)\
**Replies:** 1\
**Last updated:** [July 29, 2022, 6:37pm UTC](https://discourse.julialang.org/t/jump-2-d-array-decision-variable-computationally-expensive/84987 "2022-07-29T18:37:00Z")

</div>

I am creating a decision variable that is a 2-d array of size 25122 x 25122. (length of arcs and od are both 25122) Is there a more efficient way to initialize a multi-dimensional array as a decision variable? I have ch…

---

## [Siamese/Twin network in Flux (working with two image inputs)](https://discourse.julialang.org/t/siamese-twin-network-in-flux-working-with-two-image-inputs/84991)

<div class="topic-metadata">

**Author:** [@LetsGrillBananas](https://discourse.julialang.org/u/LetsGrillBananas)\
**Replies:** 2\
**Last updated:** [July 29, 2022, 6:19pm UTC](https://discourse.julialang.org/t/siamese-twin-network-in-flux-working-with-two-image-inputs/84991 "2022-07-29T18:19:36Z")

</div>

I’m trying Flux for the first time, so I’m porting a project I did in PyTorch earlier this year over to Flux. Essentially, the project is to classify plant leaves as healthy or diseased. I already have a baseline transfe…

---

## [When plotting Optimization.jl callbacks with Plots.jl, use just one figure in VS Code plots paned](https://discourse.julialang.org/t/when-plotting-optimization-jl-callbacks-with-plots-jl-use-just-one-figure-in-vs-code-plots-paned/84948)

<div class="topic-metadata">

**Author:** [@thompsonmj](https://discourse.julialang.org/u/thompsonmj)\
**Replies:** 3\
**Last updated:** [July 29, 2022, 2:28pm UTC](https://discourse.julialang.org/t/when-plotting-optimization-jl-callbacks-with-plots-jl-use-just-one-figure-in-vs-code-plots-paned/84948 "2022-07-29T14:28:16Z")

</div>

In the callback function below cb = function(θ,l,pred) … how can I have only the current prediction appear on a plot with each iteration to illustrate the progress of fitting? I have tried putting fig = plot() inside th…

---

## [Growing cell population using callback](https://discourse.julialang.org/t/growing-cell-population-using-callback/84970)

<div class="topic-metadata">

**Author:** [@hbooth](https://discourse.julialang.org/u/hbooth)\
**Replies:** 10\
**Last updated:** [July 29, 2022, 2:11pm UTC](https://discourse.julialang.org/t/growing-cell-population-using-callback/84970 "2022-07-29T14:11:33Z")

</div>

Hi everyone, I am building a model of a growing cell population. I want to make use of callback to evaluate the ‘division’ event, similar to the example in this tutorial - Event Handling and Callback Functions · Differe…

---

## [Plotting MCMCChains with Makie](https://discourse.julialang.org/t/plotting-mcmcchains-with-makie/84860)

<div class="topic-metadata">

**Author:** [@TheLateKronos](https://discourse.julialang.org/u/TheLateKronos)\
**Replies:** 7\
**Last updated:** [July 29, 2022, 1:13pm UTC](https://discourse.julialang.org/t/plotting-mcmcchains-with-makie/84860 "2022-07-29T13:13:06Z")

</div>

It is currently trivially simple to plot the posterior likelihood functions from a MCMC chain with StatsPlots. However, it was quite easy to make a function that allows the same automatic plotting with Makie: Setupjulia…

---

## [What is the proper way to handle a discontinuity in Differential Equations?](https://discourse.julialang.org/t/what-is-the-proper-way-to-handle-a-discontinuity-in-differential-equations/84963)

<div class="topic-metadata">

**Author:** [@Ge0rges](https://discourse.julialang.org/u/Ge0rges)\
**Replies:** 2\
**Last updated:** [July 29, 2022, 12:24pm UTC](https://discourse.julialang.org/t/what-is-the-proper-way-to-handle-a-discontinuity-in-differential-equations/84963 "2022-07-29T12:24:38Z")

</div>

Hello all, I have a system of two differential equations that I am able to solve using Rosenbrock23(). I now would like to model this system with one modification, where at time t I manually change the value of one of …

[Previous page](https://discourse.julialang.org/c/domain/10.md?page=261)

[Next page](https://discourse.julialang.org/c/domain/10.md?page=263)