# Specific Domains

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

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

**Page:** 119

---

## [Addprocs: specify both multithreading and active environment](https://discourse.julialang.org/t/addprocs-specify-both-multithreading-and-active-environment/80427)

<div class="topic-metadata">

**Author:** [@Gregstrq](https://discourse.julialang.org/u/Gregstrq)\
**Replies:** 2\
**Last updated:** [April 16, 2024, 5:11pm UTC](https://discourse.julialang.org/t/addprocs-specify-both-multithreading-and-active-environment/80427 "2024-04-16T17:11:02Z")

</div>

Let’s say I want to add a process with 8 threads. This is easy: using Distributed addprocs(1, exeflags="-t 8") One can check that @fetchfrom 2 Threads.nthreads() returns 8 as expected. But, if I try to specify the env…

---

## [GeometricMachineLearning.jl: Looking for a direct contact opportunity with persons involved](https://discourse.julialang.org/t/geometricmachinelearning-jl-looking-for-a-direct-contact-opportunity-with-persons-involved/113025)

<div class="topic-metadata">

**Author:** [@j\_u](https://discourse.julialang.org/u/j_u)\
**Replies:** 0\
**Last updated:** [April 16, 2024, 3:42pm UTC](https://discourse.julialang.org/t/geometricmachinelearning-jl-looking-for-a-direct-contact-opportunity-with-persons-involved/113025 "2024-04-16T15:42:48Z")

</div>

Hi, I am looking for a direct contact opportunity with persons involved in GeometricMachineLearning.jl. I joined Julia Community a few years ago thanks to a hobby project of mine about "Neural network analysis of uncerta…

---

## [How to specify xticks with min, max an arbitrary number?](https://discourse.julialang.org/t/how-to-specify-xticks-with-min-max-an-arbitrary-number/113021)

<div class="topic-metadata">

**Author:** [@vladdez](https://discourse.julialang.org/u/vladdez)\
**Replies:** 0\
**Last updated:** [April 16, 2024, 3:03pm UTC](https://discourse.julialang.org/t/how-to-specify-xticks-with-min-max-an-arbitrary-number/113021 "2024-04-16T15:03:13Z")

</div>

I want to specify ax.xticks with min of xvalue, 0.0, mean of xvalue, max of xvalue. How can I implement that?

---

## [GPU memory usage increasing on each epoch (Flux)](https://discourse.julialang.org/t/gpu-memory-usage-increasing-on-each-epoch-flux/112942)

<div class="topic-metadata">

**Author:** [@evanfields](https://discourse.julialang.org/u/evanfields)\
**Replies:** 5\
**Last updated:** [April 16, 2024, 2:16pm UTC](https://discourse.julialang.org/t/gpu-memory-usage-increasing-on-each-epoch-flux/112942 "2024-04-16T14:16:01Z")

</div>

I’m catching up on the last few years of Flux progress (thanks, Flux devs!) by training some toy models, and am surprised to find that GPU memory usage is increasing on each training epoch, even when all training data is…

---

## [Cubature of Singular Integrand](https://discourse.julialang.org/t/cubature-of-singular-integrand/113003)

<div class="topic-metadata">

**Author:** [@Domenico\_Lahaye](https://discourse.julialang.org/u/Domenico_Lahaye)\
**Replies:** 3\
**Last updated:** [April 16, 2024, 1:33pm UTC](https://discourse.julialang.org/t/cubature-of-singular-integrand/113003 "2024-04-16T13:33:06Z")

</div>

?quadgk provides as output !!! note Only one-dimensional integrals are provided by this function. For multi-dimensional integration (cubature), there are many different algorithms (often much better than simple nested 1…

---

## [Question regarding numerical stability among Julia versions](https://discourse.julialang.org/t/question-regarding-numerical-stability-among-julia-versions/112741)

<div class="topic-metadata">

**Author:** [@CodeLenz](https://discourse.julialang.org/u/CodeLenz)\
**Replies:** 16\
**Last updated:** [April 16, 2024, 12:54pm UTC](https://discourse.julialang.org/t/question-regarding-numerical-stability-among-julia-versions/112741 "2024-04-16T12:54:05Z")

</div>

Dear Julians. I was testing the following code (M,C and K are sparse, symmetric and positive definite matrices) Julia\> M,C,K=THS(); julia\> cholM=cholesky(M); julia\> Cb = cholM\\C; julia\> Kb = cholM\\K; julia\> F = 0.…

---

## [Solving a nonlinear equation using NonlinearSolve.jl with AutoFiniteDiff() yields errors about dual number](https://discourse.julialang.org/t/solving-a-nonlinear-equation-using-nonlinearsolve-jl-with-autofinitediff-yields-errors-about-dual-number/112988)

<div class="topic-metadata">

**Author:** [@LINGXIAO\_QIN](https://discourse.julialang.org/u/LINGXIAO_QIN)\
**Replies:** 2\
**Last updated:** [April 16, 2024, 12:13pm UTC](https://discourse.julialang.org/t/solving-a-nonlinear-equation-using-nonlinearsolve-jl-with-autofinitediff-yields-errors-about-dual-number/112988 "2024-04-16T12:13:18Z")

</div>

I am trying to solve a nonlinear equation, which is constructed using the thermophysical properties of fluids (via CoolProp package). An MWE is shown below. using NonlinearSolve, CoolProp function myfun\_P\_MWE(p, (T, i…

---

## [Speed up JuMP.jl Estimation When Jacobian and Hessian of Constraints are Sparse](https://discourse.julialang.org/t/speed-up-jump-jl-estimation-when-jacobian-and-hessian-of-constraints-are-sparse/112981)

<div class="topic-metadata">

**Author:** [@Strange\_Xue](https://discourse.julialang.org/u/Strange_Xue)\
**Replies:** 6\
**Last updated:** [April 16, 2024, 11:08am UTC](https://discourse.julialang.org/t/speed-up-jump-jl-estimation-when-jacobian-and-hessian-of-constraints-are-sparse/112981 "2024-04-16T11:08:20Z")

</div>

Hi guys, in last questions (link is shown below), I asked how to supply gradient of objective and jacobian of constraints in nonlinear optimization. But @odow told me I dont have to do that because JuMP.jl will automatic…

---

## [NeuralODE with covariates](https://discourse.julialang.org/t/neuralode-with-covariates/112999)

<div class="topic-metadata">

**Author:** [@elcup](https://discourse.julialang.org/u/elcup)\
**Replies:** 1\
**Last updated:** [April 16, 2024, 9:34am UTC](https://discourse.julialang.org/t/neuralode-with-covariates/112999 "2024-04-16T09:34:16Z")

</div>

Hi. Am an quite new with Julia and I am stuck with something. Maybe you can help me (sorry if I am trying something stupid). I am working with Neural ODEs using DiffEqFlux. I want to model the differential equation with…

---

## [Error when using MKLDirectSolver in SCS](https://discourse.julialang.org/t/error-when-using-mkldirectsolver-in-scs/112982)

<div class="topic-metadata">

**Author:** [@charshaw](https://discourse.julialang.org/u/charshaw)\
**Replies:** 2\
**Last updated:** [April 15, 2024, 11:20pm UTC](https://discourse.julialang.org/t/error-when-using-mkldirectsolver-in-scs/112982 "2024-04-15T23:20:59Z")

</div>

I am using JuMP and calling SCS as a solver. I want to use SCS.MKLDirectSolver, as discussed in an earlier thread. As discussed on the SCS.jl README page, this requires using the package SCS\_MKL\_jll. But, I’m getting an…

---

## [JuMP.jl: Supplying gradient of objective and jacobian of constraints in Nonlinear Optimization](https://discourse.julialang.org/t/jump-jl-supplying-gradient-of-objective-and-jacobian-of-constraints-in-nonlinear-optimization/106940)

<div class="topic-metadata">

**Author:** [@Strange\_Xue](https://discourse.julialang.org/u/Strange_Xue)\
**Replies:** 12\
**Last updated:** [April 15, 2024, 9:14pm UTC](https://discourse.julialang.org/t/jump-jl-supplying-gradient-of-objective-and-jacobian-of-constraints-in-nonlinear-optimization/106940 "2024-04-15T21:14:31Z")

</div>

Hi, I am faced with a troublesome optimization problem. As shown below, the non-linear optimization problem is easy to model and solve via Ipopt with JuMP.jl as interface. However, the objective function and constraints…

---

## [ForwardDiff.jacobian Method error](https://discourse.julialang.org/t/forwarddiff-jacobian-method-error/112979)

<div class="topic-metadata">

**Author:** [@Domenico\_Lahaye](https://discourse.julialang.org/u/Domenico_Lahaye)\
**Replies:** 2\
**Last updated:** [April 15, 2024, 8:59pm UTC](https://discourse.julialang.org/t/forwarddiff-jacobian-method-error/112979 "2024-04-15T20:59:41Z")

</div>

I get a method error with ForwardDiff.jacobian that I fails to understand. Below a MWE. function fooworks(cvec,xmesh) residual = cvec return residual end function foofails(cvec,xmesh) residual = zeros(len…

---

## [Mesh is displayed properly on CairoMakie, but not in GLMakie](https://discourse.julialang.org/t/mesh-is-displayed-properly-on-cairomakie-but-not-in-glmakie/112977)

<div class="topic-metadata">

**Author:** [@marcsgil](https://discourse.julialang.org/u/marcsgil)\
**Replies:** 0\
**Last updated:** [April 15, 2024, 7:58pm UTC](https://discourse.julialang.org/t/mesh-is-displayed-properly-on-cairomakie-but-not-in-glmakie/112977 "2024-04-15T19:58:24Z")

</div>

Hello! I’m running the following code snippet to produce a sphere: using GLMakie fig = Figure() ax = Axis3(fig\[1, 1\], aspect=:equal) mesh!(ax, Sphere(Point3f(0), 1), alpha=0.1) scatter!(ax,0,0,0, color=:red) fig wi…

---

## [How to extract the upper/lower triangular matrix values except the diagonal?](https://discourse.julialang.org/t/how-to-extract-the-upper-lower-triangular-matrix-values-except-the-diagonal/112974)

<div class="topic-metadata">

**Author:** [@Shayan](https://discourse.julialang.org/u/Shayan)\
**Replies:** 4\
**Last updated:** [April 15, 2024, 6:07pm UTC](https://discourse.julialang.org/t/how-to-extract-the-upper-lower-triangular-matrix-values-except-the-diagonal/112974 "2024-04-15T18:07:52Z")

</div>

Consider the following matrix: julia\> c = rand(2, 2) 2×2 Matrix{Float64}: 0.0443768 0.515352 0.531426 0.967584 Expected output: 1-element Vector{Float64}: 0.531426 I can do: julia\> using LinearAlgebra: LowerTri…

---

## [Outputting string columns using Avro.jl](https://discourse.julialang.org/t/outputting-string-columns-using-avro-jl/80112)

<div class="topic-metadata">

**Author:** [@kobusherbst](https://discourse.julialang.org/u/kobusherbst)\
**Replies:** 5\
**Last updated:** [April 15, 2024, 5:13pm UTC](https://discourse.julialang.org/t/outputting-string-columns-using-avro-jl/80112 "2024-04-15T17:13:51Z")

</div>

When I try to convert an Arrow file to Avro using this code: function t\_arrowtoavro(path, file) a = Arrow.Table(joinpath(path, "$(file).arrow")) |\> DataFrame println(Tables.schema(a)) Avro.writetable(joinpat…

---

## [ODEProblem(....) vs NeuralODE(....) for neural ODEs](https://discourse.julialang.org/t/odeproblem-vs-neuralode-for-neural-odes/112670)

<div class="topic-metadata">

**Author:** [@KianH](https://discourse.julialang.org/u/KianH)\
**Replies:** 14\
**Last updated:** [April 15, 2024, 2:16pm UTC](https://discourse.julialang.org/t/odeproblem-vs-neuralode-for-neural-odes/112670 "2024-04-15T14:16:22Z")

</div>

Hello! I am wondering about the difference when using ODEProblem(....) from DifferentialEquations.jl and NeuralODE(....) from DiffEqFlux.jl in terms of time performance. Consider the code below, thus, definining the …

---

## [Sundials KLU Linear solver](https://discourse.julialang.org/t/sundials-klu-linear-solver/112954)

<div class="topic-metadata">

**Author:** [@jespfra](https://discourse.julialang.org/u/jespfra)\
**Replies:** 1\
**Last updated:** [April 15, 2024, 10:55am UTC](https://discourse.julialang.org/t/sundials-klu-linear-solver/112954 "2024-04-15T10:55:33Z")

</div>

Hi everyone, I am trying to use the KLU method as linear solver for solving an ODE problem using Sundials. I have tried to apply it to the CVODE\_BDF solver for a simple problem, the rober equation. However, when I prov…

---

## [Plot two legends (e.g. one for color and one for shape) with \`Plots.jl\`](https://discourse.julialang.org/t/plot-two-legends-e-g-one-for-color-and-one-for-shape-with-plots-jl/40746)

<div class="topic-metadata">

**Author:** [@Krastanov](https://discourse.julialang.org/u/Krastanov)\
**Replies:** 5\
**Last updated:** [April 15, 2024, 7:04am UTC](https://discourse.julialang.org/t/plot-two-legends-e-g-one-for-color-and-one-for-shape-with-plots-jl/40746 "2024-04-15T07:04:01Z")

</div>

How can I make a plot like this in Plots.jl (this was made with matplotlib): Basically four lines: plot( data1, color=c1, linestyle=s1) plot!(data2, color=c2, linestyle=s1) plot!(data3, color=c1, linestyle=s2) plot!(…

---

## [Solving shifted linear system](https://discourse.julialang.org/t/solving-shifted-linear-system/111265)

<div class="topic-metadata">

**Author:** [@rvaradhan](https://discourse.julialang.org/u/rvaradhan)\
**Replies:** 19\
**Last updated:** [April 15, 2024, 5:52am UTC](https://discourse.julialang.org/t/solving-shifted-linear-system/111265 "2024-04-15T05:52:22Z")

</div>

I would like to solve a large shifted linear system: I have to solve (A + k^2 I)x = y, where A is symmetric PD, y is a vector, k is a scalar, and I is identity matrix. Matrix A is large O(10000 x 10000) but it is banded…

---

## [Batched LU solves (or Factorizations) with Sparse Matrices](https://discourse.julialang.org/t/batched-lu-solves-or-factorizations-with-sparse-matrices/106019)

<div class="topic-metadata">

**Author:** [@schev](https://discourse.julialang.org/u/schev)\
**Replies:** 6\
**Last updated:** [April 15, 2024, 5:37am UTC](https://discourse.julialang.org/t/batched-lu-solves-or-factorizations-with-sparse-matrices/106019 "2024-04-15T05:37:51Z")

</div>

Let’s say I have a sparse set of matrices A1, A2, ... An, each of type CuSparseMatrixCSC{Float32, Int32} along with associated RHS vectors b1, b2, ... bn. Is there some watch to do a batch LU solve on the GPU with CUDA.j…

---

## [Updating CUDA seems to have completely broken the package](https://discourse.julialang.org/t/updating-cuda-seems-to-have-completely-broken-the-package/110101)

<div class="topic-metadata">

**Author:** [@Hareruya](https://discourse.julialang.org/u/Hareruya)\
**Replies:** 3\
**Last updated:** [April 14, 2024, 11:36pm UTC](https://discourse.julialang.org/t/updating-cuda-seems-to-have-completely-broken-the-package/110101 "2024-04-14T23:36:36Z")

</div>

Salutations, I am currently a PhD student and defended my thesis last month. Most of my Master’s research and all of my PhD research were developed using julia. I decided to update julia and run my old code at which poi…

---

## [Grouping by values in either of two columns](https://discourse.julialang.org/t/grouping-by-values-in-either-of-two-columns/112499)

<div class="topic-metadata">

**Author:** [@KeepLearning](https://discourse.julialang.org/u/KeepLearning)\
**Replies:** 13\
**Last updated:** [April 14, 2024, 8:47am UTC](https://discourse.julialang.org/t/grouping-by-values-in-either-of-two-columns/112499 "2024-04-14T08:47:12Z")

</div>

Hi, I wonder if there is an efficient way of doing the following in Julia DataFrames. Given a DataFrame df = DataFrame(col1 = \["a", "a", "b", "b", "c", "d"\], col2 = \["p", "q", "q", "r", "s", "t"\]) I wa…

---

## [Optimization with LBFGS gives DimensionMismatch("dimensions must match")](https://discourse.julialang.org/t/optimization-with-lbfgs-gives-dimensionmismatch-dimensions-must-match/22167)

<div class="topic-metadata">

**Author:** [@mzhenirovskyy](https://discourse.julialang.org/u/mzhenirovskyy)\
**Replies:** 9\
**Last updated:** [April 14, 2024, 7:39pm UTC](https://discourse.julialang.org/t/optimization-with-lbfgs-gives-dimensionmismatch-dimensions-must-match/22167 "2024-04-14T19:39:36Z")

</div>

My objective function has to solve an ODE. When I use the gradient-free algorithm (e.g NelderMead) to do optimization, everything is OK. But with gradient-based (e.g. LBFGS) I get the error: DimensionMismatch(“dimension…

---

## [Dynamical system of modelling two particle collisions](https://discourse.julialang.org/t/dynamical-system-of-modelling-two-particle-collisions/112922)

<div class="topic-metadata">

**Author:** [@hisacro](https://discourse.julialang.org/u/hisacro)\
**Replies:** 4\
**Last updated:** [April 14, 2024, 5:40pm UTC](https://discourse.julialang.org/t/dynamical-system-of-modelling-two-particle-collisions/112922 "2024-04-14T17:40:05Z")

</div>

Hello, I’m trying to model a system with collision dynamics. I’m well aware of Datseris’s DynamicalBillards, contrary to that what I’m studying is collision between 2 particles rather than the fixed geometry as in billa…

---

## [Setting up multi-output Gaussian process with Surrogates.jl](https://discourse.julialang.org/t/setting-up-multi-output-gaussian-process-with-surrogates-jl/112881)

<div class="topic-metadata">

**Author:** [@epi\_hannah](https://discourse.julialang.org/u/epi_hannah)\
**Replies:** 2\
**Last updated:** [April 14, 2024, 1:39pm UTC](https://discourse.julialang.org/t/setting-up-multi-output-gaussian-process-with-surrogates-jl/112881 "2024-04-14T13:39:51Z")

</div>

Hi all, I am trying to use Surrogates to train a Gaussian process surrogate to a model with 2 different outputs (i.e., a multi-output model). According to the Surrogates documentation I can use custom kriging in Stheno…

---

## [Using \`bboptimize\` with a function containing \`JuMP\` calls](https://discourse.julialang.org/t/using-bboptimize-with-a-function-containing-jump-calls/112846)

<div class="topic-metadata">

**Author:** [@Rick](https://discourse.julialang.org/u/Rick)\
**Replies:** 2\
**Last updated:** [April 14, 2024, 6:52am UTC](https://discourse.julialang.org/t/using-bboptimize-with-a-function-containing-jump-calls/112846 "2024-04-14T06:52:53Z")

</div>

Hi everyone, I am trying blackbox optimization with a function that contains JuMP functions. So, this is like a nested optimization where the “inner problem” is evaluated using JuMP functions and the “outer problem” is …

---

## [Gtk stopped working on MacOS in March or April 2024](https://discourse.julialang.org/t/gtk-stopped-working-on-macos-in-march-or-april-2024/112904)

<div class="topic-metadata">

**Author:** [@JeffFessler](https://discourse.julialang.org/u/JeffFessler)\
**Replies:** 2\
**Last updated:** [April 13, 2024, 4:38pm UTC](https://discourse.julialang.org/t/gtk-stopped-working-on-macos-in-march-or-april-2024/112904 "2024-04-13T16:38:41Z")

</div>

Gtk stopped working for me on MacOS sometime this semester (I am teaching a music signal processing course using Julia with \>60 students). Here is a MWE of code that worked a month or so ago and now does not work on Mac…

---

## [Optimization, from MATLAB to Julia](https://discourse.julialang.org/t/optimization-from-matlab-to-julia/111212)

<div class="topic-metadata">

**Author:** [@cu2mauro](https://discourse.julialang.org/u/cu2mauro)\
**Replies:** 10\
**Last updated:** [April 13, 2024, 11:43am UTC](https://discourse.julialang.org/t/optimization-from-matlab-to-julia/111212 "2024-04-13T11:43:23Z")

</div>

Hi there! I am fairly new to Julia, I was converted after seeing some awesome performance compared to MATLAB, which I currently use. I tried rewriting in Julia an optimization problem that I had in MATLAB, and everythin…

---

## [How to avoid "Julia is not responding" when using a file dialog](https://discourse.julialang.org/t/how-to-avoid-julia-is-not-responding-when-using-a-file-dialog/112075)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 5\
**Last updated:** [March 25, 2024, 12:14pm UTC](https://discourse.julialang.org/t/how-to-avoid-julia-is-not-responding-when-using-a-file-dialog/112075 "2024-03-25T12:14:47Z")

</div>

I have a GLMakie based GUI app and want to load a file using a native file dialog. This is my code: using GLMakie, NativeFileDialog fig = Figure(size=(200, 200), backgroundcolor=RGBf(0.7, 0.8, 1)) btn\_RUN = Button(fig…

---

## [Problem with hover display in heatmap with plotlyjs backend](https://discourse.julialang.org/t/problem-with-hover-display-in-heatmap-with-plotlyjs-backend/107586)

<div class="topic-metadata">

**Author:** [@adannenberg](https://discourse.julialang.org/u/adannenberg)\
**Replies:** 7\
**Last updated:** [April 13, 2024, 1:56am UTC](https://discourse.julialang.org/t/problem-with-hover-display-in-heatmap-with-plotlyjs-backend/107586 "2024-04-13T01:56:59Z")

</div>

if i use the plotly backend to make a heatmap with a small number of categorical variables then everything is fine using Plots plotly() categories = string.(collect('a':'e')) n = length(categories) #5 heatmap(categorie…

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

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