# Latest

**URL:** https://discourse.julialang.org/latest.md?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

---

## [Custom (NumPy style) broadcasting rule that avoids iterating over elements (for GPU-acceleration)](https://discourse.julialang.org/t/custom-numpy-style-broadcasting-rule-that-avoids-iterating-over-elements-for-gpu-acceleration/131515)

<div class="topic-metadata">

**Author:** [@TimHargreaves](https://discourse.julialang.org/u/TimHargreaves)\
**Replies:** 10\
**Last updated:** [August 24, 2025, 7:48pm UTC](https://discourse.julialang.org/t/custom-numpy-style-broadcasting-rule-that-avoids-iterating-over-elements-for-gpu-acceleration/131515 "2025-08-24T19:48:57Z")

</div>

Minimal Example As a minimal example, suppose that I have the following two structs defined, struct Foo{T, M\<:AbstractMatrix{T}} A::M B::M end struct Bar{T, M\<:AbstractMatrix{T}} C::M foo::Foo{T, M} end…

---

## [JuliaPackageComparisons is looking for a new maintainer](https://discourse.julialang.org/t/juliapackagecomparisons-is-looking-for-a-new-maintainer/130667)

<div class="topic-metadata">

**Author:** [@TheLateKronos](https://discourse.julialang.org/u/TheLateKronos)\
**Replies:** 11\
**Last updated:** [August 24, 2025, 3:48pm UTC](https://discourse.julialang.org/t/juliapackagecomparisons-is-looking-for-a-new-maintainer/130667 "2025-08-24T15:48:10Z")

</div>

Hi everyone! Introduction For those who have not heard about JuliaPackageComparisons (which I will reference as JPC), I made a pre-announcement about it here in september 2023. Since then a lot has happened, but not eno…

---

## [\[ANN\] StippleUsers](https://discourse.julialang.org/t/ann-stippleusers/131799)

<div class="topic-metadata">

**Author:** [@hhaensel](https://discourse.julialang.org/u/hhaensel)\
**Replies:** 1\
**Last updated:** [August 24, 2025, 7:22pm UTC](https://discourse.julialang.org/t/ann-stippleusers/131799 "2025-08-24T19:22:36Z")

</div>

I’m happy to announce StippleUsers.jl - a middle layer for user authorisation in the GenieFramework. What does it do? It predefines and exports functions get\_user(), get\_roles(), is\_authorised() which can be called in a…

---

## [Null initialisation for a loop that appends vectors](https://discourse.julialang.org/t/null-initialisation-for-a-loop-that-appends-vectors/131754)

<div class="topic-metadata">

**Author:** [@rcomm](https://discourse.julialang.org/u/rcomm)\
**Replies:** 14\
**Last updated:** [August 24, 2025, 5:35pm UTC](https://discourse.julialang.org/t/null-initialisation-for-a-loop-that-appends-vectors/131754 "2025-08-24T17:35:48Z")

</div>

In Matlab, it is common to initialise a loop like this: x = ; for i = 1:n; x = \[x; another numerical row vector (of fixed length L)\]; end This creates a matrix of size \[n,L\]. It works regardless of chosen length L. W…

---

## [Does \`@label\`/\`@goto\` introduce scope like loop blocks?](https://discourse.julialang.org/t/does-label-goto-introduce-scope-like-loop-blocks/129607)

<div class="topic-metadata">

**Author:** [@Benny](https://discourse.julialang.org/u/Benny)\
**Replies:** 5\
**Last updated:** [August 24, 2025, 4:17pm UTC](https://discourse.julialang.org/t/does-label-goto-introduce-scope-like-loop-blocks/129607 "2025-08-24T16:17:45Z")

</div>

Variables assigned in a loop block scope are different across iterations while if statements don’t introduce new local scope at all, but what’s the variable scoping behavior for @label-@goto jumps? Initially I’d assume t…

---

## [Close instead of wait if TCP socket is never drained](https://discourse.julialang.org/t/close-instead-of-wait-if-tcp-socket-is-never-drained/131727)

<div class="topic-metadata">

**Author:** [@ettersi](https://discourse.julialang.org/u/ettersi)\
**Replies:** 1\
**Last updated:** [August 24, 2025, 3:23pm UTC](https://discourse.julialang.org/t/close-instead-of-wait-if-tcp-socket-is-never-drained/131727 "2025-08-24T15:23:50Z")

</div>

I was wondering how the triplet of OS, libuv and Julia handles the scenario where a TCP socket is constantly written to but never read from on the peer’s side, so I came up with this little test script. s = listen(1234)…

---

## [Going off the rails wrangling Python's type system](https://discourse.julialang.org/t/going-off-the-rails-wrangling-pythons-type-system/131755)

<div class="topic-metadata">

**Author:** [@mrufsvold](https://discourse.julialang.org/u/mrufsvold)\
**Replies:** 5\
**Last updated:** [August 24, 2025, 12:52pm UTC](https://discourse.julialang.org/t/going-off-the-rails-wrangling-pythons-type-system/131755 "2025-08-24T12:52:04Z")

</div>

I love types… maybe too much. It’s definitely my favorite thing about Julia. I recently made the type checker much stricter on my Python projects trying to get better clarity about my Python code. But that means repeati…

---

## [PresetTimeCallback with array in affect function](https://discourse.julialang.org/t/presettimecallback-with-array-in-affect-function/131806)

<div class="topic-metadata">

**Author:** [@Marco\_Masi](https://discourse.julialang.org/u/Marco_Masi)\
**Replies:** 1\
**Last updated:** [August 24, 2025, 11:57am UTC](https://discourse.julialang.org/t/presettimecallback-with-array-in-affect-function/131806 "2025-08-24T11:57:08Z")

</div>

I’m trying to understand how to use the PresetTimeCallback(tstops,affect!) where tstops is an array of time steps and the affect function should also affect the integration each time with different values. For example, …

---

## [How to make a series optional in AlgebraOfGraphics.jl](https://discourse.julialang.org/t/how-to-make-a-series-optional-in-algebraofgraphics-jl/131798)

<div class="topic-metadata">

**Author:** [@Abhro](https://discourse.julialang.org/u/Abhro)\
**Replies:** 2\
**Last updated:** [August 24, 2025, 9:23am UTC](https://discourse.julialang.org/t/how-to-make-a-series-optional-in-algebraofgraphics-jl/131798 "2025-08-24T09:23:10Z")

</div>

I have a “wide” DataFrame where one column is a shared vector of x-values, and two different columns are separate vectors of y-values. For example, df = let t = range(0, 9, length = 1000) mass\_lab\_frame = @. cos…

---

## [MTK - solve: weird error message](https://discourse.julialang.org/t/mtk-solve-weird-error-message/131638)

<div class="topic-metadata">

**Author:** [@BLI](https://discourse.julialang.org/u/BLI)\
**Replies:** 7\
**Last updated:** [August 24, 2025, 7:28am UTC](https://discourse.julialang.org/t/mtk-solve-weird-error-message/131638 "2025-08-24T07:28:05Z")

</div>

I’m on Julia v1.11.6 on Windows 11 64bit, using MTK v10.19.0, OrdinaryDiffEq v6.101.0. I try to use Clapeyron v0.6.15 with the code. My code worked without a flaw (except some warnings related to arrays and Clapeyron) i…

---

## [Using functions loaded via pycall (kwargs)](https://discourse.julialang.org/t/using-functions-loaded-via-pycall-kwargs/131775)

<div class="topic-metadata">

**Author:** [@Geegee](https://discourse.julialang.org/u/Geegee)\
**Replies:** 10\
**Last updated:** [August 24, 2025, 5:36am UTC](https://discourse.julialang.org/t/using-functions-loaded-via-pycall-kwargs/131775 "2025-08-24T05:36:33Z")

</div>

I think this might be an easy thing but currently I don’t know how to do this. I managed to import my python code and can call it but part of the functions require me to do something like f(x=a, y=c, \*\*rest) and as Ju…

---

## [\[ANN\] SAMP.jl: Simple Application Messaging Protocol](https://discourse.julialang.org/t/ann-samp-jl-simple-application-messaging-protocol/131796)

<div class="topic-metadata">

**Author:** [@Marco\_Lombardi](https://discourse.julialang.org/u/Marco_Lombardi)\
**Replies:** 0\
**Last updated:** [August 23, 2025, 9:43pm UTC](https://discourse.julialang.org/t/ann-samp-jl-simple-application-messaging-protocol/131796 "2025-08-23T21:43:46Z")

</div>

I am happy to share a still unregistered package that provides a Julia implementation of the Simple Application Messaging Protocol (SAMP). SAMP is a messaging protocol that enables astronomy software tools to interopera…

---

## [Hidden version of searchsortedfirst](https://discourse.julialang.org/t/hidden-version-of-searchsortedfirst/131776)

<div class="topic-metadata">

**Author:** [@alan](https://discourse.julialang.org/u/alan)\
**Replies:** 4\
**Last updated:** [August 23, 2025, 9:05pm UTC](https://discourse.julialang.org/t/hidden-version-of-searchsortedfirst/131776 "2025-08-23T21:05:09Z")

</div>

In the Julia source code, there is an undocumented version of searchsortedfirst: function searchsortedfirst(v::AbstractVector, x, lo::T, hi::T, o::Ordering)::keytype(v) where T\<:Integer I need a function like this, in …

---

## [How to use Autodiff in Gridap.jl](https://discourse.julialang.org/t/how-to-use-autodiff-in-gridap-jl/131652)

<div class="topic-metadata">

**Author:** [@Sigmund](https://discourse.julialang.org/u/Sigmund)\
**Replies:** 4\
**Last updated:** [August 23, 2025, 5:45pm UTC](https://discourse.julialang.org/t/how-to-use-autodiff-in-gridap-jl/131652 "2025-08-23T17:45:57Z")

</div>

TLDR: I want to solve an optimal control problem using automatic differentiation (Gridap uses ForwardDiff.jl), but i get an error when i try to find the gradient of the objective function. The objective function depends …

---

## [I want advice - Struggling to Connect Python with Julia](https://discourse.julialang.org/t/i-want-advice-struggling-to-connect-python-with-julia/131788)

<div class="topic-metadata">

**Author:** [@Olivia](https://discourse.julialang.org/u/Olivia)\
**Replies:** 1\
**Last updated:** [August 23, 2025, 3:44pm UTC](https://discourse.julialang.org/t/i-want-advice-struggling-to-connect-python-with-julia/131788 "2025-08-23T15:44:52Z")

</div>

Hi everyone I have been experimenting with Julia for a while now mainly because I such as its speed & how well it handles numerical tasks. I wanted to integrate some of my existing Python scripts into Julia but I am faci…

---

## [Discrepancies of Errors caused by \`collect(Base.OneTo(N))\`](https://discourse.julialang.org/t/discrepancies-of-errors-caused-by-collect-base-oneto-n/131787)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 1\
**Last updated:** [August 23, 2025, 11:17am UTC](https://discourse.julialang.org/t/discrepancies-of-errors-caused-by-collect-base-oneto-n/131787 "2025-08-23T11:17:29Z")

</div>

I discover this (using Version 1.11.5 (2025-04-14) ) julia\> collect(Base.OneTo(2^59)) ERROR: OutOfMemoryError() Stacktrace: \[1\] GenericMemory @ .\\boot.jl:516 \[inlined\] \[2\] Array @ .\\boot.jl:578 \[inlined\] \[3\] Ar…

---

## [cuSOLVER: two calls to cusolverDnDgesvdj\_bufferSize, one via Juila, the other via CUDA yield (very) different results](https://discourse.julialang.org/t/cusolver-two-calls-to-cusolverdndgesvdj-buffersize-one-via-juila-the-other-via-cuda-yield-very-different-results/125904)

<div class="topic-metadata">

**Author:** [@FHulsemann](https://discourse.julialang.org/u/FHulsemann)\
**Replies:** 2\
**Last updated:** [August 22, 2025, 5:09pm UTC](https://discourse.julialang.org/t/cusolver-two-calls-to-cusolverdndgesvdj-buffersize-one-via-juila-the-other-via-cuda-yield-very-different-results/125904 "2025-08-22T17:09:21Z")

</div>

Hello. I was wondering whether someone could explain to me why the working space for an SVD (gesvdj to be precise) differs wildly between a call to cusolverDnDgesvdj\_bufferSize via the Julia interface and a call via CU…

---

## [3D isosurface plotting in Makie](https://discourse.julialang.org/t/3d-isosurface-plotting-in-makie/131470)

<div class="topic-metadata">

**Author:** [@mreichMPI-BGC](https://discourse.julialang.org/u/mreichMPI-BGC)\
**Replies:** 14\
**Last updated:** [August 22, 2025, 2:38pm UTC](https://discourse.julialang.org/t/3d-isosurface-plotting-in-makie/131470 "2025-08-22T14:38:11Z")

</div>

I am using the (semi-)documented contour function for 3D-isosurfaces, cf. volume | Makie, first example. One cannot define vector of alphas (which would be nice because outer surfaces (shells) “should” be more transpare…

---

## [Bring Intel x86 simd sort library to Julia](https://discourse.julialang.org/t/bring-intel-x86-simd-sort-library-to-julia/131719)

<div class="topic-metadata">

**Author:** [@gitboy16](https://discourse.julialang.org/u/gitboy16)\
**Replies:** 35\
**Last updated:** [August 22, 2025, 2:23pm UTC](https://discourse.julialang.org/t/bring-intel-x86-simd-sort-library-to-julia/131719 "2025-08-22T14:23:13Z")

</div>

Hi, Intel x86 simd sort library can be used via C++ template. If I wanted to make this library available to the maximum of Julia users (a bit like numpy did) what would be the best steps to follow? Should I use Binary…

---

## [Aqua Compare Project.toml and test/Project.toml fails](https://discourse.julialang.org/t/aqua-compare-project-toml-and-test-project-toml-fails/131760)

<div class="topic-metadata">

**Author:** [@Jake](https://discourse.julialang.org/u/Jake)\
**Replies:** 3\
**Last updated:** [August 22, 2025, 11:47am UTC](https://discourse.julialang.org/t/aqua-compare-project-toml-and-test-project-toml-fails/131760 "2025-08-22T11:47:16Z")

</div>

I am trying out Aqua on some packages and getting an error I cannot see to figure out. julia\> Aqua.test\_all(Winding;stale\_deps=true) Test Summary: | Pass Total Time Method ambiguity | 1 1 2.5s Test Summary…

---

## [\[Pre-ANN\] ARS.jl](https://discourse.julialang.org/t/pre-ann-ars-jl/131766)

<div class="topic-metadata">

**Author:** [@Eliassj](https://discourse.julialang.org/u/Eliassj)\
**Replies:** 4\
**Last updated:** [August 22, 2025, 10:31am UTC](https://discourse.julialang.org/t/pre-ann-ars-jl/131766 "2025-08-22T10:31:08Z")

</div>

Repo: GitHub - Eliassj/ARS.jl Docs: eliassj.github.io/ARS.jl/dev/ While working on another project I found that AdaptiveRejectionSampling.jl, while very good in other ways, did not allow for other numeric types than Fl…

---

## [Is there an "Algebra of basis-functions"](https://discourse.julialang.org/t/is-there-an-algebra-of-basis-functions/131750)

<div class="topic-metadata">

**Author:** [@JADekker](https://discourse.julialang.org/u/JADekker)\
**Replies:** 8\
**Last updated:** [August 22, 2025, 9:23am UTC](https://discourse.julialang.org/t/is-there-an-algebra-of-basis-functions/131750 "2025-08-22T09:23:50Z")

</div>

In my research, I often work with the following kind of problem that essentially reduces to approximating a conditional expectation, say \\mathbb{E}\[g(X)|X=x\], by a linear combination of basis-functions, \\mathbb{E}\[g(X)|X…

---

## [LocalPreferences not picked up in test when set globally](https://discourse.julialang.org/t/localpreferences-not-picked-up-in-test-when-set-globally/131748)

<div class="topic-metadata">

**Author:** [@luraess](https://discourse.julialang.org/u/luraess)\
**Replies:** 5\
**Last updated:** [August 22, 2025, 8:20am UTC](https://discourse.julialang.org/t/localpreferences-not-picked-up-in-test-when-set-globally/131748 "2025-08-22T08:20:08Z")

</div>

Hi, I am working with various packages such as Chmy.jl, JustRelax.jl, among others, where we support both CPU, GPU and distributed computing backends. Mostly for GPU backends and MPI support on servers such as HPC clus…

---

## [Parsing one line try catch](https://discourse.julialang.org/t/parsing-one-line-try-catch/131756)

<div class="topic-metadata">

**Author:** [@GHTaarn](https://discourse.julialang.org/u/GHTaarn)\
**Replies:** 5\
**Last updated:** [August 22, 2025, 4:03am UTC](https://discourse.julialang.org/t/parsing-one-line-try-catch/131756 "2025-08-22T04:03:01Z")

</div>

I was surprised that the following line does not parse: \_ \_ \_ \_(\_)\_ | Documentation: https://docs.julialang.org (\_) | (\_) (\_) | \_ \_ \_| |\_ \_\_ \_ | Type "?" for help, "\]?" for…

---

## [Finding package that requires admin permission to install](https://discourse.julialang.org/t/finding-package-that-requires-admin-permission-to-install/131758)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 3\
**Last updated:** [August 21, 2025, 8:01pm UTC](https://discourse.julialang.org/t/finding-package-that-requires-admin-permission-to-install/131758 "2025-08-21T20:01:51Z")

</div>

I gave a short-course this last weekend and the participants reported a pop-up window asking for admin permission while installing GeoStats.jl on Windows machines in their lab. How can I identify the indirect dependenc…

---

## [Efficient representation of a sparse "pentadiagonal" matrix?](https://discourse.julialang.org/t/efficient-representation-of-a-sparse-pentadiagonal-matrix/131744)

<div class="topic-metadata">

**Author:** [@bremez](https://discourse.julialang.org/u/bremez)\
**Replies:** 8\
**Last updated:** [August 21, 2025, 6:53pm UTC](https://discourse.julialang.org/t/efficient-representation-of-a-sparse-pentadiagonal-matrix/131744 "2025-08-21T18:53:13Z")

</div>

I am interested in the eigenspectrum of Hermitian matrices that have a “pentadiagonal” structure, with A\[m,n\]=0 if abs(m-n)\>2. I am most interested in eigenstates with values close to zero. At present I am using Arpack.j…

---

## [Putting non linear optimization in a real-time closed-loop](https://discourse.julialang.org/t/putting-non-linear-optimization-in-a-real-time-closed-loop/131753)

<div class="topic-metadata">

**Author:** [@daro.slai](https://discourse.julialang.org/u/daro.slai)\
**Replies:** 1\
**Last updated:** [August 21, 2025, 6:20pm UTC](https://discourse.julialang.org/t/putting-non-linear-optimization-in-a-real-time-closed-loop/131753 "2025-08-21T18:20:35Z")

</div>

Hi, I’m trying to implement an Energy Management Algorithm (QCP) I built using InfiniteOpt.jl. The algo gets measurements from a WebSocket connection. My main scripting function goes as: function oneStepOpt\_test1(set…

---

## [Pkg Manager Issues due to VPN/Windows 11?](https://discourse.julialang.org/t/pkg-manager-issues-due-to-vpn-windows-11/131633)

<div class="topic-metadata">

**Author:** [@co1emi11er](https://discourse.julialang.org/u/co1emi11er)\
**Replies:** 9\
**Last updated:** [August 21, 2025, 4:06pm UTC](https://discourse.julialang.org/t/pkg-manager-issues-due-to-vpn-windows-11/131633 "2025-08-21T16:06:10Z")

</div>

I use julia on a work computer and I don’t recall this being an issue in the past. Anytime I try to install or update packages, it pops up a message that looks something like this: (@v1.11) pkg\> registry up General …

---

## [Catalyst.jl - what's the best way to extract interpolation from the reaction\_network](https://discourse.julialang.org/t/catalyst-jl-whats-the-best-way-to-extract-interpolation-from-the-reaction-network/130171)

<div class="topic-metadata">

**Author:** [@Mindy](https://discourse.julialang.org/u/Mindy)\
**Replies:** 7\
**Last updated:** [August 21, 2025, 3:22pm UTC](https://discourse.julialang.org/t/catalyst-jl-whats-the-best-way-to-extract-interpolation-from-the-reaction-network/130171 "2025-08-21T15:22:11Z")

</div>

Hi! If I want to extract and plot the interpolating variable - $d\_rate (temporal dynamic profile) from this example code, does catalyst.jl have any built-in function so that I can extract the results directly, instead o…

---

## [Dagger not fully utilizing CPU cores](https://discourse.julialang.org/t/dagger-not-fully-utilizing-cpu-cores/131733)

<div class="topic-metadata">

**Author:** [@Miroboru](https://discourse.julialang.org/u/Miroboru)\
**Replies:** 5\
**Last updated:** [August 21, 2025, 1:05pm UTC](https://discourse.julialang.org/t/dagger-not-fully-utilizing-cpu-cores/131733 "2025-08-21T13:05:43Z")

</div>

I have a time consuming simulation that is run over a set of signal to noise ratios, and thus is trivial to parallelize. I am currently using Dagger.jl to start 24 processes (I am on a 24 core threadripper system with 25…

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

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