# Specific Domains

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

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

**Page:** 209

---

## [Makie and X11 forwarding](https://discourse.julialang.org/t/makie-and-x11-forwarding/91975)

<div class="topic-metadata">

**Author:** [@wsshin](https://discourse.julialang.org/u/wsshin)\
**Replies:** 7\
**Last updated:** [February 27, 2023, 9:46pm UTC](https://discourse.julialang.org/t/makie-and-x11-forwarding/91975 "2023-02-27T21:46:33Z")

</div>

I am developing in a remote server. I am connected to the server from my local machine via SSH. To enable X11 forwarding, I connect using the -Y option (i.e., ssh -Y). This allows me to use the GUI interface of the app…

---

## [Error in using KNITRO with JuMP](https://discourse.julialang.org/t/error-in-using-knitro-with-jump/95951)

<div class="topic-metadata">

**Author:** [@Rick](https://discourse.julialang.org/u/Rick)\
**Replies:** 6\
**Last updated:** [March 12, 2023, 9:51pm UTC](https://discourse.julialang.org/t/error-in-using-knitro-with-jump/95951 "2023-03-12T21:51:31Z")

</div>

Hi everyone, I have just updated my packages, then I am unable to use KNITRO.jl with JuMP anymore. The code I run is: using KNITRO, JuMP model = Model(KNITRO.Optimizer) Then, I get the following error: The provided …

---

## [Flux with Matrix input dimensions](https://discourse.julialang.org/t/flux-with-matrix-input-dimensions/43610)

<div class="topic-metadata">

**Author:** [@Ruan\_Stander](https://discourse.julialang.org/u/Ruan_Stander)\
**Replies:** 2\
**Last updated:** [March 12, 2023, 6:30pm UTC](https://discourse.julialang.org/t/flux-with-matrix-input-dimensions/43610 "2023-03-12T18:30:15Z")

</div>

Hi Just started using flux, cant get on top of input/layer dimensions where input is a matrix Adapting https://github.com/FluxML/model-zoo/blob/master/vision/mnist/mlp.jl to my input space but its not working My Input…

---

## [What happened to XLA.jl](https://discourse.julialang.org/t/what-happened-to-xla-jl/88088)

<div class="topic-metadata">

**Author:** [@Honza9723](https://discourse.julialang.org/u/Honza9723)\
**Replies:** 16\
**Last updated:** [March 12, 2023, 3:24pm UTC](https://discourse.julialang.org/t/what-happened-to-xla-jl/88088 "2023-03-12T15:24:22Z")

</div>

Dear All, I would like to ask, what happened to XLA.jl. When I checked the repository, it says it was archived. Is the project of building the Julia-to-XLA interface dead? Best, Jan

---

## [How to implement a model with missing observations in Turing?](https://discourse.julialang.org/t/how-to-implement-a-model-with-missing-observations-in-turing/95572)

<div class="topic-metadata">

**Author:** [@AdamR](https://discourse.julialang.org/u/AdamR)\
**Replies:** 6\
**Last updated:** [March 12, 2023, 12:04pm UTC](https://discourse.julialang.org/t/how-to-implement-a-model-with-missing-observations-in-turing/95572 "2023-03-12T12:04:13Z")

</div>

I would like to write in Julia the following toy Bernoulli model written in JAGS: model{ p0 ~ dbeta(1,1) for(i in 1:n){ y\[i\] ~ dbern(p0) } } Here’s the complete driver true\_p0 \<- 0.2 n \<- 600 y \<- sampl…

---

## [Best practice for element-wise vector value update without memory allocation](https://discourse.julialang.org/t/best-practice-for-element-wise-vector-value-update-without-memory-allocation/94864)

<div class="topic-metadata">

**Author:** [@Charles](https://discourse.julialang.org/u/Charles)\
**Replies:** 3\
**Last updated:** [March 12, 2023, 10:14am UTC](https://discourse.julialang.org/t/best-practice-for-element-wise-vector-value-update-without-memory-allocation/94864 "2023-03-12T10:14:12Z")

</div>

I am implementing Gaussian Elimination using Julia. The important thing is, when doing the elimination the row vectors are kept updating row by row. While writing a for loop to iterate over the row vector and assigning n…

---

## [ArchGDAL.read error: GDALError (CE\_Failure, code 4):not recognized as a supported file format](https://discourse.julialang.org/t/archgdal-read-error-gdalerror-ce-failure-code-4-not-recognized-as-a-supported-file-format/95959)

<div class="topic-metadata">

**Author:** [@sunayana](https://discourse.julialang.org/u/sunayana)\
**Replies:** 4\
**Last updated:** [March 12, 2023, 10:08am UTC](https://discourse.julialang.org/t/archgdal-read-error-gdalerror-ce-failure-code-4-not-recognized-as-a-supported-file-format/95959 "2023-03-12T10:08:44Z")

</div>

Hi, I am currently stuck with a strange issue in GDAL and do not understand how to solve it. I have a method inside an IO module: read\_image\_file(filename::AbstractPath) = AG.read(string(filename)) This method is bein…

---

## [Problem Using annotate! to add info to an animated plot but it overlaps and cannot be read](https://discourse.julialang.org/t/problem-using-annotate-to-add-info-to-an-animated-plot-but-it-overlaps-and-cannot-be-read/95935)

<div class="topic-metadata">

**Author:** [@J\_Braulio\_Ramz\_G](https://discourse.julialang.org/u/J_Braulio_Ramz_G)\
**Replies:** 1\
**Last updated:** [March 11, 2023, 9:05pm UTC](https://discourse.julialang.org/t/problem-using-annotate-to-add-info-to-an-animated-plot-but-it-overlaps-and-cannot-be-read/95935 "2023-03-11T21:05:39Z")

</div>

Hi I am triying to use !annotate() to add text info like current iteration, aproximation and relative error to an animation of Newton-Raphson method but the previous info is accumulating in the new frames and overlaps.…

---

## [PCA on one time series variable into components](https://discourse.julialang.org/t/pca-on-one-time-series-variable-into-components/95915)

<div class="topic-metadata">

**Author:** [@BMval](https://discourse.julialang.org/u/BMval)\
**Replies:** 3\
**Last updated:** [March 11, 2023, 8:54pm UTC](https://discourse.julialang.org/t/pca-on-one-time-series-variable-into-components/95915 "2023-03-11T20:54:54Z")

</div>

I tried to use PCA from MultivariateStats, and I got a little bit stuck. I have a one-time series variable - Y. I need to decompose it into orthogonal factors (by using PCA). It looks like a number of factors can be l…

---

## [Maximizing a term with barrier in JuMP](https://discourse.julialang.org/t/maximizing-a-term-with-barrier-in-jump/95884)

<div class="topic-metadata">

**Author:** [@Shayan](https://discourse.julialang.org/u/Shayan)\
**Replies:** 7\
**Last updated:** [March 11, 2023, 5:39pm UTC](https://discourse.julialang.org/t/maximizing-a-term-with-barrier-in-jump/95884 "2023-03-11T17:39:01Z")

</div>

I’m sorry that the title isn’t clear enough; I promise to provide a suitable description with a minimal reproducible example. I want to optimize the following model to achieve the optimal set of w\_j (which is a vector): …

---

## [How do I avoid atomics in this kind of code on GPU?](https://discourse.julialang.org/t/how-do-i-avoid-atomics-in-this-kind-of-code-on-gpu/95881)

<div class="topic-metadata">

**Author:** [@Ahmed\_Salih](https://discourse.julialang.org/u/Ahmed_Salih)\
**Replies:** 2\
**Last updated:** [March 11, 2023, 8:39am UTC](https://discourse.julialang.org/t/how-do-i-avoid-atomics-in-this-kind-of-code-on-gpu/95881 "2023-03-11T08:39:14Z")

</div>

I want to stress that this is a toy example. It is not a real calculation. It is a toy example, intended to be a minimal working example. It is showing some approach I want to apply to a problem of mine Hello! I am wri…

---

## [Interpolations of a function that has NaN](https://discourse.julialang.org/t/interpolations-of-a-function-that-has-nan/95704)

<div class="topic-metadata">

**Author:** [@roi.holtzman](https://discourse.julialang.org/u/roi.holtzman)\
**Replies:** 9\
**Last updated:** [March 11, 2023, 1:29am UTC](https://discourse.julialang.org/t/interpolations-of-a-function-that-has-nan/95704 "2023-03-11T01:29:11Z")

</div>

I am trying to interpolate a 2d function that is defined only on a partial set of my parameters. Where it is not defined, my 2d array has NaNs. Note that in my example the range of where the function is defined is not a…

---

## [Lorentz Vector reusability for QED.jl (event generator)](https://discourse.julialang.org/t/lorentz-vector-reusability-for-qed-jl-event-generator/95871)

<div class="topic-metadata">

**Author:** [@jling](https://discourse.julialang.org/u/jling)\
**Replies:** 2\
**Last updated:** [March 10, 2023, 11:54pm UTC](https://discourse.julialang.org/t/lorentz-vector-reusability-for-qed-jl-event-generator/95871 "2023-03-10T23:54:13Z")

</div>

I was watching in awe of https://www.youtube.com/watch?v=BqRps2dwLxE by @szabo137. Afterwards @peremato asked one of the questions I had which is why not re-use some of the existing JuliaHEP stuff, simplest example being…

---

## [(help request) using ArchGDAL.gdalwarp() for resampling a raster dataset](https://discourse.julialang.org/t/help-request-using-archgdal-gdalwarp-for-resampling-a-raster-dataset/47039)

<div class="topic-metadata">

**Author:** [@vlandau](https://discourse.julialang.org/u/vlandau)\
**Replies:** 9\
**Last updated:** [March 10, 2023, 9:00pm UTC](https://discourse.julialang.org/t/help-request-using-archgdal-gdalwarp-for-resampling-a-raster-dataset/47039 "2023-03-10T21:00:08Z")

</div>

Hi all, I’m struggling to figure out how to use ArchGDAL.gdalwarp (currently only unsafe\_gdalwarp is documented). ArchGDAL.gdalwarp’s first parameter being a function is a bit confusing as well. Here’s the GDAL CLI vers…

---

## [Constriant in ADNLPModels package](https://discourse.julialang.org/t/constriant-in-adnlpmodels-package/95832)

<div class="topic-metadata">

**Author:** [@X\_Chen](https://discourse.julialang.org/u/X_Chen)\
**Replies:** 4\
**Last updated:** [March 10, 2023, 8:29pm UTC](https://discourse.julialang.org/t/constriant-in-adnlpmodels-package/95832 "2023-03-10T20:29:09Z")

</div>

Hi everyone, I am trying to do a nonlinear optimization with six parameters using the ADNLModles package. For example: function constriants(para) cons=1-para\[1\]^2-para\[2\]^2-para\[3\]^2+2\*para\[1\]\*para\[2\]\*para\[3\] r…

---

## [JuMP/CPLEX: seemingly not all lazy cuts are respected](https://discourse.julialang.org/t/jump-cplex-seemingly-not-all-lazy-cuts-are-respected/95798)

<div class="topic-metadata">

**Author:** [@mike\_k](https://discourse.julialang.org/u/mike_k)\
**Replies:** 4\
**Last updated:** [March 10, 2023, 5:44pm UTC](https://discourse.julialang.org/t/jump-cplex-seemingly-not-all-lazy-cuts-are-respected/95798 "2023-03-10T17:44:41Z")

</div>

Dear folks, it seems that CPLEX sometimes does not respect all lazy cuts in the final solution. I know that there is no guarantee that it respects submitted cuts immediately, but it should do in the final solution afaik…

---

## [Cannot show that the matrix form of my SPD preconditioner is SPD](https://discourse.julialang.org/t/cannot-show-that-the-matrix-form-of-my-spd-preconditioner-is-spd/95869)

<div class="topic-metadata">

**Author:** [@nvenkov1](https://discourse.julialang.org/u/nvenkov1)\
**Replies:** 4\
**Last updated:** [March 10, 2023, 4:12pm UTC](https://discourse.julialang.org/t/cannot-show-that-the-matrix-form-of-my-spd-preconditioner-is-spd/95869 "2023-03-10T16:12:50Z")

</div>

I define a symmetric positive definite (SPD) matrix A. Then, based on this matrix, I define two preconditioners, which should also be SPD, namely an AMG and a Cholesky factorization. Both preconditioners are matrix-free.…

---

## [ControlSystemIdentification](https://discourse.julialang.org/t/controlsystemidentification/95794)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 12\
**Last updated:** [March 10, 2023, 1:18pm UTC](https://discourse.julialang.org/t/controlsystemidentification/95794 "2023-03-10T13:18:07Z")

</div>

Hello, I am trying out the package ControlSystemIdentification. It looks very good, each step very well explained! One question: In the example Temperature control · ControlSystemIdentification Documentation you are u…

---

## [How to create SecondOrderODEProblem from ODESystem](https://discourse.julialang.org/t/how-to-create-secondorderodeproblem-from-odesystem/95847)

<div class="topic-metadata">

**Author:** [@TheLateKronos](https://discourse.julialang.org/u/TheLateKronos)\
**Replies:** 2\
**Last updated:** [March 10, 2023, 8:55am UTC](https://discourse.julialang.org/t/how-to-create-secondorderodeproblem-from-odesystem/95847 "2023-03-10T08:55:10Z")

</div>

I want to simulate a single body orbiting (0, 0) by VelocityVerlet() and using a system defined through ModelingToolkit. Functional code without modeling toolkit is as follows: using DifferentialEquations function ddu\_k…

---

## [Fit Gaussian process to training points](https://discourse.julialang.org/t/fit-gaussian-process-to-training-points/90841)

<div class="topic-metadata">

**Author:** [@jar1](https://discourse.julialang.org/u/jar1)\
**Replies:** 7\
**Last updated:** [January 14, 2023, 7:01pm UTC](https://discourse.julialang.org/t/fit-gaussian-process-to-training-points/90841 "2023-01-14T19:01:22Z")

</div>

For self-study I am trying to fit a GP according to the formula in the GPML book I am able to draw some random functions from the prior distribution, but when I try to condition on the data it seems to make no differ…

---

## [Light position for 3D surface plots with bands](https://discourse.julialang.org/t/light-position-for-3d-surface-plots-with-bands/95823)

<div class="topic-metadata">

**Author:** [@yhchang96](https://discourse.julialang.org/u/yhchang96)\
**Replies:** 1\
**Last updated:** [March 10, 2023, 4:57am UTC](https://discourse.julialang.org/t/light-position-for-3d-surface-plots-with-bands/95823 "2023-03-10T04:57:10Z")

</div>

I tried following the example in this tutorial to create a 3D surface plot with a ‘side band’: using GLMakie; GLMakie.activate!() # v.0.8.2 cmap = Reverse(:Spectral\_4) lightpos1 = Vec3f(4, -5, 5) set\_theme!() update…

---

## [SciMLSensitivity.jl : adjoint\_sensitivities and CVODE\_BDF crashes for medium sized stiff ODE model](https://discourse.julialang.org/t/scimlsensitivity-jl-adjoint-sensitivities-and-cvode-bdf-crashes-for-medium-sized-stiff-ode-model/95667)

<div class="topic-metadata">

**Author:** [@sebapersson](https://discourse.julialang.org/u/sebapersson)\
**Replies:** 11\
**Last updated:** [March 9, 2023, 11:35pm UTC](https://discourse.julialang.org/t/scimlsensitivity-jl-adjoint-sensitivities-and-cvode-bdf-crashes-for-medium-sized-stiff-ode-model/95667 "2023-03-09T23:35:43Z")

</div>

We have recently done some benchmarks to test the performance of DifferentialEquations.jl for parameter estimation for ODE models in the field of system biology. Overall, DifferentialEquations has performed great, but we…

---

## [Dot broadcasting of an assignment to a tuple](https://discourse.julialang.org/t/dot-broadcasting-of-an-assignment-to-a-tuple/95829)

<div class="topic-metadata">

**Author:** [@dmbates](https://discourse.julialang.org/u/dmbates)\
**Replies:** 1\
**Last updated:** [March 9, 2023, 11:09pm UTC](https://discourse.julialang.org/t/dot-broadcasting-of-an-assignment-to-a-tuple/95829 "2023-03-09T23:09:20Z")

</div>

The context of this question is an algorithm called iteratively re-weighted least squares (IRLS) used to fit generalized linear models (GLMs). At each iteration vectors of the mean responses, μ, the unit deviances, dev, …

---

## [How to make a contour plot from a parametric surface in Makie?](https://discourse.julialang.org/t/how-to-make-a-contour-plot-from-a-parametric-surface-in-makie/95818)

<div class="topic-metadata">

**Author:** [@leogabac](https://discourse.julialang.org/u/leogabac)\
**Replies:** 1\
**Last updated:** [March 9, 2023, 7:14pm UTC](https://discourse.julialang.org/t/how-to-make-a-contour-plot-from-a-parametric-surface-in-makie/95818 "2023-03-09T19:14:40Z")

</div>

Hello! Recently for a project I was wondering if there is any way to visualize the contour plot of a parametric surface in Makie. Typically, if I have the explicit expression f(x,y), I can just evaluate in a domain and …

---

## [Holomorphic function not being recognized as such by ChainRules.jl](https://discourse.julialang.org/t/holomorphic-function-not-being-recognized-as-such-by-chainrules-jl/95806)

<div class="topic-metadata">

**Author:** [@Hareruya](https://discourse.julialang.org/u/Hareruya)\
**Replies:** 2\
**Last updated:** [March 9, 2023, 5:11pm UTC](https://discourse.julialang.org/t/holomorphic-function-not-being-recognized-as-such-by-chainrules-jl/95806 "2023-03-09T17:11:56Z")

</div>

Greetings. I am currently trying to implement some code using @scalar\_rule for a custom gradient for a very important group of functions in Physics and Engineering, the spherical bessel functions, in this case particula…

---

## [Help with nested optimization](https://discourse.julialang.org/t/help-with-nested-optimization/95453)

<div class="topic-metadata">

**Author:** [@nilshg](https://discourse.julialang.org/u/nilshg)\
**Replies:** 18\
**Last updated:** [March 9, 2023, 10:28am UTC](https://discourse.julialang.org/t/help-with-nested-optimization/95453 "2023-03-09T10:28:11Z")

</div>

After a few weeks of thinking about this on and off with a lot of help from the good folks in Slack’s #math-optimization channel I’ve decided I don’t trust myself to come up with an optimal (no pun intended) solution to …

---

## [Rolling/running functions with complex output on multiple variables](https://discourse.julialang.org/t/rolling-running-functions-with-complex-output-on-multiple-variables/95723)

<div class="topic-metadata">

**Author:** [@mreichMPI-BGC](https://discourse.julialang.org/u/mreichMPI-BGC)\
**Replies:** 6\
**Last updated:** [March 9, 2023, 9:55am UTC](https://discourse.julialang.org/t/rolling-running-functions-with-complex-output-on-multiple-variables/95723 "2023-03-09T09:55:01Z")

</div>

Dear Julia community, I want to run a rolling function with multiple output on multiple columns of a DataFrame. RollingFunctions.jl seems to allow only scalar output. ImageFiltering.mapwindow only allows as input a sin…

---

## [How to assign a device for each worker correctly in a multi gpu multi node scenario?](https://discourse.julialang.org/t/how-to-assign-a-device-for-each-worker-correctly-in-a-multi-gpu-multi-node-scenario/95715)

<div class="topic-metadata">

**Author:** [@hunKing](https://discourse.julialang.org/u/hunKing)\
**Replies:** 6\
**Last updated:** [March 9, 2023, 1:56am UTC](https://discourse.julialang.org/t/how-to-assign-a-device-for-each-worker-correctly-in-a-multi-gpu-multi-node-scenario/95715 "2023-03-09T01:56:48Z")

</div>

Hello folks, Does someone know how to correctly set the GPUs in each node correctly to individual workers in those nodes? For instance, imagine I have 2 nodes, with each having 2 GPUs. If I try to use the same procedure…

---

## [Catalyst.jl: Algebraically relate two species](https://discourse.julialang.org/t/catalyst-jl-algebraically-relate-two-species/95762)

<div class="topic-metadata">

**Author:** [@Theodore\_Rieger](https://discourse.julialang.org/u/Theodore_Rieger)\
**Replies:** 2\
**Last updated:** [March 9, 2023, 1:45am UTC](https://discourse.julialang.org/t/catalyst-jl-algebraically-relate-two-species/95762 "2023-03-09T01:45:40Z")

</div>

Using Catalyst.jl (DSL) to develop a model with pharmacology. The drug concentration is tracked in milligrams, but the pharmacology requires an algebraic transformation to ng/mL-unbound (multiplication by a few parameter…

---

## [JuMP and Ipopt optimization: How to reduce run time for iterative optimization?](https://discourse.julialang.org/t/jump-and-ipopt-optimization-how-to-reduce-run-time-for-iterative-optimization/95684)

<div class="topic-metadata">

**Author:** [@codingfriend](https://discourse.julialang.org/u/codingfriend)\
**Replies:** 3\
**Last updated:** [March 8, 2023, 6:48pm UTC](https://discourse.julialang.org/t/jump-and-ipopt-optimization-how-to-reduce-run-time-for-iterative-optimization/95684 "2023-03-08T18:48:53Z")

</div>

Hi everyone! I am seeking help to speed up my code using JuMP and Ipopt. My current approach is too slow for my needs. :grinning: My original code looks something like this: ‘’‘julia for i = 1:N optimizer = optimizer…

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

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