# Specific Domains

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

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

**Page:** 37

---

## [SDDP model sometimes fails during training, or simulation](https://discourse.julialang.org/t/sddp-model-sometimes-fails-during-training-or-simulation/130139)

<div class="topic-metadata">

**Author:** [@slwu89](https://discourse.julialang.org/u/slwu89)\
**Replies:** 2\
**Last updated:** [June 23, 2025, 9:49pm UTC](https://discourse.julialang.org/t/sddp-model-sometimes-fails-during-training-or-simulation/130139 "2025-06-23T21:49:33Z")

</div>

Hi all. I’m seeking help to diagnose a simple inventory problem over a finite horizon I’ve set up with SDDP. Currently, my model is such that the stochasticity is demand, for which we assume we know the intensity of the …

---

## [ExponentialUtilities.jl expv!() allocations](https://discourse.julialang.org/t/exponentialutilities-jl-expv-allocations/129917)

<div class="topic-metadata">

**Author:** [@random\_guy](https://discourse.julialang.org/u/random_guy)\
**Replies:** 11\
**Last updated:** [June 23, 2025, 4:56pm UTC](https://discourse.julialang.org/t/exponentialutilities-jl-expv-allocations/129917 "2025-06-23T16:56:05Z")

</div>

Hi everyone, I’m trying to use Julia’s ExponentialUtilities.jl package to perform matrix-exponential-matrix multiplication. I have a sparse NxN Hamiltonian H, a small time step dt, and an NxN matrix A. The goal is to co…

---

## [Nested loop in @mtkmodel / Matrix manipulation](https://discourse.julialang.org/t/nested-loop-in-mtkmodel-matrix-manipulation/130133)

<div class="topic-metadata">

**Author:** [@madtorresga](https://discourse.julialang.org/u/madtorresga)\
**Replies:** 2\
**Last updated:** [June 23, 2025, 4:56pm UTC](https://discourse.julialang.org/t/nested-loop-in-mtkmodel-matrix-manipulation/130133 "2025-06-23T16:56:09Z")

</div>

Hi :slight\_smile: I’m trying to automate the battery cells’ creation from an already defined ECM model: const Qs = \[4.8 4.4 4.2 ; 3.7 3.5 3.3\] @mtkmodel Series\_Connection begin @variables begin v\_max …

---

## [How to make a Proximal Point using Julia](https://discourse.julialang.org/t/how-to-make-a-proximal-point-using-julia/130132)

<div class="topic-metadata">

**Author:** [@Deyvid\_Andrade](https://discourse.julialang.org/u/Deyvid_Andrade)\
**Replies:** 2\
**Last updated:** [June 23, 2025, 2:28pm UTC](https://discourse.julialang.org/t/how-to-make-a-proximal-point-using-julia/130132 "2025-06-23T14:28:00Z")

</div>

I am new to julia and I’d like create an iterative function that can to modeling the proximal point method. And, apply to solve an equilibrium problems, i.e., Given function f(x,y) convex only in x, the algorithm is def…

---

## [Modeling piecewise linear diodes (discontinuous behaviors) with Dyad?](https://discourse.julialang.org/t/modeling-piecewise-linear-diodes-discontinuous-behaviors-with-dyad/130058)

<div class="topic-metadata">

**Author:** [@pierre-haessig](https://discourse.julialang.org/u/pierre-haessig)\
**Replies:** 5\
**Last updated:** [June 23, 2025, 9:45am UTC](https://discourse.julialang.org/t/modeling-piecewise-linear-diodes-discontinuous-behaviors-with-dyad/130058 "2025-06-23T09:45:47Z")

</div>

Hello, In the light of Dyad announcement, and perhaps as a continuation of the Ideal Diode and Ideal Switch instability problem or Diode model not working as expected with ModelingToolkit discussions, I’d like to know h…

---

## [Sensitivity of an ODEProblem defined by another package (QuantumToolbox)](https://discourse.julialang.org/t/sensitivity-of-an-odeproblem-defined-by-another-package-quantumtoolbox/129654)

<div class="topic-metadata">

**Author:** [@Jeff\_Wack](https://discourse.julialang.org/u/Jeff_Wack)\
**Replies:** 6\
**Last updated:** [June 23, 2025, 9:08am UTC](https://discourse.julialang.org/t/sensitivity-of-an-odeproblem-defined-by-another-package-quantumtoolbox/129654 "2025-06-23T09:08:12Z")

</div>

I am interested in calculating derivatives of a quantum state defined by QuantumToolbox with respect to system parameters in order to calculate quantum fisher information. I have found initial success with a finite diffe…

---

## [Piecewise linear for a supply function](https://discourse.julialang.org/t/piecewise-linear-for-a-supply-function/130121)

<div class="topic-metadata">

**Author:** [@nemoincognito](https://discourse.julialang.org/u/nemoincognito)\
**Replies:** 1\
**Last updated:** [June 23, 2025, 8:08am UTC](https://discourse.julialang.org/t/piecewise-linear-for-a-supply-function/130121 "2025-06-23T08:08:52Z")

</div>

I am modelling costs for mines that have a fixed cost and a variable cost to slice the mine into realistic average costs for various levels of output. mine\_name = prod\_capa\_cost\_df\[100,:mine\_name\] capacity = prod\_capa\_c…

---

## [Multicommodity Armington style adjustments](https://discourse.julialang.org/t/multicommodity-armington-style-adjustments/130099)

<div class="topic-metadata">

**Author:** [@nemoincognito](https://discourse.julialang.org/u/nemoincognito)\
**Replies:** 3\
**Last updated:** [June 23, 2025, 7:38am UTC](https://discourse.julialang.org/t/multicommodity-armington-style-adjustments/130099 "2025-06-23T07:38:10Z")

</div>

After some generous help on this forum we have solved a very chunky multicommodity network flow problem. The question we have now is that there is some discrepancy with our “backcast” using optimal flow and what is recor…

---

## [Catalyst.jl - what is the appropriate way to model energy?](https://discourse.julialang.org/t/catalyst-jl-what-is-the-appropriate-way-to-model-energy/130088)

<div class="topic-metadata">

**Author:** [@spackman](https://discourse.julialang.org/u/spackman)\
**Replies:** 6\
**Last updated:** [June 23, 2025, 1:33am UTC](https://discourse.julialang.org/t/catalyst-jl-what-is-the-appropriate-way-to-model-energy/130088 "2025-06-23T01:33:40Z")

</div>

Hi there! I’m interested in tracking the thermodynamics of a chemical system along with the kinetics, but I don’t see any example documentation that would indicate the best/native way to do this. I think in an older ve…

---

## [Not achieving convergence after prolonged training](https://discourse.julialang.org/t/not-achieving-convergence-after-prolonged-training/130093)

<div class="topic-metadata">

**Author:** [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)\
**Replies:** 3\
**Last updated:** [June 22, 2025, 8:03pm UTC](https://discourse.julialang.org/t/not-achieving-convergence-after-prolonged-training/130093 "2025-06-22T20:03:21Z")

</div>

I am dealing with an infinite-horizon multi-stage mixed-integer linear stochastic global supply chain optimization problem and using SDDP.jl for its solution. In this regard, after a significant amount of training, the …

---

## [Issue when using previous solution for initial guess with MIRK Boundary Value Problem solver](https://discourse.julialang.org/t/issue-when-using-previous-solution-for-initial-guess-with-mirk-boundary-value-problem-solver/129948)

<div class="topic-metadata">

**Author:** [@ThermionicsGuy](https://discourse.julialang.org/u/ThermionicsGuy)\
**Replies:** 3\
**Last updated:** [June 21, 2025, 7:26pm UTC](https://discourse.julialang.org/t/issue-when-using-previous-solution-for-initial-guess-with-mirk-boundary-value-problem-solver/129948 "2025-06-21T19:26:57Z")

</div>

I’ve found an odd issue while using a previous solution as the initial guess in solving a boundary value problem using the SciML MIRK solvers and TwoPointBVProblem (part of the BoundaryValueDiffEq.jl package. According t…

---

## [Build a JuMP bilinear expression with a function](https://discourse.julialang.org/t/build-a-jump-bilinear-expression-with-a-function/130091)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 2\
**Last updated:** [June 21, 2025, 2:20pm UTC](https://discourse.julialang.org/t/build-a-jump-bilinear-expression-with-a-function/130091 "2025-06-21T14:20:35Z")

</div>

I’m writing some algorithms recently and I find this bilinear structure is particularly useful, because its construction steps can be reused in the 2 models. Is there any advisable idea? out = JuMP.Model() # outer mast…

---

## [Minimum and maximum for OnlineStats.OrderStats](https://discourse.julialang.org/t/minimum-and-maximum-for-onlinestats-orderstats/129989)

<div class="topic-metadata">

**Author:** [@mkitti](https://discourse.julialang.org/u/mkitti)\
**Replies:** 4\
**Last updated:** [June 20, 2025, 11:00pm UTC](https://discourse.julialang.org/t/minimum-and-maximum-for-onlinestats-orderstats/129989 "2025-06-20T23:00:30Z")

</div>

Would these definitions make sense for OnlineStats.OrderStats from OnlineStats.jl? Base.minimum(o::OrderStats) = minimum(o.ex) Base.maximum(o::OrderStats) = maximum(o.ex) Base.extrema(o::OrderStats) = extrema(o.ex)

---

## [Building Catalyst.jl reaction network from a .csv file - error in construction?](https://discourse.julialang.org/t/building-catalyst-jl-reaction-network-from-a-csv-file-error-in-construction/130043)

<div class="topic-metadata">

**Author:** [@spackman](https://discourse.julialang.org/u/spackman)\
**Replies:** 6\
**Last updated:** [June 20, 2025, 8:51pm UTC](https://discourse.julialang.org/t/building-catalyst-jl-reaction-network-from-a-csv-file-error-in-construction/130043 "2025-06-20T20:51:51Z")

</div>

Hey Julia crew! I am trying to import a reaction network from a .csv file to use with Catalyst.jl via Reaction NetworkImporters.jl. I have been able to successfully construct a ReactionSystem object which can print its…

---

## [Column generation requiring a lot of memory?](https://discourse.julialang.org/t/column-generation-requiring-a-lot-of-memory/130056)

<div class="topic-metadata">

**Author:** [@mposs](https://discourse.julialang.org/u/mposs)\
**Replies:** 4\
**Last updated:** [June 20, 2025, 12:08pm UTC](https://discourse.julialang.org/t/column-generation-requiring-a-lot-of-memory/130056 "2025-06-20T12:08:11Z")

</div>

Hi I’m solving a (rather large) LP through two approaches. One involves solving the full LP. The second consists in doing column generation, gradually adding variables through @variable and set\_normalized\_coefficient. I…

---

## [Continuous event detection with @mtkmodel - multiple events](https://discourse.julialang.org/t/continuous-event-detection-with-mtkmodel-multiple-events/130037)

<div class="topic-metadata">

**Author:** [@RThul](https://discourse.julialang.org/u/RThul)\
**Replies:** 2\
**Last updated:** [June 20, 2025, 8:11am UTC](https://discourse.julialang.org/t/continuous-event-detection-with-mtkmodel-multiple-events/130037 "2025-06-20T08:11:57Z")

</div>

I’m trying to detect multiple events in a system of ODEs. Basically, as soon as any variable crosses a threshold, the value of that variable is changed. I know I can achieve this by explicitly coding up these conditions …

---

## [CairoMakie: how to set appropriate animation margins?](https://discourse.julialang.org/t/cairomakie-how-to-set-appropriate-animation-margins/130049)

<div class="topic-metadata">

**Author:** [@fusion809](https://discourse.julialang.org/u/fusion809)\
**Replies:** 4\
**Last updated:** [June 20, 2025, 5:08am UTC](https://discourse.julialang.org/t/cairomakie-how-to-set-appropriate-animation-margins/130049 "2025-06-20T05:08:28Z")

</div>

I created the below animation of the elastic pendulum using CairoMakie and this code. My only qualm with it is the plot margins; they are too wide. I would like the Cartesian plane shown to span most of the anima…

---

## [Repeated high (6–7) dimensional interpolation](https://discourse.julialang.org/t/repeated-high-6-7-dimensional-interpolation/130042)

<div class="topic-metadata">

**Author:** [@JADekker](https://discourse.julialang.org/u/JADekker)\
**Replies:** 3\
**Last updated:** [June 19, 2025, 9:21pm UTC](https://discourse.julialang.org/t/repeated-high-6-7-dimensional-interpolation/130042 "2025-06-19T21:21:51Z")

</div>

Thank you, I’ll check it out! I can sketch the problem: We are interested in approximating a (costly) function f(x\_1, ..., x\_p) by simulating a grid of points (x\_1^{(n)},...,x\_p^{(n)}),~n = 1,...,N, evaluating the funct…

---

## [SlurmClusterManagers timeout](https://discourse.julialang.org/t/slurmclustermanagers-timeout/130038)

<div class="topic-metadata">

**Author:** [@marcobonici](https://discourse.julialang.org/u/marcobonici)\
**Replies:** 2\
**Last updated:** [June 19, 2025, 6:59pm UTC](https://discourse.julialang.org/t/slurmclustermanagers-timeout/130038 "2025-06-19T18:59:30Z")

</div>

Hello, I have been trying to run some jobs on a Slurm cluster (specifically, the canadian Narval Cluster of the Research Digital Alliance). Possibly due to my inexperience with it, I encountered some issues (I am not o…

---

## [GeoMakie GeoAxis Subticks?](https://discourse.julialang.org/t/geomakie-geoaxis-subticks/130024)

<div class="topic-metadata">

**Author:** [@consumere](https://discourse.julialang.org/u/consumere)\
**Replies:** 2\
**Last updated:** [June 19, 2025, 1:16pm UTC](https://discourse.julialang.org/t/geomakie-geoaxis-subticks/130024 "2025-06-19T13:16:34Z")

</div>

If i run this code, i get wrong axis labels. using GeoMakie set\_theme!(theme\_latexfonts()) fig = Figure(fontsize=12) lonlims = (10,11) latlims = (49,51) ax = GeoAxis(fig\[1, 1\], limits=(lonlims, latlims)) fig Q…

---

## [Local geobox in Albers Equal Area projection in GeoMakie](https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975)

<div class="topic-metadata">

**Author:** [@yum-yab](https://discourse.julialang.org/u/yum-yab)\
**Replies:** 1\
**Last updated:** [June 19, 2025, 12:08pm UTC](https://discourse.julialang.org/t/local-geobox-in-albers-equal-area-projection-in-geomakie/114975 "2025-06-19T12:08:27Z")

</div>

Hi everyone, i want to do something very simple: I want to cut out a local map of europe and the atlantic ocaen, but have troubles doing so in GeoMakie and the Documentation is lacking any good examples that do that. P…

---

## [Solver attributes and set\_optimizer with ParametricOptInterface.jl and JuMP.jl](https://discourse.julialang.org/t/solver-attributes-and-set-optimizer-with-parametricoptinterface-jl-and-jump-jl/129935)

<div class="topic-metadata">

**Author:** [@dcelisgarza](https://discourse.julialang.org/u/dcelisgarza)\
**Replies:** 9\
**Last updated:** [June 19, 2025, 11:44am UTC](https://discourse.julialang.org/t/solver-attributes-and-set-optimizer-with-parametricoptinterface-jl-and-jump-jl/129935 "2025-06-19T11:44:28Z")

</div>

Hi, I could really use ParametricOptInterface.jl (POI) for a library I’m developing (I’ll make an announcement at some point), but I’ve run into a few issues. I don’t know how to provide solver-specific attributes, an…

---

## [Discrepancy of ODE Sensitivity Analysis paper results with benchmarks](https://discourse.julialang.org/t/discrepancy-of-ode-sensitivity-analysis-paper-results-with-benchmarks/129853)

<div class="topic-metadata">

**Author:** [@mra](https://discourse.julialang.org/u/mra)\
**Replies:** 31\
**Last updated:** [June 18, 2025, 2:01pm UTC](https://discourse.julialang.org/t/discrepancy-of-ode-sensitivity-analysis-paper-results-with-benchmarks/129853 "2025-06-18T14:01:19Z")

</div>

I have tried to reproduce results from the paper “A Comparison of Automatic Differentiation and Continuous Sensitivity Analysis for Derivatives of Differential Equation Solutions” (https://arxiv.org/pdf/1812.01892, Fig. …

---

## [Clarification on save/load\_checkpoint() and populate\_from\_csv!()](https://discourse.julialang.org/t/clarification-on-save-load-checkpoint-and-populate-from-csv/129990)

<div class="topic-metadata">

**Author:** [@elisadonati23](https://discourse.julialang.org/u/elisadonati23)\
**Replies:** 0\
**Last updated:** [June 18, 2025, 1:32pm UTC](https://discourse.julialang.org/t/clarification-on-save-load-checkpoint-and-populate-from-csv/129990 "2025-06-18T13:32:23Z")

</div>

Dear all, I’m working on a relatively large agent-based model (ABM) for fish populations, without explicit spatial structure, using Agents.jl. The model includes individual stochasticity in parameter values and uses mul…

---

## [How to get upper triangular Cholesky factor of sparse matrix?](https://discourse.julialang.org/t/how-to-get-upper-triangular-cholesky-factor-of-sparse-matrix/129979)

<div class="topic-metadata">

**Author:** [@slwu89](https://discourse.julialang.org/u/slwu89)\
**Replies:** 3\
**Last updated:** [June 18, 2025, 1:22pm UTC](https://discourse.julialang.org/t/how-to-get-upper-triangular-cholesky-factor-of-sparse-matrix/129979 "2025-06-18T13:22:34Z")

</div>

Hi all. I read the documentation here Linear Algebra · The Julia Language for what is returned when you call cholesky on a SparseMatrixCSC object but I still can’t figure out how to get back the upper triangular factor o…

---

## [A custom type-stable regularization function in Lux & Enzyme](https://discourse.julialang.org/t/a-custom-type-stable-regularization-function-in-lux-enzyme/129984)

<div class="topic-metadata">

**Author:** [@user22](https://discourse.julialang.org/u/user22)\
**Replies:** 0\
**Last updated:** [June 18, 2025, 9:44am UTC](https://discourse.julialang.org/t/a-custom-type-stable-regularization-function-in-lux-enzyme/129984 "2025-06-18T09:44:39Z")

</div>

We are trying to develop a custom regularization function for the biases and weights of a simple dense neural network in Lux. However, we have not found a way to establish type-stable gradient computation for Lux and Enz…

---

## [Nyquist-Shannon sampling theorem question](https://discourse.julialang.org/t/nyquist-shannon-sampling-theorem-question/129892)

<div class="topic-metadata">

**Author:** [@Ju\_ska](https://discourse.julialang.org/u/Ju_ska)\
**Replies:** 7\
**Last updated:** [June 18, 2025, 9:07am UTC](https://discourse.julialang.org/t/nyquist-shannon-sampling-theorem-question/129892 "2025-06-18T09:07:37Z")

</div>

Hi, I always thought that sampling at Fs \>= 2Fmax is the thing to do when sampling a signal. But I came across a school document that says if Fmin is also given, the sampling rate should be Fs \>= 2(Fmax - Fmin). It pu…

---

## [Loading resampled data for PINN in Reactant + Enzyme + Lux](https://discourse.julialang.org/t/loading-resampled-data-for-pinn-in-reactant-enzyme-lux/129974)

<div class="topic-metadata">

**Author:** [@NoFishLikeIan](https://discourse.julialang.org/u/NoFishLikeIan)\
**Replies:** 2\
**Last updated:** [June 17, 2025, 10:58pm UTC](https://discourse.julialang.org/t/loading-resampled-data-for-pinn-in-reactant-enzyme-lux/129974 "2025-06-17T22:58:15Z")

</div>

I am implementing a neural network using Lux and Reactant + Enzyme. At each epoch, I need to resample the data, divide it in batches, and take a training step. For now, I am working on a CPU. What is the best way to do t…

---

## [Rank-deficient jacobian error in Newton-Raphson method for unstable limit cycle refinement](https://discourse.julialang.org/t/rank-deficient-jacobian-error-in-newton-raphson-method-for-unstable-limit-cycle-refinement/129975)

<div class="topic-metadata">

**Author:** [@R.Kanakov](https://discourse.julialang.org/u/R.Kanakov)\
**Replies:** 3\
**Last updated:** [June 17, 2025, 10:20pm UTC](https://discourse.julialang.org/t/rank-deficient-jacobian-error-in-newton-raphson-method-for-unstable-limit-cycle-refinement/129975 "2025-06-17T22:20:57Z")

</div>

I am trying to refine the position of an unstable limit cycle within a chaotic attractor using the Newton-Raphson method, but I am running into a confusing error related to a rank-deficient matrix. To avoid degeneracy o…

---

## [GLMakie displays false colours in WSL2](https://discourse.julialang.org/t/glmakie-displays-false-colours-in-wsl2/129590)

<div class="topic-metadata">

**Author:** [@aris](https://discourse.julialang.org/u/aris)\
**Replies:** 12\
**Last updated:** [June 17, 2025, 2:19pm UTC](https://discourse.julialang.org/t/glmakie-displays-false-colours-in-wsl2/129590 "2025-06-17T14:19:04Z")

</div>

Hello everyone, In WSL2 ubuntu, GLMakie gives me weird colours for some reason. For example, using GLMakie lines(rand(10)) results in: If I save the plot, the colours in the .png image are fine, its just the dis…

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

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