# Specific Domains

**URL:** https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=119

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

**Page:** 120

---

## [How to do "reactive conditional plotting" in Makie recipes?](https://discourse.julialang.org/t/how-to-do-reactive-conditional-plotting-in-makie-recipes/112871)

<div class="topic-metadata">

**Author:** [@eliascarv](https://discourse.julialang.org/u/eliascarv)\
**Replies:** 6\
**Last updated:** [April 12, 2024, 8:00pm UTC](https://discourse.julialang.org/t/how-to-do-reactive-conditional-plotting-in-makie-recipes/112871 "2024-04-12T20:00:26Z")

</div>

I want to show a plot in the recipe if any condition is true and do something else otherwise, maybe another plot, or do nothing in the case of my MWE. How can I do this in Makie recipes? MWE: using Makie using GLMakie…

---

## [How to test if a function is continuous](https://discourse.julialang.org/t/how-to-test-if-a-function-is-continuous/112737)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 18\
**Last updated:** [April 12, 2024, 4:54pm UTC](https://discourse.julialang.org/t/how-to-test-if-a-function-is-continuous/112737 "2024-04-12T16:54:10Z")

</div>

I am writing some tests for a package (not yet public) where I want to test that a function f(x::Real)::Real is continuous at some x. Is there a numerical way of doing this? A solution I thought of is finding the left-…

---

## [Computing matrix derivatives for eigenvaule sensitivity analysis](https://discourse.julialang.org/t/computing-matrix-derivatives-for-eigenvaule-sensitivity-analysis/112845)

<div class="topic-metadata">

**Author:** [@jwalker](https://discourse.julialang.org/u/jwalker)\
**Replies:** 6\
**Last updated:** [April 12, 2024, 3:29pm UTC](https://discourse.julialang.org/t/computing-matrix-derivatives-for-eigenvaule-sensitivity-analysis/112845 "2024-04-12T15:29:01Z")

</div>

I have a problem where I would like to apply eigenvalue sensitivity analysis, using Enzyme to compute the derivatives of the matrix A(\\mathbf{p}), where \\mathbf{p} = (p\_1,\\ldots,p\_n). Per @stevengj’s notes, if I want the…

---

## [Testing multimodal distribution](https://discourse.julialang.org/t/testing-multimodal-distribution/44346)

<div class="topic-metadata">

**Author:** [@Dario\_Sarra](https://discourse.julialang.org/u/Dario_Sarra)\
**Replies:** 5\
**Last updated:** [April 12, 2024, 12:21pm UTC](https://discourse.julialang.org/t/testing-multimodal-distribution/44346 "2024-04-12T12:21:42Z")

</div>

I am trying to find implementations of any statistical method to test whether a sample is unimodal or multimodal distributed. From Wikipedia Multimodal Distribution, you can find multiple example of such tests. In parti…

---

## [Colours from a \[0,1\] float?](https://discourse.julialang.org/t/colours-from-a-0-1-float/112866)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 4\
**Last updated:** [April 12, 2024, 11:10am UTC](https://discourse.julialang.org/t/colours-from-a-0-1-float/112866 "2024-04-12T11:10:00Z")

</div>

I have some data that I want to represent in a histogram with a different colour of the bin that depends on a float in the 0,1 range. I can see I can pass to the histogram call a colour parameter with the vector of colo…

---

## [Solving type instability when using Graphs.jl](https://discourse.julialang.org/t/solving-type-instability-when-using-graphs-jl/112862)

<div class="topic-metadata">

**Author:** [@alemelis](https://discourse.julialang.org/u/alemelis)\
**Replies:** 4\
**Last updated:** [April 12, 2024, 10:50am UTC](https://discourse.julialang.org/t/solving-type-instability-when-using-graphs-jl/112862 "2024-04-12T10:50:56Z")

</div>

This is the flame plot I get when profiling my iterative solver. You can see it’s pointed by tiny red blocks, which I’m trying to address one by one. In particular I’m now looking at the one highlighted by the red a…

---

## [Attaching leads to honeycomb lattice](https://discourse.julialang.org/t/attaching-leads-to-honeycomb-lattice/111304)

<div class="topic-metadata">

**Author:** [@alex123](https://discourse.julialang.org/u/alex123)\
**Replies:** 6\
**Last updated:** [April 12, 2024, 9:30am UTC](https://discourse.julialang.org/t/attaching-leads-to-honeycomb-lattice/111304 "2024-04-12T09:30:41Z")

</div>

Hi, I am trying to attach leads with a honeycomb basis to a honeycomb lattice, and can’t seem to understand why my code doesn’t work, I don’t really understand the logic of attaching the lead (and how the r variable is …

---

## [Implementation of algorithms ODE Problems differentialequations.jl](https://discourse.julialang.org/t/implementation-of-algorithms-ode-problems-differentialequations-jl/112840)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 2\
**Last updated:** [April 12, 2024, 5:45am UTC](https://discourse.julialang.org/t/implementation-of-algorithms-ode-problems-differentialequations-jl/112840 "2024-04-12T05:45:01Z")

</div>

Hello! Where i can find implementation of algorithms ODE Problems differentialequations.jl in repository? For example Vern9. Thank you for your help!

---

## [Allocation-free higher-order derivatives of R^n-\>R](https://discourse.julialang.org/t/allocation-free-higher-order-derivatives-of-r-n-r/112829)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 4\
**Last updated:** [April 12, 2024, 3:26am UTC](https://discourse.julialang.org/t/allocation-free-higher-order-derivatives-of-r-n-r/112829 "2024-04-12T03:26:16Z")

</div>

This question describes a particular AD application, I am either looking for a solution, someone telling me why it is not a good idea, or potential collaborators who are interested in the same problem, which comes up fre…

---

## [ExaModels.jl and MadNLP.jl on GPUs](https://discourse.julialang.org/t/examodels-jl-and-madnlp-jl-on-gpus/111377)

<div class="topic-metadata">

**Author:** [@sshin23](https://discourse.julialang.org/u/sshin23)\
**Replies:** 0\
**Last updated:** [March 8, 2024, 8:09pm UTC](https://discourse.julialang.org/t/examodels-jl-and-madnlp-jl-on-gpus/111377 "2024-03-08T20:09:16Z")

</div>

ExaModels and MadNLP: Solving Large-Scale Nonlinear Programs on GPUs We are pleased to announce that ExaModels v0.6 and MadNLP v0.8 now have full support for solving large, sparse, constrained nonlinear programs on NVIDI…

---

## [Calculating duals with SOS2 constraints](https://discourse.julialang.org/t/calculating-duals-with-sos2-constraints/112831)

<div class="topic-metadata">

**Author:** [@RuaridhMacd](https://discourse.julialang.org/u/RuaridhMacd)\
**Replies:** 2\
**Last updated:** [April 12, 2024, 2:24am UTC](https://discourse.julialang.org/t/calculating-duals-with-sos2-constraints/112831 "2024-04-12T02:24:51Z")

</div>

Hi, I am working with a model in JuMP and Gurobi which is usually a LP or MILP. I have added a set of SOS2 constraints to create a piecewise linear function. Gurobi successfully reformulates the constraints to avoid nee…

---

## [Count Binary Variables in JuMP Model](https://discourse.julialang.org/t/count-binary-variables-in-jump-model/87232)

<div class="topic-metadata">

**Author:** [@chelseas](https://discourse.julialang.org/u/chelseas)\
**Replies:** 4\
**Last updated:** [April 12, 2024, 12:30am UTC](https://discourse.julialang.org/t/count-binary-variables-in-jump-model/87232 "2024-04-12T00:30:24Z")

</div>

What is the proper way to count binary variables in my jump model? I have been using: MathOptInterface = MOI const\_types = list\_of\_constraint\_types(model) n\_bin = 0 n\_lin = 0 for i = 1:length(const\_typ…

---

## [Newbie question regarding GPU with metal.jl](https://discourse.julialang.org/t/newbie-question-regarding-gpu-with-metal-jl/112850)

<div class="topic-metadata">

**Author:** [@cloner0110](https://discourse.julialang.org/u/cloner0110)\
**Replies:** 0\
**Last updated:** [April 11, 2024, 10:39pm UTC](https://discourse.julialang.org/t/newbie-question-regarding-gpu-with-metal-jl/112850 "2024-04-11T22:39:14Z")

</div>

Hi Everyone, Im a bit new to julia and GPU programming, Im about to solve a system of equation using metal.jl on a GPU and see the performance, as a starting point, does any of you have any suggestions? Cheers,

---

## [Pre allocating interpolant with Interpolations.jl](https://discourse.julialang.org/t/pre-allocating-interpolant-with-interpolations-jl/112712)

<div class="topic-metadata">

**Author:** [@Iddingsite](https://discourse.julialang.org/u/Iddingsite)\
**Replies:** 5\
**Last updated:** [April 11, 2024, 6:49pm UTC](https://discourse.julialang.org/t/pre-allocating-interpolant-with-interpolations-jl/112712 "2024-04-11T18:49:05Z")

</div>

Hi, I have a loop where I need to create new interpolants each time. I would like to know if I can somehow pre allocate the interpolant in Interpolations.jl so I don’t have to create new memory allocation at each iterat…

---

## [Cairomakie inset plot at specific (x,y) coordinates](https://discourse.julialang.org/t/cairomakie-inset-plot-at-specific-x-y-coordinates/84797)

<div class="topic-metadata">

**Author:** [@Salmon](https://discourse.julialang.org/u/Salmon)\
**Replies:** 3\
**Last updated:** [April 11, 2024, 12:15pm UTC](https://discourse.julialang.org/t/cairomakie-inset-plot-at-specific-x-y-coordinates/84797 "2024-04-11T12:15:54Z")

</div>

Hi Guys, by using BBox() one can create an inset plot. The location is determined by the arguments, which are given as pixels. This works mostly but for certain use cases, it would be nice to specify an absolute locati…

---

## [How to differentiate using Zygote in the presence of missing values](https://discourse.julialang.org/t/how-to-differentiate-using-zygote-in-the-presence-of-missing-values/112822)

<div class="topic-metadata">

**Author:** [@Nikos\_Gianniotis](https://discourse.julialang.org/u/Nikos_Gianniotis)\
**Replies:** 0\
**Last updated:** [April 11, 2024, 10:05am UTC](https://discourse.julialang.org/t/how-to-differentiate-using-zygote-in-the-presence-of-missing-values/112822 "2024-04-11T10:05:18Z")

</div>

I have encountered a problem in my code. While I cannot easily write a MWE for my particular case, I thought that perhaps dealing with the following issue would help me further. Let us say, that we have some data y and…

---

## [RQ decomposition](https://discourse.julialang.org/t/rq-decomposition/112795)

<div class="topic-metadata">

**Author:** [@RJDennis](https://discourse.julialang.org/u/RJDennis)\
**Replies:** 12\
**Last updated:** [April 10, 2024, 11:13pm UTC](https://discourse.julialang.org/t/rq-decomposition/112795 "2024-04-10T23:13:04Z")

</div>

Hi, can anyone point me to a Julia implementation of the RQ decomposition. Thanks.

---

## [DiffEq: Discrete callback causes dt\<=dtmin](https://discourse.julialang.org/t/diffeq-discrete-callback-causes-dt-dtmin/104679)

<div class="topic-metadata">

**Author:** [@aaliqureshi](https://discourse.julialang.org/u/aaliqureshi)\
**Replies:** 6\
**Last updated:** [April 10, 2024, 4:07pm UTC](https://discourse.julialang.org/t/diffeq-discrete-callback-causes-dt-dtmin/104679 "2024-04-10T16:07:50Z")

</div>

Hello! I am simulating a DAE model as an ODEProblem using a constant mass matrix. A model parameter p changes from 0 to 1 at time t=t1 using DiscreteCallback. The solver returns dt\<=dtmin at t=t1 when p changes. Howeve…

---

## [Bushes instead of text](https://discourse.julialang.org/t/bushes-instead-of-text/112663)

<div class="topic-metadata">

**Author:** [@qqs\_qqs](https://discourse.julialang.org/u/qqs_qqs)\
**Replies:** 4\
**Last updated:** [April 10, 2024, 3:54pm UTC](https://discourse.julialang.org/t/bushes-instead-of-text/112663 "2024-04-10T15:54:05Z")

</div>

I have this simple (or so I thought) code in which I wanted to generate random texts in files: z = \[“b”, “5”, “x”\] for k in 1:10 how much = rand(3:50) filename = “job1\_$k.txt” file = open(filename, “in”) for i in 1…

---

## [I cannot use Julia kernel in Jupyter Notebook/Jupyter lab after updating to Julia 1.9.0](https://discourse.julialang.org/t/i-cannot-use-julia-kernel-in-jupyter-notebook-jupyter-lab-after-updating-to-julia-1-9-0/99649)

<div class="topic-metadata">

**Author:** [@arinbasu](https://discourse.julialang.org/u/arinbasu)\
**Replies:** 15\
**Last updated:** [April 10, 2024, 3:46pm UTC](https://discourse.julialang.org/t/i-cannot-use-julia-kernel-in-jupyter-notebook-jupyter-lab-after-updating-to-julia-1-9-0/99649 "2024-04-10T15:46:19Z")

</div>

I have installed Jupyter Lab, and IJulia. Now if I do $ jupyter lab, I get to see the Jupyter Lab web interface, where if I select “Julia-1.9.0 kernel”, the kernel does not load, and in the terminal I see the followin…

---

## [JuMP, large expression, and hessian "limited-memory" option in Ipopt](https://discourse.julialang.org/t/jump-large-expression-and-hessian-limited-memory-option-in-ipopt/112609)

<div class="topic-metadata">

**Author:** [@niro](https://discourse.julialang.org/u/niro)\
**Replies:** 8\
**Last updated:** [April 9, 2024, 10:04pm UTC](https://discourse.julialang.org/t/jump-large-expression-and-hessian-limited-memory-option-in-ipopt/112609 "2024-04-09T22:04:09Z")

</div>

Hi all, I have a question regarding JuMP, which I try to use to solve an optimal control problem. For this problem, I have to generate a quite complicated function using symbolic computations. Then, I build a nonlinear @…

---

## [Custom recorder on stages with different constraints](https://discourse.julialang.org/t/custom-recorder-on-stages-with-different-constraints/112718)

<div class="topic-metadata">

**Author:** [@FelixSchmidt](https://discourse.julialang.org/u/FelixSchmidt)\
**Replies:** 2\
**Last updated:** [April 9, 2024, 9:19pm UTC](https://discourse.julialang.org/t/custom-recorder-on-stages-with-different-constraints/112718 "2024-04-09T21:19:59Z")

</div>

Hi Oscar, I’ve been having issues with creating a custom recorder in SDDP.simulate. I have a capacity expansion problem where stage 1 covers the investment decision and stages 2,…T are monthly dispatch problems at hourl…

---

## [Accessing dual variables](https://discourse.julialang.org/t/accessing-dual-variables/112719)

<div class="topic-metadata">

**Author:** [@heiwie](https://discourse.julialang.org/u/heiwie)\
**Replies:** 2\
**Last updated:** [April 9, 2024, 7:44pm UTC](https://discourse.julialang.org/t/accessing-dual-variables/112719 "2024-04-09T19:44:54Z")

</div>

I am wondering how I can name constraints I define in a for-loop so I later could access the dual variables. My code is the following: for (m, parentnode) in enumerate(lattice\_edges\[2\]) for (n, childnode) in enumera…

---

## [Generalized eigenvalue problem in ApproxFun.jl](https://discourse.julialang.org/t/generalized-eigenvalue-problem-in-approxfun-jl/112441)

<div class="topic-metadata">

**Author:** [@jwalker](https://discourse.julialang.org/u/jwalker)\
**Replies:** 2\
**Last updated:** [April 9, 2024, 7:26pm UTC](https://discourse.julialang.org/t/generalized-eigenvalue-problem-in-approxfun-jl/112441 "2024-04-09T19:26:02Z")

</div>

Looking through the examples in ApproxFun.jl, I didn’t see one of a generalized eigenvalue problem. Is this possible? Is there a short example?

---

## [Root finding / NonLinear problem](https://discourse.julialang.org/t/root-finding-nonlinear-problem/112720)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 12\
**Last updated:** [April 9, 2024, 6:18pm UTC](https://discourse.julialang.org/t/root-finding-nonlinear-problem/112720 "2024-04-09T18:18:38Z")

</div>

Question: I am trying to solve a root-finding problem on a vector of 14 elements, but each function evaluation needs about 7 seconds. What is the best solver for this? And is there a way to declare a valid range of th…

---

## [Can I define a type-unstable EnzymeRule?](https://discourse.julialang.org/t/can-i-define-a-type-unstable-enzymerule/112732)

<div class="topic-metadata">

**Author:** [@gdalle](https://discourse.julialang.org/u/gdalle)\
**Replies:** 1\
**Last updated:** [April 9, 2024, 6:09pm UTC](https://discourse.julialang.org/t/can-i-define-a-type-unstable-enzymerule/112732 "2024-04-09T18:09:34Z")

</div>

I am trying to define custom Enzyme rules for my package ImplicitDifferentiation.jl, but I’m held back by type inference failures inside my rules. The typical error message looks like Enzyme execution failed. Enzyme: …

---

## [Problem with observables and animations](https://discourse.julialang.org/t/problem-with-observables-and-animations/112731)

<div class="topic-metadata">

**Author:** [@oddsen](https://discourse.julialang.org/u/oddsen)\
**Replies:** 2\
**Last updated:** [April 9, 2024, 2:45pm UTC](https://discourse.julialang.org/t/problem-with-observables-and-animations/112731 "2024-04-09T14:45:45Z")

</div>

Hi, Im trying to animate something with points In this code, it is supposed to show points moving for each iteration. The problem is that the points from the previous iteration aren’t being removed, resulting in an a…

---

## [Matrix exponential \* vector for various scales](https://discourse.julialang.org/t/matrix-exponential-vector-for-various-scales/112022)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 38\
**Last updated:** [April 9, 2024, 2:43pm UTC](https://discourse.julialang.org/t/matrix-exponential-vector-for-various-scales/112022 "2024-04-09T14:43:13Z")

</div>

Given a matrix Q which is a continuous-time transition matrix, ie Q\_{i,j} \\ge 0 for j \\ne i, \\sum\_j Q\_{i,j} = 0 \\ \\forall i I would like to calculate \\exp(Qt)v for multiple t's (fixed v, Q). Potentially Q is around 5…

---

## [Ordinal probit model with Turing - difficulty sampling](https://discourse.julialang.org/t/ordinal-probit-model-with-turing-difficulty-sampling/112528)

<div class="topic-metadata">

**Author:** [@opera\_malenky](https://discourse.julialang.org/u/opera_malenky)\
**Replies:** 4\
**Last updated:** [April 9, 2024, 1:43pm UTC](https://discourse.julialang.org/t/ordinal-probit-model-with-turing-difficulty-sampling/112528 "2024-04-09T13:43:58Z")

</div>

I am trying to use Turing to estimate an ordinal multinomial model with partially fixed thresholds/cutpoints, but I’m wondering if I’m doing something wrong: metropolis-hastings sampling works great, but the autodiff-bas…

---

## [Problem with JuMP + Gurobi.jl](https://discourse.julialang.org/t/problem-with-jump-gurobi-jl/111505)

<div class="topic-metadata">

**Author:** [@njozefow](https://discourse.julialang.org/u/njozefow)\
**Replies:** 7\
**Last updated:** [April 9, 2024, 1:36pm UTC](https://discourse.julialang.org/t/problem-with-jump-gurobi-jl/111505 "2024-04-09T13:36:03Z")

</div>

Hello, In continuation of the thread : Poor time performance in modifying parameters of JuMP models when using Gurobi - #10 by odow Here is a code that show the behavior (it emulates the behavior of a column generation…

[Previous page](https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=118)

[Next page](https://discourse.julialang.org/c/domain/10.md?no_subcategories=false&page=120)