# Specific Domains

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

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

**Page:** 36

---

## [Is there a standard way of building packages with graphical extension](https://discourse.julialang.org/t/is-there-a-standard-way-of-building-packages-with-graphical-extension/130254)

<div class="topic-metadata">

**Author:** [@Herve](https://discourse.julialang.org/u/Herve)\
**Replies:** 3\
**Last updated:** [July 1, 2025, 1:35am UTC](https://discourse.julialang.org/t/is-there-a-standard-way-of-building-packages-with-graphical-extension/130254 "2025-07-01T01:35:13Z")

</div>

Hi. I often bump into the following problem, related to graphics libraries. I have worked with Matlab for more than 20 years, so I had my habits (I nearly stopped since I adopted Julia, and have been all the more happy)…

---

## [Gurobi's Root relaxation objective in MIP Logging is Not a natural LP relaxation one](https://discourse.julialang.org/t/gurobis-root-relaxation-objective-in-mip-logging-is-not-a-natural-lp-relaxation-one/130337)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 3\
**Last updated:** [June 30, 2025, 10:57pm UTC](https://discourse.julialang.org/t/gurobis-root-relaxation-objective-in-mip-logging-is-not-a-natural-lp-relaxation-one/130337 "2025-06-30T22:57:11Z")

</div>

I’m not sure whether this is reasonable, but it is somewhat misleading. Compare the 2 lines where I draw 3 orange balls. julia\> import JuMP, Gurobi; GRB\_ENV = Gurobi.Env(); Set parameter Username Set parameter LicenseI…

---

## [Relation between KernelAbstractions and Adapt](https://discourse.julialang.org/t/relation-between-kernelabstractions-and-adapt/130326)

<div class="topic-metadata">

**Author:** [@bmxam](https://discourse.julialang.org/u/bmxam)\
**Replies:** 1\
**Last updated:** [June 30, 2025, 7:28pm UTC](https://discourse.julialang.org/t/relation-between-kernelabstractions-and-adapt/130326 "2025-06-30T19:28:17Z")

</div>

Hi everyone, I’m trying to move some operations of an existing package on GPU using KernelAbstractions. Since my package defines custom structs, I had to Adapt them, specializing adapt\_structure. Everything works fine …

---

## [Cannot manage to use CUDA.atomic\_add!](https://discourse.julialang.org/t/cannot-manage-to-use-cuda-atomic-add/130328)

<div class="topic-metadata">

**Author:** [@Ludovic\_Dumoulin](https://discourse.julialang.org/u/Ludovic_Dumoulin)\
**Replies:** 4\
**Last updated:** [June 30, 2025, 7:24pm UTC](https://discourse.julialang.org/t/cannot-manage-to-use-cuda-atomic-add/130328 "2025-06-30T19:24:40Z")

</div>

Hello, I want to place particles into a 10x10 grid. To do this, I would like to use CUDA.atomic\_add!() to count how many particles land in each cell. Here is a minimal example: using CUDA function kernel\_test\_atomicad…

---

## [How to use Lux's \`apply\` function?](https://discourse.julialang.org/t/how-to-use-luxs-apply-function/130339)

<div class="topic-metadata">

**Author:** [@Torkel](https://discourse.julialang.org/u/Torkel)\
**Replies:** 2\
**Last updated:** [June 30, 2025, 5:01pm UTC](https://discourse.julialang.org/t/how-to-use-luxs-apply-function/130339 "2025-06-30T17:01:08Z")

</div>

I am trying to get Lux to work, but am stumbling at a relatively early stage of the tutorials. I am trying to use a combination of the quickstart and getting started pages (Getting Started | Lux.jl Docs, Getting Started …

---

## [Deploying an interactive WGLMakie.jl example on the internet](https://discourse.julialang.org/t/deploying-an-interactive-wglmakie-jl-example-on-the-internet/130313)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 4\
**Last updated:** [June 30, 2025, 4:16pm UTC](https://discourse.julialang.org/t/deploying-an-interactive-wglmakie-jl-example-on-the-internet/130313 "2025-06-30T16:16:21Z")

</div>

I have a simple code (see OptimalFlightPaths.jl) that displays a plot with some sliders, and you can change the plot using the sliders. How can I display this on the internet in a safe way?

---

## [Solve linear system under Hermitian constraint](https://discourse.julialang.org/t/solve-linear-system-under-hermitian-constraint/130158)

<div class="topic-metadata">

**Author:** [@araujoms](https://discourse.julialang.org/u/araujoms)\
**Replies:** 25\
**Last updated:** [June 30, 2025, 1:26pm UTC](https://discourse.julialang.org/t/solve-linear-system-under-hermitian-constraint/130158 "2025-06-30T13:26:05Z")

</div>

Is there a nice way to solve the linear system XA = B under the constraint X = X^\\dagger? Here A and B are complex matrices that are not necessarily square. What I came up with is to solve YA = B and A^\\dagger Y = B^\\da…

---

## [PackageCompiler fails to create app for MadNLPGPU + ExaModels (CUDSS linear solver)](https://discourse.julialang.org/t/packagecompiler-fails-to-create-app-for-madnlpgpu-examodels-cudss-linear-solver/125998)

<div class="topic-metadata">

**Author:** [@KSepetanc](https://discourse.julialang.org/u/KSepetanc)\
**Replies:** 11\
**Last updated:** [June 30, 2025, 7:26am UTC](https://discourse.julialang.org/t/packagecompiler-fails-to-create-app-for-madnlpgpu-examodels-cudss-linear-solver/125998 "2025-06-30T07:26:00Z")

</div>

I can not create “app” with PackageCompiler even for simplest GPU optimization with ExaModels+MadNLPGPU running cuDSS. I’ve made full MWE project with package to be built in attempt to help this to debug. Compiling does …

---

## [What is the intended use for \`Turing.returned(model)\`?](https://discourse.julialang.org/t/what-is-the-intended-use-for-turing-returned-model/130105)

<div class="topic-metadata">

**Author:** [@lunafazio](https://discourse.julialang.org/u/lunafazio)\
**Replies:** 4\
**Last updated:** [June 30, 2025, 7:16am UTC](https://discourse.julialang.org/t/what-is-the-intended-use-for-turing-returned-model/130105 "2025-06-30T07:16:24Z")

</div>

There are no examples in the documentation, but based on its description, Return a model wrapper indicating that it is a model over its return-values. I assumed that this would allow me to do something like chain = …

---

## [Gurobi's trial solution is numerically weird on a primitive CTPLN model](https://discourse.julialang.org/t/gurobis-trial-solution-is-numerically-weird-on-a-primitive-ctpln-model/130294)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 1\
**Last updated:** [June 29, 2025, 11:11pm UTC](https://discourse.julialang.org/t/gurobis-trial-solution-is-numerically-weird-on-a-primitive-ctpln-model/130294 "2025-06-29T23:11:27Z")

</div>

Although this might not necessarily be deemed a bug, it’s certainly a bit weird. Why? julia\> import JuMP, Gurobi; GRB\_ENV = Gurobi.Env(); Set parameter Username Set parameter LicenseID to value 2663413 Academic license …

---

## [Optim.jl : user defined Manifold, error in the definition of the retraction](https://discourse.julialang.org/t/optim-jl-user-defined-manifold-error-in-the-definition-of-the-retraction/130297)

<div class="topic-metadata">

**Author:** [@aitzkora](https://discourse.julialang.org/u/aitzkora)\
**Replies:** 2\
**Last updated:** [June 29, 2025, 8:49pm UTC](https://discourse.julialang.org/t/optim-jl-user-defined-manifold-error-in-the-definition-of-the-retraction/130297 "2025-06-29T20:49:42Z")

</div>

Hi, I try to optimize a simple quadratic function on a user-defined manifold (the hyperplan orthogonal to ones(n)) with the following code using Optim const c = ones(2) ./ sqrt(2) f(x) = 0.5 \* sum((x - c).^2) g!(g,x) =…

---

## [Get Turing's model predictions after optimization](https://discourse.julialang.org/t/get-turings-model-predictions-after-optimization/130275)

<div class="topic-metadata">

**Author:** [@khoffie](https://discourse.julialang.org/u/khoffie)\
**Replies:** 6\
**Last updated:** [June 29, 2025, 7:01pm UTC](https://discourse.julialang.org/t/get-turings-model-predictions-after-optimization/130275 "2025-06-29T19:01:03Z")

</div>

Hi, since Turing 0.35 generated\_quantities() is deprecated. So how do I now get my model’s predictions? Before 0.35 I did something like this function predict(mdl, mles) ## mdl is an instantiated Turing model pa…

---

## [Morris Screening Method reduces the number of parameters](https://discourse.julialang.org/t/morris-screening-method-reduces-the-number-of-parameters/130258)

<div class="topic-metadata">

**Author:** [@AlirezaMiraliakbar](https://discourse.julialang.org/u/AlirezaMiraliakbar)\
**Replies:** 1\
**Last updated:** [June 27, 2025, 11:09am UTC](https://discourse.julialang.org/t/morris-screening-method-reduces-the-number-of-parameters/130258 "2025-06-27T11:09:22Z")

</div>

Hello Everyone, I am new to GlobalSensitivity.jl. I am trying to find the most sensitive parameters in my large ODE system which has 37 parameters. I gave all parameter ranges and set p\_steps 10 for all parameters in a …

---

## [How to get value and matrix index at a given (x;y) point?](https://discourse.julialang.org/t/how-to-get-value-and-matrix-index-at-a-given-x-y-point/129692)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 9\
**Last updated:** [June 27, 2025, 5:04am UTC](https://discourse.julialang.org/t/how-to-get-value-and-matrix-index-at-a-given-x-y-point/129692 "2025-06-27T05:04:10Z")

</div>

I have a set of rasters with the same grid and a large number of points and need to create a CZV with these points + information from the rasters (so my questions is similar to this one, but I am using the Rasters.jl pac…

---

## [Sersic profile](https://discourse.julialang.org/t/sersic-profile/129894)

<div class="topic-metadata">

**Author:** [@Cyan](https://discourse.julialang.org/u/Cyan)\
**Replies:** 1\
**Last updated:** [June 27, 2025, 1:57am UTC](https://discourse.julialang.org/t/sersic-profile/129894 "2025-06-27T01:57:46Z")

</div>

Does anyone know whether there is package with Sersic profile and its derivative/Hession profile?

---

## [JuMP diet tutorial issue?](https://discourse.julialang.org/t/jump-diet-tutorial-issue/130231)

<div class="topic-metadata">

**Author:** [@Vee\_Kalkunte](https://discourse.julialang.org/u/Vee_Kalkunte)\
**Replies:** 2\
**Last updated:** [June 26, 2025, 9:22pm UTC](https://discourse.julialang.org/t/jump-diet-tutorial-issue/130231 "2025-06-26T21:22:06Z")

</div>

Hi! I’m following the JuMP documentation’s Diet problem code, and I’m running into a constraint error (“ERROR: At code.jl:66: @constraint(model, row.min \<= sum(foods\[!, row.nutrient\] .\* foods.x) \<= row.max): Unrecognize…

---

## [Custom differentiation with Flux](https://discourse.julialang.org/t/custom-differentiation-with-flux/130249)

<div class="topic-metadata">

**Author:** [@Fabrice\_Rosay](https://discourse.julialang.org/u/Fabrice_Rosay)\
**Replies:** 2\
**Last updated:** [June 26, 2025, 7:50pm UTC](https://discourse.julialang.org/t/custom-differentiation-with-flux/130249 "2025-06-26T19:50:32Z")

</div>

Currently the derivative of roundis 0. I would like it to be 1 instead and also that it works when broadcasted. So i wrote a custom my\_roundand tried to write the new rule with adjoint and @scalar\_rule but it does not fu…

---

## [Slow Turing.jl sampling compared to python pymc](https://discourse.julialang.org/t/slow-turing-jl-sampling-compared-to-python-pymc/130228)

<div class="topic-metadata">

**Author:** [@mahe](https://discourse.julialang.org/u/mahe)\
**Replies:** 7\
**Last updated:** [June 26, 2025, 2:59pm UTC](https://discourse.julialang.org/t/slow-turing-jl-sampling-compared-to-python-pymc/130228 "2025-06-26T14:59:02Z")

</div>

I have a spatial model on a (currently) 5x5 grid, for two variables (total size 3737). And 483 observations for one variables. In this simple setting, I want to work in full space (no space reduction), and model my syste…

---

## [Copying a JuMP model is slower than building it from scratch](https://discourse.julialang.org/t/copying-a-jump-model-is-slower-than-building-it-from-scratch/130190)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 6\
**Last updated:** [June 26, 2025, 2:40pm UTC](https://discourse.julialang.org/t/copying-a-jump-model-is-slower-than-building-it-from-scratch/130190 "2025-06-26T14:40:39Z")

</div>

As the title. Then what is the point of the Base.copy(model::AbstractModel) function? julia\> Min\_ȷ\_fixed\_part Min\_ȷ\_fixed\_part (generic function with 1 method) julia\> t = time(); model = Min\_ȷ\_fixed\_part(); time() - t …

---

## [Solving a large sparse ODE using Exponential methods](https://discourse.julialang.org/t/solving-a-large-sparse-ode-using-exponential-methods/125159)

<div class="topic-metadata">

**Author:** [@albertomercurio](https://discourse.julialang.org/u/albertomercurio)\
**Replies:** 8\
**Last updated:** [June 25, 2025, 7:49pm UTC](https://discourse.julialang.org/t/solving-a-large-sparse-ode-using-exponential-methods/125159 "2025-06-25T19:49:29Z")

</div>

Hello, I need to integrate a differential equation of the form \\dot{u} = L(t) u \\, , where u is a vector, L(t) = L\_0 + c(t) \* L\_1 \\, , L\_0 and L\_1 large sparse matrices, and c(t) a time dependent coefficient. I k…

---

## [Heterogeneous random seeding](https://discourse.julialang.org/t/heterogeneous-random-seeding/130196)

<div class="topic-metadata">

**Author:** [@0samuraiE](https://discourse.julialang.org/u/0samuraiE)\
**Replies:** 1\
**Last updated:** [June 25, 2025, 3:00pm UTC](https://discourse.julialang.org/t/heterogeneous-random-seeding/130196 "2025-06-25T15:00:20Z")

</div>

Hello all, Is there a unified way to set the random seed across heterogeneous devices, such as CPU and GPU? I’d be happy if GPUArrays.jl or KernelAbstractions.jl provides such a method. julia\> using CUDA julia\> using…

---

## [Bifurcation diagram 2D reaction diffusion type system on GPU](https://discourse.julialang.org/t/bifurcation-diagram-2d-reaction-diffusion-type-system-on-gpu/128490)

<div class="topic-metadata">

**Author:** [@mellever](https://discourse.julialang.org/u/mellever)\
**Replies:** 5\
**Last updated:** [June 25, 2025, 9:48am UTC](https://discourse.julialang.org/t/bifurcation-diagram-2d-reaction-diffusion-type-system-on-gpu/128490 "2025-06-25T09:48:56Z")

</div>

Hi all, This is my first time posting, so forgive me if I am not following the proper etiquette. I am currently working on computing bifurcation diagrams of 2D reaction diffusion type models (e.g. Brusselator model) wi…

---

## [IPOPT suppress output file generation](https://discourse.julialang.org/t/ipopt-suppress-output-file-generation/127595)

<div class="topic-metadata">

**Author:** [@hjalgra](https://discourse.julialang.org/u/hjalgra)\
**Replies:** 22\
**Last updated:** [June 25, 2025, 7:03am UTC](https://discourse.julialang.org/t/ipopt-suppress-output-file-generation/127595 "2025-06-25T07:03:16Z")

</div>

I am using IPOPT to solve a sequence of optimization problems (interfaced using SNOW). In order to prevent the creation of an output file I am initializing IPOPT using the following options ip\_options = Dict( …

---

## [Agent based network model using EventQueueABM](https://discourse.julialang.org/t/agent-based-network-model-using-eventqueueabm/129889)

<div class="topic-metadata">

**Author:** [@zoek](https://discourse.julialang.org/u/zoek)\
**Replies:** 15\
**Last updated:** [June 24, 2025, 7:33pm UTC](https://discourse.julialang.org/t/agent-based-network-model-using-eventqueueabm/129889 "2025-06-24T19:33:02Z")

</div>

Hello, I’m trying to use EventQueueABM from agents.jl on a completely connected network to model disease spread (susceptible, exposed, infected and recovered states). Since it’s completely connected, I believe I should b…

---

## [Using JuMP for solving a symbolic algebraic system](https://discourse.julialang.org/t/using-jump-for-solving-a-symbolic-algebraic-system/130153)

<div class="topic-metadata">

**Author:** [@jonlym](https://discourse.julialang.org/u/jonlym)\
**Replies:** 2\
**Last updated:** [June 24, 2025, 7:00pm UTC](https://discourse.julialang.org/t/using-jump-for-solving-a-symbolic-algebraic-system/130153 "2025-06-24T19:00:05Z")

</div>

Hi All, I’d like to solve a system of algebraic equations with the flexibility to specify one variable as independent and one as the dependent variable. In the toy example below, I am calculating the Reynolds number, Re …

---

## [KNITRO license path](https://discourse.julialang.org/t/knitro-license-path/130145)

<div class="topic-metadata">

**Author:** [@CeterisPartybus](https://discourse.julialang.org/u/CeterisPartybus)\
**Replies:** 4\
**Last updated:** [June 24, 2025, 6:42pm UTC](https://discourse.julialang.org/t/knitro-license-path/130145 "2025-06-24T18:42:53Z")

</div>

From where does Julia in VS Code source environmental variables setting the KNITRO license path ENV\[“ARTELYS\_LICENSE”\]? The issue: I updated my KNITRO license file and its location and set the environment variable in .…

---

## [Solvers for Block-Tridiagonal](https://discourse.julialang.org/t/solvers-for-block-tridiagonal/130013)

<div class="topic-metadata">

**Author:** [@rsenne](https://discourse.julialang.org/u/rsenne)\
**Replies:** 14\
**Last updated:** [June 24, 2025, 6:41pm UTC](https://discourse.julialang.org/t/solvers-for-block-tridiagonal/130013 "2025-06-24T18:41:55Z")

</div>

Block Tridiagonal Linear System Performance: Custom vs Package Implementations I recently submitted an issue to BlockBandedMatrices.jl about this, but thought it might be interesting to the broader community. Problem I …

---

## [Matrix multiplication and vector multiplication give different results!](https://discourse.julialang.org/t/matrix-multiplication-and-vector-multiplication-give-different-results/130177)

<div class="topic-metadata">

**Author:** [@tomtom](https://discourse.julialang.org/u/tomtom)\
**Replies:** 4\
**Last updated:** [June 24, 2025, 5:04pm UTC](https://discourse.julialang.org/t/matrix-multiplication-and-vector-multiplication-give-different-results/130177 "2025-06-24T17:04:11Z")

</div>

I found that (slightly) different results coming from matrix-matrix vs. matrix-vector multiplication! julia\> m1 = randn(200, 300); julia\> m2 = randn(300, 400); julia\> m12 = m1 \* m2; # matrix-matrix multiplicati…

---

## [Training Neural ODEs - Advice](https://discourse.julialang.org/t/training-neural-odes-advice/130141)

<div class="topic-metadata">

**Author:** [@julianewbsourcream](https://discourse.julialang.org/u/julianewbsourcream)\
**Replies:** 2\
**Last updated:** [June 24, 2025, 1:49pm UTC](https://discourse.julialang.org/t/training-neural-odes-advice/130141 "2025-06-24T13:49:26Z")

</div>

Hi, I am looking for some advice about training neural ODE models. I am looking to train a neural ODE model on a sample of ~25000, 10x4000 matrices. The columns represent time so I was looking to at use an increasing am…

---

## [Sorting when the comparison operator depends on the entire input vector](https://discourse.julialang.org/t/sorting-when-the-comparison-operator-depends-on-the-entire-input-vector/130152)

<div class="topic-metadata">

**Author:** [@colintbowers](https://discourse.julialang.org/u/colintbowers)\
**Replies:** 6\
**Last updated:** [June 24, 2025, 6:04am UTC](https://discourse.julialang.org/t/sorting-when-the-comparison-operator-depends-on-the-entire-input-vector/130152 "2025-06-24T06:04:22Z")

</div>

Hi all. Given a parameter Vector{Float64} called x, I want to find x\* such that f(x\*) \> f(x) for all possible permutations of x, where f maps x to a single Float64. I could (and have successfully) used a general optimiza…

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

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