# Specific Domains

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

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

**Page:** 33

---

## [Numerical artifact in solution of system of ODEs](https://discourse.julialang.org/t/numerical-artifact-in-solution-of-system-of-odes/130941)

<div class="topic-metadata">

**Author:** [@Gako](https://discourse.julialang.org/u/Gako)\
**Replies:** 12\
**Last updated:** [July 24, 2025, 10:09am UTC](https://discourse.julialang.org/t/numerical-artifact-in-solution-of-system-of-odes/130941 "2025-07-24T10:09:46Z")

</div>

I am currently solving the following system of ODEs: \\begin{cases} \\ddot\\phi=-3H\\dot{\\phi}-\\phi\\\\ \\dot{a}=Ha\\\\ \\ddot\\sigma=-H\\dot{\\sigma}-\\frac{m\_{Eff}^2}{a^2}\\sigma\\\\ \\ddot{\\rho}\_0=\\frac{4\\rho\_2}{a^2}-H\\dot{\\rho}\_{0}-…

---

## [Ideal Diode and Ideal Switch instability problem](https://discourse.julialang.org/t/ideal-diode-and-ideal-switch-instability-problem/107769)

<div class="topic-metadata">

**Author:** [@luca.andriotto](https://discourse.julialang.org/u/luca.andriotto)\
**Replies:** 11\
**Last updated:** [July 24, 2025, 9:44am UTC](https://discourse.julialang.org/t/ideal-diode-and-ideal-switch-instability-problem/107769 "2025-07-24T09:44:11Z")

</div>

Greetings, I am new to Julia community and, in particular, to the use of ModelingToolkit package. I would like to implement an IdealDiode component to be replaced by a short circuit when the controlled variable (i.e. i…

---

## [Second-order cone constraint in JuMP](https://discourse.julialang.org/t/second-order-cone-constraint-in-jump/130927)

<div class="topic-metadata">

**Author:** [@Z\_M](https://discourse.julialang.org/u/Z_M)\
**Replies:** 9\
**Last updated:** [July 24, 2025, 3:14am UTC](https://discourse.julialang.org/t/second-order-cone-constraint-in-jump/130927 "2025-07-24T03:14:58Z")

</div>

Hi guys! How to solve this constraint in Julia in JuMP, where the Q is a 2×2 Hessian matrix of a convex function, so it is a positive semidefinite matrix. I know it can be reformulated as a second-order cone constraint…

---

## [Sparsity of SimpleChains Errors](https://discourse.julialang.org/t/sparsity-of-simplechains-errors/130947)

<div class="topic-metadata">

**Author:** [@johnb](https://discourse.julialang.org/u/johnb)\
**Replies:** 6\
**Last updated:** [July 23, 2025, 8:31pm UTC](https://discourse.julialang.org/t/sparsity-of-simplechains-errors/130947 "2025-07-23T20:31:58Z")

</div>

I’ve been playing around with SimpleChains.jl for potential use in a larger sparse model, but I’m having some issues figuring out how to obtain a sparsity pattern. After some testing with both SparseConnectivityTracer.jl…

---

## [Set distribution for function parameters](https://discourse.julialang.org/t/set-distribution-for-function-parameters/130963)

<div class="topic-metadata">

**Author:** [@Cyan](https://discourse.julialang.org/u/Cyan)\
**Replies:** 2\
**Last updated:** [July 23, 2025, 2:18pm UTC](https://discourse.julialang.org/t/set-distribution-for-function-parameters/130963 "2025-07-23T14:18:24Z")

</div>

Hi, all. I am stucked by the problems of building a @model for a parameters of a function when using Turing. Such as for function f(x; b, c), I want to obtain a model like @model fuction(x,y) b ~ truncated(Normal(1., 1)…

---

## [An infeasibility that seems to be dependent on the state of HiGHS (or JuMP)](https://discourse.julialang.org/t/an-infeasibility-that-seems-to-be-dependent-on-the-state-of-highs-or-jump/130969)

<div class="topic-metadata">

**Author:** [@BdeKoning](https://discourse.julialang.org/u/BdeKoning)\
**Replies:** 5\
**Last updated:** [July 23, 2025, 12:52pm UTC](https://discourse.julialang.org/t/an-infeasibility-that-seems-to-be-dependent-on-the-state-of-highs-or-jump/130969 "2025-07-23T12:52:12Z")

</div>

Hi all, I’m working on a water allocation optimization code with JuMP + HiGHS, and I came across an infeasibility. The weird thing is: HiGHS concludes that my model is infeasible when it runs embedded in my code, but whe…

---

## [Annotation (annotation!() and text!()) not working as expected for longitude and latitude input in GeoMakie geoaxis with Tyler map tiles](https://discourse.julialang.org/t/annotation-annotation-and-text-not-working-as-expected-for-longitude-and-latitude-input-in-geomakie-geoaxis-with-tyler-map-tiles/130971)

<div class="topic-metadata">

**Author:** [@ewillcox](https://discourse.julialang.org/u/ewillcox)\
**Replies:** 3\
**Last updated:** [July 23, 2025, 11:02am UTC](https://discourse.julialang.org/t/annotation-annotation-and-text-not-working-as-expected-for-longitude-and-latitude-input-in-geomakie-geoaxis-with-tyler-map-tiles/130971 "2025-07-23T11:02:02Z")

</div>

Problem description I am attempting to annotate some points in a scatter plot in a GeoMakie.jl geoaxis that contains a Tyler.jl maptile. None of my efforts so far turned out as expected. I suspect I need to find the appr…

---

## [Window title for Plots.jl](https://discourse.julialang.org/t/window-title-for-plots-jl/130929)

<div class="topic-metadata">

**Author:** [@GHTaarn](https://discourse.julialang.org/u/GHTaarn)\
**Replies:** 2\
**Last updated:** [July 23, 2025, 9:27am UTC](https://discourse.julialang.org/t/window-title-for-plots-jl/130929 "2025-07-23T09:27:25Z")

</div>

By default, the window title of my plot seems to be “GKS QtTerm”, is there a way to set this to something that I choose? I tried using Plots plot(; window\_title="Abc") but this does not change the anything that I can s…

---

## [Rigorous establishment of the feasibility of the trained policies in SDDP.jl using constraint violation penalties in the objective function](https://discourse.julialang.org/t/rigorous-establishment-of-the-feasibility-of-the-trained-policies-in-sddp-jl-using-constraint-violation-penalties-in-the-objective-function/130961)

<div class="topic-metadata">

**Author:** [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)\
**Replies:** 10\
**Last updated:** [July 23, 2025, 8:32am UTC](https://discourse.julialang.org/t/rigorous-establishment-of-the-feasibility-of-the-trained-policies-in-sddp-jl-using-constraint-violation-penalties-in-the-objective-function/130961 "2025-07-23T08:32:16Z")

</div>

I am working with a multi-stage mixed-integer stochastic global supply chain optimization problem. In this regard, under realistic parameter settings, some constraints are making the instances violate the relatively com…

---

## [3D volume plot for spherical coordinate grid points](https://discourse.julialang.org/t/3d-volume-plot-for-spherical-coordinate-grid-points/129668)

<div class="topic-metadata">

**Author:** [@raman\_kumar](https://discourse.julialang.org/u/raman_kumar)\
**Replies:** 31\
**Last updated:** [July 7, 2025, 12:50pm UTC](https://discourse.julialang.org/t/3d-volume-plot-for-spherical-coordinate-grid-points/129668 "2025-07-07T12:50:47Z")

</div>

How to get volume plot for the following spherical coordinate positions? :child: I have a lot of blocks each having such r, θ and ϕ spherical coordinate position vectors and i want to create 3D volume plot for each bloc…

---

## [DiffEqGPU Trajectory Failure Handling and Heterogeneous Trajectories](https://discourse.julialang.org/t/diffeqgpu-trajectory-failure-handling-and-heterogeneous-trajectories/129962)

<div class="topic-metadata">

**Author:** [@tosm01](https://discourse.julialang.org/u/tosm01)\
**Replies:** 4\
**Last updated:** [July 22, 2025, 12:06pm UTC](https://discourse.julialang.org/t/diffeqgpu-trajectory-failure-handling-and-heterogeneous-trajectories/129962 "2025-07-22T12:06:50Z")

</div>

Hi, I have a general question about the suitability of using GPUs for the kind of problem I’m working on. So far, I’ve had good success using DiffEqGPU for parameter sweeps, specifically in cases where the parameters d…

---

## [Getting results until the point of nonconvergence in DAE system](https://discourse.julialang.org/t/getting-results-until-the-point-of-nonconvergence-in-dae-system/130897)

<div class="topic-metadata">

**Author:** [@total](https://discourse.julialang.org/u/total)\
**Replies:** 7\
**Last updated:** [July 22, 2025, 10:21am UTC](https://discourse.julialang.org/t/getting-results-until-the-point-of-nonconvergence-in-dae-system/130897 "2025-07-22T10:21:12Z")

</div>

Hello, I’m trying to simulate a DAE system with Sundials IDA solver. I implemented my controlvector with callbacks and it works just fine. However, depending on the controls I run into infeasible regions and the solver…

---

## [Symbolics.jl - Retrieving Differential Variables](https://discourse.julialang.org/t/symbolics-jl-retrieving-differential-variables/130891)

<div class="topic-metadata">

**Author:** [@Gray\_GALLANT](https://discourse.julialang.org/u/Gray_GALLANT)\
**Replies:** 2\
**Last updated:** [July 21, 2025, 8:18pm UTC](https://discourse.julialang.org/t/symbolics-jl-retrieving-differential-variables/130891 "2025-07-21T20:18:07Z")

</div>

Is there a way to retrieve differential variables from an expression using Symbolics.jl? For example, if I use get\_variables on Differential(t)(x(t)) it only returns x(t), but I need to get the differential part as well …

---

## [KKT condition not satisfied in DCOPF](https://discourse.julialang.org/t/kkt-condition-not-satisfied-in-dcopf/130902)

<div class="topic-metadata">

**Author:** [@simoni0](https://discourse.julialang.org/u/simoni0)\
**Replies:** 4\
**Last updated:** [July 21, 2025, 6:33pm UTC](https://discourse.julialang.org/t/kkt-condition-not-satisfied-in-dcopf/130902 "2025-07-21T18:33:41Z")

</div>

I wanted to evaluate whether the KKT condition \\lambda\_i - \\lambda\_j = \\mu^{+}\_{ij} - \\mu^{-}\_{ij} holds for standard DC-OPF. I made a small test on the 5-bus case. using PowerModels using GLPK result = PowerModels.sol…

---

## [How to get rid of noise in phase isosurface plots? Very short MWE included](https://discourse.julialang.org/t/how-to-get-rid-of-noise-in-phase-isosurface-plots-very-short-mwe-included/130890)

<div class="topic-metadata">

**Author:** [@Deceneu](https://discourse.julialang.org/u/Deceneu)\
**Replies:** 5\
**Last updated:** [July 21, 2025, 4:24pm UTC](https://discourse.julialang.org/t/how-to-get-rid-of-noise-in-phase-isosurface-plots-very-short-mwe-included/130890 "2025-07-21T16:24:06Z")

</div>

To keep it brief, I want to plot complex phase (angle) of a complex field. It seemed like Makie.jl’s contour or volume, should do the trick, but I get very noisy results for some reason? using GLMakie ph = \[(r = hypot(…

---

## [Control systems, pid, and high-order transfer functions](https://discourse.julialang.org/t/control-systems-pid-and-high-order-transfer-functions/130837)

<div class="topic-metadata">

**Author:** [@romain-chiap](https://discourse.julialang.org/u/romain-chiap)\
**Replies:** 7\
**Last updated:** [July 21, 2025, 10:26am UTC](https://discourse.julialang.org/t/control-systems-pid-and-high-order-transfer-functions/130837 "2025-07-21T10:26:06Z")

</div>

Hello, It’s specific to control systems and this is my first post, so I still decided to post under “new to Julia” as I might be asking beginner questions. I’m making a 1-dimension model of the PX4 cascaded control loo…

---

## [Set default Makie theme](https://discourse.julialang.org/t/set-default-makie-theme/130893)

<div class="topic-metadata">

**Author:** [@jar1](https://discourse.julialang.org/u/jar1)\
**Replies:** 2\
**Last updated:** [July 21, 2025, 7:09am UTC](https://discourse.julialang.org/t/set-default-makie-theme/130893 "2025-07-21T07:09:09Z")

</div>

I would like to set a default Makie theme for all my plots, but not actually load Makie in every session (because that takes a while). Just whenever I do load Makie, I want it to use my theme. Is there a way to do that?

---

## [Very high memory use with DelayDifferentialEquations](https://discourse.julialang.org/t/very-high-memory-use-with-delaydifferentialequations/99862)

<div class="topic-metadata">

**Author:** [@jarias9](https://discourse.julialang.org/u/jarias9)\
**Replies:** 12\
**Last updated:** [July 21, 2025, 6:05am UTC](https://discourse.julialang.org/t/very-high-memory-use-with-delaydifferentialequations/99862 "2025-07-21T06:05:22Z")

</div>

I have a large PDE system that I originally was solving with an explicit ODE solver via ODEProblem with DiffEq.jl. I have added a component to that system that requires the time derivative info so I changed it to a DDE a…

---

## [Makie max pan zoom](https://discourse.julialang.org/t/makie-max-pan-zoom/114105)

<div class="topic-metadata">

**Author:** [@bjarthur](https://discourse.julialang.org/u/bjarthur)\
**Replies:** 3\
**Last updated:** [July 20, 2025, 6:31pm UTC](https://discourse.julialang.org/t/makie-max-pan-zoom/114105 "2025-07-20T18:31:13Z")

</div>

i’d like to put a maximum on how much a user can pan or zoom out on a plot. but i could not find a function in the makie docs to set an upper bound on the axis limits. so i hacked it like this: using GLMakie f,a,s = s…

---

## [Unhelpful interpolation for solution from differential equation](https://discourse.julialang.org/t/unhelpful-interpolation-for-solution-from-differential-equation/130857)

<div class="topic-metadata">

**Author:** [@miles-1](https://discourse.julialang.org/u/miles-1)\
**Replies:** 3\
**Last updated:** [July 20, 2025, 1:40am UTC](https://discourse.julialang.org/t/unhelpful-interpolation-for-solution-from-differential-equation/130857 "2025-07-20T01:40:23Z")

</div>

I am trying to use the solution from my differential equation to find local maxima. I am following the example found here, but am running into an issue. I have found that the interpolation of my solution is very poor. He…

---

## [Does AMDGPU.jl support integrated graphics?](https://discourse.julialang.org/t/does-amdgpu-jl-support-integrated-graphics/130675)

<div class="topic-metadata">

**Author:** [@Mend-Amar\_Badral](https://discourse.julialang.org/u/Mend-Amar_Badral)\
**Replies:** 3\
**Last updated:** [July 19, 2025, 10:21pm UTC](https://discourse.julialang.org/t/does-amdgpu-jl-support-integrated-graphics/130675 "2025-07-19T22:21:44Z")

</div>

I have ubuntu 24.04 and my gpu1 is nvidia gpu and gpu0 is AMD ATI Radeon RX Vega 6. I was wondering if it is possible to use AMDGPU.jl to use integrated gpu.

---

## [Kernel with dynamic parallelism seems to be calling CPU functions](https://discourse.julialang.org/t/kernel-with-dynamic-parallelism-seems-to-be-calling-cpu-functions/130829)

<div class="topic-metadata">

**Author:** [@victor\_vhrn](https://discourse.julialang.org/u/victor_vhrn)\
**Replies:** 4\
**Last updated:** [July 19, 2025, 6:19pm UTC](https://discourse.julialang.org/t/kernel-with-dynamic-parallelism-seems-to-be-calling-cpu-functions/130829 "2025-07-19T18:19:25Z")

</div>

I’m trying to implement the adjoint\_v1\_custom\_kernel! using dynamic parallelism (See code below). I’m not sure why, but when trying to run it I get A LOT of unsupported call errors, after a dwarf warning. I tried to fol…

---

## [User defined ELBO in AdvancedVI or Turing](https://discourse.julialang.org/t/user-defined-elbo-in-advancedvi-or-turing/130833)

<div class="topic-metadata">

**Author:** [@Garren\_Hermanus](https://discourse.julialang.org/u/Garren_Hermanus)\
**Replies:** 1\
**Last updated:** [July 19, 2025, 4:40pm UTC](https://discourse.julialang.org/t/user-defined-elbo-in-advancedvi-or-turing/130833 "2025-07-19T16:40:57Z")

</div>

Hi all, I am fairly new to ADVI and Turing (and sub-packages in general). I get the general idea of ADVI in that we make a set of parameter transformations to obtain an estimator to the intractable expectation \\mathbb{…

---

## [GLMakie -- text versus annotation -- annotation not working?](https://discourse.julialang.org/t/glmakie-text-versus-annotation-annotation-not-working/130854)

<div class="topic-metadata">

**Author:** [@diadora77](https://discourse.julialang.org/u/diadora77)\
**Replies:** 5\
**Last updated:** [July 19, 2025, 11:55am UTC](https://discourse.julialang.org/t/glmakie-text-versus-annotation-annotation-not-working/130854 "2025-07-19T11:55:45Z")

</div>

Hi- Tried using annotation according to the docs and receiving a recipe error. A minimum viable example is provide below which is very similar to the docs example. ts = rand(10,2) lbs = string.(1:10) f = Figure() ax…

---

## [Active Bruno.jl maintainers](https://discourse.julialang.org/t/active-bruno-jl-maintainers/130864)

<div class="topic-metadata">

**Author:** [@BioinfoBuggaBoo](https://discourse.julialang.org/u/BioinfoBuggaBoo)\
**Replies:** 0\
**Last updated:** [July 19, 2025, 8:37am UTC](https://discourse.julialang.org/t/active-bruno-jl-maintainers/130864 "2025-07-19T08:37:10Z")

</div>

Just posting here to see if any maintainers of Bruno.jl are active here? I’m using the package regularly at present and am keen to discuss Issue 45 and possibly any work they may have been half way through completing be…

---

## [Why are the entries identical in \`zeros(JuMP.AffExpr, N)\`](https://discourse.julialang.org/t/why-are-the-entries-identical-in-zeros-jump-affexpr-n/130860)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 9\
**Last updated:** [July 19, 2025, 8:04am UTC](https://discourse.julialang.org/t/why-are-the-entries-identical-in-zeros-jump-affexpr-n/130860 "2025-07-19T08:04:07Z")

</div>

My understanding about the zeros function in julia is that it has independent entries. But today I find a weird behavior which I cannot explain why. Compare the following unexpected outcome 🔴 and the expected outcome ✅ …

---

## [Getting warning, but model works well](https://discourse.julialang.org/t/getting-warning-but-model-works-well/130850)

<div class="topic-metadata">

**Author:** [@Bart\_van\_de\_Lint](https://discourse.julialang.org/u/Bart_van_de_Lint)\
**Replies:** 5\
**Last updated:** [July 19, 2025, 2:11am UTC](https://discourse.julialang.org/t/getting-warning-but-model-works-well/130850 "2025-07-19T02:11:41Z")

</div>

I am getting this warning during mtkcompile, while I was not getting it before. It might be related to an update in StructuralTransformations. The model still works fine. How can I disable the warning? ┌ Warning: Intern…

---

## [Energy conserving (or symplectic) solvers for differential equations?](https://discourse.julialang.org/t/energy-conserving-or-symplectic-solvers-for-differential-equations/130840)

<div class="topic-metadata">

**Author:** [@Flow](https://discourse.julialang.org/u/Flow)\
**Replies:** 2\
**Last updated:** [July 18, 2025, 10:52am UTC](https://discourse.julialang.org/t/energy-conserving-or-symplectic-solvers-for-differential-equations/130840 "2025-07-18T10:52:41Z")

</div>

Hello there, I am relatively new to Julia and am currently trying to calculate trajectories of a physical system on GPUs, which in generally works like a charm, as long as I do not include an minimum coupling term in my…

---

## [JuMP : Vector version of set\_lower\_bound and set\_upper\_bound?](https://discourse.julialang.org/t/jump-vector-version-of-set-lower-bound-and-set-upper-bound/130703)

<div class="topic-metadata">

**Author:** [@LaurentPlagne](https://discourse.julialang.org/u/LaurentPlagne)\
**Replies:** 13\
**Last updated:** [July 18, 2025, 10:05am UTC](https://discourse.julialang.org/t/jump-vector-version-of-set-lower-bound-and-set-upper-bound/130703 "2025-07-18T10:05:14Z")

</div>

Hi, I work on a project involving the resolution of numerous (10K) linear optimization problems. JuMP is very expressive and convenient allowing to work with different backend solvers : many thanks to the JuMP devs for…

---

## [How to define a variable with specific values from a set using JuMP](https://discourse.julialang.org/t/how-to-define-a-variable-with-specific-values-from-a-set-using-jump/130808)

<div class="topic-metadata">

**Author:** [@Luiz\_Monteiro](https://discourse.julialang.org/u/Luiz_Monteiro)\
**Replies:** 3\
**Last updated:** [July 18, 2025, 1:21am UTC](https://discourse.julialang.org/t/how-to-define-a-variable-with-specific-values-from-a-set-using-jump/130808 "2025-07-18T01:21:56Z")

</div>

Hi, I’m using JuMP and I would like to define a variable that can only assume specific values from a predefined set, for example: allowed\_values = \[10.0, 10.5, 11.0, 11.5, 12.0\] Is there a way to define a variable in …

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

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