# Latest

**URL:** https://discourse.julialang.org/latest.md?page=392

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

**Page:** 393

---

## [Sampling predictive with MvNormal in Turing](https://discourse.julialang.org/t/sampling-predictive-with-mvnormal-in-turing/115292)

<div class="topic-metadata">

**Author:** [@dae](https://discourse.julialang.org/u/dae)\
**Replies:** 2\
**Last updated:** [June 7, 2024, 1:26am UTC](https://discourse.julialang.org/t/sampling-predictive-with-mvnormal-in-turing/115292 "2024-06-07T01:26:02Z")

</div>

tldr; none of these seem to work, and give different errors: Vector{Float64}(undef, n) ~ MvNormal(yhat, Diagonal(fill(1.0, n))) Vector{Real}(undef, n) ~ MvNormal(yhat, Diagonal(fill(1.0, n))) Vector{Missing}(missing, n)…

---

## [VSCode Debugger not hitting breakpoint in include file](https://discourse.julialang.org/t/vscode-debugger-not-hitting-breakpoint-in-include-file/46014)

<div class="topic-metadata">

**Author:** [@Brad\_Carman](https://discourse.julialang.org/u/Brad_Carman)\
**Replies:** 6\
**Last updated:** [June 7, 2024, 1:05am UTC](https://discourse.julialang.org/t/vscode-debugger-not-hitting-breakpoint-in-include-file/46014 "2024-06-07T01:05:15Z")

</div>

If I set a breakpoint in an included file, the breakpoint is not getting hit. Below shows my setup. I simply click the debug button to run test1.jl, the debugger runs and never stops at the break point. Running Julia …

---

## [JuMP write\_to\_file: Writing legacy model to file is much smaller](https://discourse.julialang.org/t/jump-write-to-file-writing-legacy-model-to-file-is-much-smaller/115299)

<div class="topic-metadata">

**Author:** [@Stuart\_Rogers](https://discourse.julialang.org/u/Stuart_Rogers)\
**Replies:** 7\
**Last updated:** [June 7, 2024, 12:57am UTC](https://discourse.julialang.org/t/jump-write-to-file-writing-legacy-model-to-file-is-much-smaller/115299 "2024-06-07T00:57:40Z")

</div>

I construct a model for a nonlinear optimization problem in the Julia JuMP code below, using either the non-legacy (use\_legacy\_model = false) or legacy (use\_legacy\_model = true) model interface and save the model to a .n…

---

## [How to determine which variables are fixed to 0 or 1 in a node of a branch and bound tree using CPLEX.jl?](https://discourse.julialang.org/t/how-to-determine-which-variables-are-fixed-to-0-or-1-in-a-node-of-a-branch-and-bound-tree-using-cplex-jl/115028)

<div class="topic-metadata">

**Author:** [@ASUKaiwenFang](https://discourse.julialang.org/u/ASUKaiwenFang)\
**Replies:** 10\
**Last updated:** [June 7, 2024, 12:56am UTC](https://discourse.julialang.org/t/how-to-determine-which-variables-are-fixed-to-0-or-1-in-a-node-of-a-branch-and-bound-tree-using-cplex-jl/115028 "2024-06-07T00:56:04Z")

</div>

For example, at node A, x\_1 = 0, x\_2 = 1, x\_3 = 1, which API should I use to query this information?

---

## [A good way of running Julia MPI jobs with slurm on HPC platforms](https://discourse.julialang.org/t/a-good-way-of-running-julia-mpi-jobs-with-slurm-on-hpc-platforms/115244)

<div class="topic-metadata">

**Author:** [@PetrKryslUCSD](https://discourse.julialang.org/u/PetrKryslUCSD)\
**Replies:** 7\
**Last updated:** [June 6, 2024, 4:53pm UTC](https://discourse.julialang.org/t/a-good-way-of-running-julia-mpi-jobs-with-slurm-on-hpc-platforms/115244 "2024-06-06T16:53:07Z")

</div>

For some reason the Julia MPI jobs that I submit end in error due to a lack of instantiation / activation. I do run it as mpiexec julia --project=~/a64fx/FinEtoolsDDParallel.jl/examples ~/a64fx/FinEtoolsDDParallel.jl/e…

---

## [Is there a way to change the theta tick label rotation for a PolarAxis? (Makie.jl)](https://discourse.julialang.org/t/is-there-a-way-to-change-the-theta-tick-label-rotation-for-a-polaraxis-makie-jl/115303)

<div class="topic-metadata">

**Author:** [@BenjaminDoran](https://discourse.julialang.org/u/BenjaminDoran)\
**Replies:** 0\
**Last updated:** [June 6, 2024, 11:14pm UTC](https://discourse.julialang.org/t/is-there-a-way-to-change-the-theta-tick-label-rotation-for-a-polaraxis-makie-jl/115303 "2024-06-06T23:14:15Z")

</div>

It looks like there is no thetaticklabelrotation option for PolarAxis is there another option to get rotated tick labels for polar axis plots?

---

## [NLP: Attribute ScalarNonlinearFunction is not supported by the model](https://discourse.julialang.org/t/nlp-attribute-scalarnonlinearfunction-is-not-supported-by-the-model/110444)

<div class="topic-metadata">

**Author:** [@Shuhua](https://discourse.julialang.org/u/Shuhua)\
**Replies:** 21\
**Last updated:** [June 6, 2024, 10:13pm UTC](https://discourse.julialang.org/t/nlp-attribute-scalarnonlinearfunction-is-not-supported-by-the-model/110444 "2024-06-06T22:13:50Z")

</div>

Hi, I am working on an optimization problem as described in the PDF document in the GitHub repository. The nonlinearity is introduced mainly the equation G in (16). I implemented the model following closely the PDF desc…

---

## [Cost model of weak references](https://discourse.julialang.org/t/cost-model-of-weak-references/115207)

<div class="topic-metadata">

**Author:** [@cstjean](https://discourse.julialang.org/u/cstjean)\
**Replies:** 1\
**Last updated:** [June 6, 2024, 9:49pm UTC](https://discourse.julialang.org/t/cost-model-of-weak-references/115207 "2024-06-06T21:49:21Z")

</div>

A long time ago, on another system (SBCL), I’ve been bitten by weak references having a quadratic O(n²) GC cost. I’m curious about Julia. Is there any admonition against using WeakRef and WeakDict a lot? The sweep phase…

---

## [\[ANN\] CliffordNumbers.jl: fast, static geometric algebra with multivectors as Number subtypes](https://discourse.julialang.org/t/ann-cliffordnumbers-jl-fast-static-geometric-algebra-with-multivectors-as-number-subtypes/114073)

<div class="topic-metadata">

**Author:** [@brainandforce](https://discourse.julialang.org/u/brainandforce)\
**Replies:** 11\
**Last updated:** [June 6, 2024, 9:33pm UTC](https://discourse.julialang.org/t/ann-cliffordnumbers-jl-fast-static-geometric-algebra-with-multivectors-as-number-subtypes/114073 "2024-06-06T21:33:01Z")

</div>

CliffordNumbers.jl is a new package for working with Clifford algebras (or geometric algebras) and their elements. Design philosophy AbstractCliffordNumber{Q,T\<:Union{Real,Complex}} type subtypes Number, just as Complex …

---

## [Help resolving perturbation confusion in \`ForwardDiff.jl\` with manual higher-order derivatives](https://discourse.julialang.org/t/help-resolving-perturbation-confusion-in-forwarddiff-jl-with-manual-higher-order-derivatives/115297)

<div class="topic-metadata">

**Author:** [@cgeoga](https://discourse.julialang.org/u/cgeoga)\
**Replies:** 0\
**Last updated:** [June 6, 2024, 8:25pm UTC](https://discourse.julialang.org/t/help-resolving-perturbation-confusion-in-forwarddiff-jl-with-manual-higher-order-derivatives/115297 "2024-06-06T20:25:39Z")

</div>

Hey ForwardDiffers, I’m wondering if I could get some help with resolving a perturbation confusing issue I’m having with some code. I use the follow setup to compute higher-order derivatives with ForwardDiff in a fast a…

---

## [Lock free sorted linked list](https://discourse.julialang.org/t/lock-free-sorted-linked-list/115157)

<div class="topic-metadata">

**Author:** [@melevy](https://discourse.julialang.org/u/melevy)\
**Replies:** 18\
**Last updated:** [June 6, 2024, 7:29pm UTC](https://discourse.julialang.org/t/lock-free-sorted-linked-list/115157 "2024-06-06T19:29:10Z")

</div>

Hi! I’m trying to implement a lock free sorted linked list as shown in 9.8 Nonblocking synchronization of the book titled The art of multiprocessor programming. I need something like AtomicMarkableReference or maybe…

---

## [Trouble with CPU+GPU adjoints for UDEs](https://discourse.julialang.org/t/trouble-with-cpu-gpu-adjoints-for-udes/115290)

<div class="topic-metadata">

**Author:** [@19fitzna](https://discourse.julialang.org/u/19fitzna)\
**Replies:** 1\
**Last updated:** [June 6, 2024, 5:49pm UTC](https://discourse.julialang.org/t/trouble-with-cpu-gpu-adjoints-for-udes/115290 "2024-06-06T17:49:02Z")

</div>

Hi all! I am currently trying to create UDEs using a GPU backend for the neural network only. I have been following the examples here, but have not been able to get anything to work. I was wondering if anyone might know …

---

## [Slicing a CuArray in a kernel](https://discourse.julialang.org/t/slicing-a-cuarray-in-a-kernel/115285)

<div class="topic-metadata">

**Author:** [@peteole](https://discourse.julialang.org/u/peteole)\
**Replies:** 1\
**Last updated:** [June 6, 2024, 4:21pm UTC](https://discourse.julialang.org/t/slicing-a-cuarray-in-a-kernel/115285 "2024-06-06T16:21:35Z")

</div>

Hi, I would like to take a slice from a CuArray in a kernel: using CUDA using StaticArrays using LinearAlgebra function run() a::Float64=5.0 n=3 vectors = CUDA.rand(Float64, 100,n) #vectors=CuArray(\[@S…

---

## [Solving linear system BC = X where B is a symmetric positive semi-definite matrix](https://discourse.julialang.org/t/solving-linear-system-bc-x-where-b-is-a-symmetric-positive-semi-definite-matrix/115280)

<div class="topic-metadata">

**Author:** [@Uranium238](https://discourse.julialang.org/u/Uranium238)\
**Replies:** 7\
**Last updated:** [June 6, 2024, 4:19pm UTC](https://discourse.julialang.org/t/solving-linear-system-bc-x-where-b-is-a-symmetric-positive-semi-definite-matrix/115280 "2024-06-06T16:19:18Z")

</div>

I am trying to solve for X in BC=X where B is a symmetric positive semi-definite matrix and X is simply a vector with 1 in all of it’s entries. As far as I understand, this can be solved using C=B\\X , but I was wondering…

---

## [Workarounds for splatting in index sets?](https://discourse.julialang.org/t/workarounds-for-splatting-in-index-sets/115212)

<div class="topic-metadata">

**Author:** [@solliolli](https://discourse.julialang.org/u/solliolli)\
**Replies:** 10\
**Last updated:** [June 6, 2024, 7:19am UTC](https://discourse.julialang.org/t/workarounds-for-splatting-in-index-sets/115212 "2024-06-06T07:19:40Z")

</div>

Hi! I ran into a small issue when writing a MILP model, namely the fact that the splatting operator cannot be used in the definition of an index set. What I’m trying to achieve is to define JuMP variables with dimensions…

---

## [Constructing a diagonal 2D CuArray while preserving Zygote's gradient functionality on the diagonal line](https://discourse.julialang.org/t/constructing-a-diagonal-2d-cuarray-while-preserving-zygotes-gradient-functionality-on-the-diagonal-line/115246)

<div class="topic-metadata">

**Author:** [@hxjz233](https://discourse.julialang.org/u/hxjz233)\
**Replies:** 3\
**Last updated:** [June 6, 2024, 3:12pm UTC](https://discourse.julialang.org/t/constructing-a-diagonal-2d-cuarray-while-preserving-zygotes-gradient-functionality-on-the-diagonal-line/115246 "2024-06-06T15:12:42Z")

</div>

I am a beginner in the use of CUDA and Zygote, and is currently carrying my code from a CPU to a GPU version. In my calculation, I need to put a 1D CuArray on the diagonal line of a 2D CuArray, and do some matrix multip…

---

## [Mul! with 3-dimensional array](https://discourse.julialang.org/t/mul-with-3-dimensional-array/115263)

<div class="topic-metadata">

**Author:** [@Sihyun\_Kim](https://discourse.julialang.org/u/Sihyun_Kim)\
**Replies:** 3\
**Last updated:** [June 6, 2024, 3:08pm UTC](https://discourse.julialang.org/t/mul-with-3-dimensional-array/115263 "2024-06-06T15:08:09Z")

</div>

I am new to Julia. I need to multiply 3 dimensional array by 2 dimensional array. The simplified form is here: vR = fill(1.0, (251, 51, 21)); Π = fill(1.0, (21, 21)); function mat\_mul(A, B) C = similar(A) @inboun…

---

## [Logging from a Pluto notebook into the REPL from which it launched](https://discourse.julialang.org/t/logging-from-a-pluto-notebook-into-the-repl-from-which-it-launched/115269)

<div class="topic-metadata">

**Author:** [@gdalle](https://discourse.julialang.org/u/gdalle)\
**Replies:** 4\
**Last updated:** [June 6, 2024, 2:42pm UTC](https://discourse.julialang.org/t/logging-from-a-pluto-notebook-into-the-repl-from-which-it-launched/115269 "2024-06-06T14:42:01Z")

</div>

Is there a way for a Pluto notebook to print/log into the REPL from which I launched it? I would need that to track progress of notebook CI runs

---

## [Computing fidelity of quantum state](https://discourse.julialang.org/t/computing-fidelity-of-quantum-state/115264)

<div class="topic-metadata">

**Author:** [@fqguo](https://discourse.julialang.org/u/fqguo)\
**Replies:** 6\
**Last updated:** [June 6, 2024, 2:24pm UTC](https://discourse.julialang.org/t/computing-fidelity-of-quantum-state/115264 "2024-06-06T14:24:15Z")

</div>

I write my program, however, i cannot plot the correct figure. When i use the ket psi\_0 as initial state, i can plot the correct figure same Matlab. However, when i use the density psi\_0⊗dagger(psi\_0), i cannot plot the…

---

## [Pkg Extension loading errors only on HPC filesystem](https://discourse.julialang.org/t/pkg-extension-loading-errors-only-on-hpc-filesystem/115035)

<div class="topic-metadata">

**Author:** [@jling](https://discourse.julialang.org/u/jling)\
**Replies:** 3\
**Last updated:** [June 6, 2024, 2:21pm UTC](https://discourse.julialang.org/t/pkg-extension-loading-errors-only-on-hpc-filesystem/115035 "2024-06-06T14:21:40Z")

</div>

Status \`/tmp/jiling/jl\_ZPW1yy/Project.toml\` \[8197267c\] IntervalSets v0.7.10 \[9a3f8284\] Random julia\> using IntervalSets julia\> Base.retry\_load\_extensions() ┌ Error: Error during loading of extension IntervalSetsRan…

---

## [How to import submodules](https://discourse.julialang.org/t/how-to-import-submodules/115277)

<div class="topic-metadata">

**Author:** [@The\_Mastermage](https://discourse.julialang.org/u/The_Mastermage)\
**Replies:** 1\
**Last updated:** [June 6, 2024, 2:20pm UTC](https://discourse.julialang.org/t/how-to-import-submodules/115277 "2024-06-06T14:20:08Z")

</div>

So I want to write a package for Querying Astronomy data and I want to make extensive use of submodules to seperate different queryable services. An example is in the following code: module AstronomyQuery module Sim…

---

## [How to disable automatic precompilation of packages in Julia 1.6](https://discourse.julialang.org/t/how-to-disable-automatic-precompilation-of-packages-in-julia-1-6/57053)

<div class="topic-metadata">

**Author:** [@amrods](https://discourse.julialang.org/u/amrods)\
**Replies:** 2\
**Last updated:** [June 6, 2024, 1:58pm UTC](https://discourse.julialang.org/t/how-to-disable-automatic-precompilation-of-packages-in-julia-1-6/57053 "2024-06-06T13:58:44Z")

</div>

How can I prevent automatic precompilation after adding a package in Julia 1.6?

---

## [Memory buffer type?](https://discourse.julialang.org/t/memory-buffer-type/15345)

<div class="topic-metadata">

**Author:** [@loisel](https://discourse.julialang.org/u/loisel)\
**Replies:** 9\
**Last updated:** [June 6, 2024, 1:16pm UTC](https://discourse.julialang.org/t/memory-buffer-type/15345 "2024-06-06T13:16:54Z")

</div>

I thought I saw somewhere that multidimensional arrays ought to be built on top of arrays of bytes. This would be useful to me for a thing. Is this happening?

---

## [Best practice to choose epoch or number of features from a recursive feature elimination process?](https://discourse.julialang.org/t/best-practice-to-choose-epoch-or-number-of-features-from-a-recursive-feature-elimination-process/115281)

<div class="topic-metadata">

**Author:** [@liuyxpp](https://discourse.julialang.org/u/liuyxpp)\
**Replies:** 0\
**Last updated:** [June 6, 2024, 1:13pm UTC](https://discourse.julialang.org/t/best-practice-to-choose-epoch-or-number-of-features-from-a-recursive-feature-elimination-process/115281 "2024-06-06T13:13:23Z")

</div>

In general, how to find the smallest number which corresponds to best score in a noisy curve? Below is an example result form a recursive feature elimination process. The x axis is the number of features selected and the…

---

## [Pkg behaving differently for Package from Registry vs the direct URL](https://discourse.julialang.org/t/pkg-behaving-differently-for-package-from-registry-vs-the-direct-url/115228)

<div class="topic-metadata">

**Author:** [@Simon\_Exner](https://discourse.julialang.org/u/Simon_Exner)\
**Replies:** 2\
**Last updated:** [June 6, 2024, 12:50pm UTC](https://discourse.julialang.org/t/pkg-behaving-differently-for-package-from-registry-vs-the-direct-url/115228 "2024-06-06T12:50:48Z")

</div>

Hello, I have an issue with the behavior of pkg: add FMIZoo results in FMIZoo.jl beeing installed and along with it Interpolations.jl 0.14.0 (which is in FMIZoo.jls dependencys) add "https://github.com/ThummeTo/FMIZoo…

---

## [Maximum of a continous function](https://discourse.julialang.org/t/maximum-of-a-continous-function/115258)

<div class="topic-metadata">

**Author:** [@cmichelenstrofer](https://discourse.julialang.org/u/cmichelenstrofer)\
**Replies:** 2\
**Last updated:** [June 6, 2024, 12:16pm UTC](https://discourse.julialang.org/t/maximum-of-a-continous-function/115258 "2024-06-06T12:16:48Z")

</div>

Packages like Plots and Integrals are able to determine point spacing by themselves (I imagine it is some sort of adaptive algorithm based on how quickly the function is changing). I would imagine that something similar …

---

## [LoadError: Cannot locate artifact 'kaleido' for aarch64-apple-darwin…](https://discourse.julialang.org/t/loaderror-cannot-locate-artifact-kaleido-for-aarch64-apple-darwin/115245)

<div class="topic-metadata">

**Author:** [@artbx25](https://discourse.julialang.org/u/artbx25)\
**Replies:** 1\
**Last updated:** [June 6, 2024, 11:19am UTC](https://discourse.julialang.org/t/loaderror-cannot-locate-artifact-kaleido-for-aarch64-apple-darwin/115245 "2024-06-06T11:19:41Z")

</div>

Hi, I’m new to Julia and trying to run a project on a mac M1 with Julia v1.10 that needs the packages: ✗ PlotlyBase ✗ ORCA ✗ PlotlyJS However when trying to add them I get an error (they also can’t be precompiled). …

---

## [TODAY JuliaUserGroupMunich: Jay Pottharst - Probabilistic Computation and Gen Primer](https://discourse.julialang.org/t/today-juliausergroupmunich-jay-pottharst-probabilistic-computation-and-gen-primer/115101)

<div class="topic-metadata">

**Author:** [@schlichtanders](https://discourse.julialang.org/u/schlichtanders)\
**Replies:** 3\
**Last updated:** [June 6, 2024, 9:32am UTC](https://discourse.julialang.org/t/today-juliausergroupmunich-jay-pottharst-probabilistic-computation-and-gen-primer/115101 "2024-06-06T09:32:53Z")

</div>

:round\_pushpin:online, organized via meetup :spiral\_calendar: TODAY, Monday 3. June 2024 :clock630: CEST 18:30 Probabilistic Programming with Gen.jl ▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔▔ Finally again some probabilities! Julia has a f…

---

## [ANN: JuMP-dev 2024: tickets now on sale](https://discourse.julialang.org/t/ann-jump-dev-2024-tickets-now-on-sale/111168)

<div class="topic-metadata">

**Author:** [@odow](https://discourse.julialang.org/u/odow)\
**Replies:** 2\
**Last updated:** [June 6, 2024, 9:30am UTC](https://discourse.julialang.org/t/ann-jump-dev-2024-tickets-now-on-sale/111168 "2024-06-06T09:30:18Z")

</div>

Dear all, Tickets for JuMP-dev 2024 are now on sale. The early bird cost is CAD50 before June 19, and CAD75 after. To sign up directly, go to https://symposia.gerad.ca/JumpDev-2024/register. For more details, includi…

---

## [Intersection of level sets](https://discourse.julialang.org/t/intersection-of-level-sets/114913)

<div class="topic-metadata">

**Author:** [@cmarcotte](https://discourse.julialang.org/u/cmarcotte)\
**Replies:** 12\
**Last updated:** [June 6, 2024, 9:26am UTC](https://discourse.julialang.org/t/intersection-of-level-sets/114913 "2024-06-06T09:26:57Z")

</div>

Hello! I am computing the intersections of two level sets computed from some input arrays. The idea is that, having retrieved the contours from Makie we can compute the intersections of the curves by iterating over sequ…

[Previous page](https://discourse.julialang.org/latest.md?page=391)

[Next page](https://discourse.julialang.org/latest.md?page=393)