# Latest

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

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

**Page:** 452

---

## [\[ANN\] MultiPrecisionArrays.jl: Algorithms and data structures for iterative refinement](https://discourse.julialang.org/t/ann-multiprecisionarrays-jl-algorithms-and-data-structures-for-iterative-refinement/111993)

<div class="topic-metadata">

**Author:** [@ctkelley](https://discourse.julialang.org/u/ctkelley)\
**Replies:** 0\
**Last updated:** [March 22, 2024, 8:35pm UTC](https://discourse.julialang.org/t/ann-multiprecisionarrays-jl-algorithms-and-data-structures-for-iterative-refinement/111993 "2024-03-22T20:35:35Z")

</div>

I just tagged v0.1.0 of MultiPrecisionArrays.jl. This package provides data structures and solvers for several variants of iterative refinement (IR) for solving linear systems A x= b. IR reduces matrix factorization c…

---

## [Programmatically adding multiple colums to a dataframe](https://discourse.julialang.org/t/programmatically-adding-multiple-colums-to-a-dataframe/111983)

<div class="topic-metadata">

**Author:** [@pascal](https://discourse.julialang.org/u/pascal)\
**Replies:** 4\
**Last updated:** [March 22, 2024, 7:10pm UTC](https://discourse.julialang.org/t/programmatically-adding-multiple-colums-to-a-dataframe/111983 "2024-03-22T19:10:39Z")

</div>

Suppose I have the following code: n = \["a\_" \* string(i) for i in 1:3\] v = \[1,2,3\] df = DataFrame(a = \[1, 2\], b = \[3, 4\]) df\_new = @chain df begin @transform $(n\[1\]) = v\[1\] end Can I change the @transform to add 3 …

---

## [Does calling a function with no memory allocations create an allocation?](https://discourse.julialang.org/t/does-calling-a-function-with-no-memory-allocations-create-an-allocation/111987)

<div class="topic-metadata">

**Author:** [@OmerPaz](https://discourse.julialang.org/u/OmerPaz)\
**Replies:** 5\
**Last updated:** [March 22, 2024, 7:54pm UTC](https://discourse.julialang.org/t/does-calling-a-function-with-no-memory-allocations-create-an-allocation/111987 "2024-03-22T19:54:47Z")

</div>

Hi all, I am running a numerical simulation in which I document my results in a pre-allocated structure, Obs: mutable struct Obs{T} name::String data::Array{T} num\_of\_measure::Int64 c\_num\_of\_measure::In…

---

## [Are there any plans for Julia to adopt an effects system?](https://discourse.julialang.org/t/are-there-any-plans-for-julia-to-adopt-an-effects-system/110729)

<div class="topic-metadata">

**Author:** [@xiaoxi](https://discourse.julialang.org/u/xiaoxi)\
**Replies:** 7\
**Last updated:** [March 22, 2024, 6:54pm UTC](https://discourse.julialang.org/t/are-there-any-plans-for-julia-to-adopt-an-effects-system/110729 "2024-03-22T18:54:21Z")

</div>

It seems that industry programming languages such as Rust or Scala are studying the implementation of an effects system. Is there any Julep about it? For reference: Rust: Extending Rust's Effect System Scala: CanThr…

---

## [GLMakie crash when using PrecompileTools](https://discourse.julialang.org/t/glmakie-crash-when-using-precompiletools/111974)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 3\
**Last updated:** [March 22, 2024, 5:19pm UTC](https://discourse.julialang.org/t/glmakie-crash-when-using-precompiletools/111974 "2024-03-22T17:19:49Z")

</div>

I observe a strange crash when I execute in the terminal: mkdir tmp cd tmp julia --project="." and then in Julia using Pkg pkg"add KiteViewers" pkg"add KiteUtils" using KiteViewers, KiteUtils viewer=Viewer3D(true); s…

---

## [Looking for a quadrature rule](https://discourse.julialang.org/t/looking-for-a-quadrature-rule/111923)

<div class="topic-metadata">

**Author:** [@lrnv](https://discourse.julialang.org/u/lrnv)\
**Replies:** 17\
**Last updated:** [March 22, 2024, 4:52pm UTC](https://discourse.julialang.org/t/looking-for-a-quadrature-rule/111923 "2024-03-22T16:52:48Z")

</div>

Hi, I need compute thousands of integrals of the form \\int\_{0}^{t} \\frac{e^{at}}{\\sum\_{i=1}^n b\_i e^{c\_i t}} dt. I know that : the time frame \[0,t\] is typically small (say t \\le 1) rates a,c\_1,...c\_n are of order 1e…

---

## [Flux: Dimension mismatch error](https://discourse.julialang.org/t/flux-dimension-mismatch-error/111981)

<div class="topic-metadata">

**Author:** [@Longblackcoffee](https://discourse.julialang.org/u/Longblackcoffee)\
**Replies:** 1\
**Last updated:** [March 22, 2024, 4:57pm UTC](https://discourse.julialang.org/t/flux-dimension-mismatch-error/111981 "2024-03-22T16:57:09Z")

</div>

Hi everyone, I’m a beginner with Flux, and I want to use it to approximate a high-dimensional input function. However, I’m running into a dimension mismatch error. Here’s a MWE, where I’m just trying to approximate f(x,…

---

## [KernelAbstractions for splines](https://discourse.julialang.org/t/kernelabstractions-for-splines/104312)

<div class="topic-metadata">

**Author:** [@weymouth](https://discourse.julialang.org/u/weymouth)\
**Replies:** 6\
**Last updated:** [March 22, 2024, 2:58pm UTC](https://discourse.julialang.org/t/kernelabstractions-for-splines/104312 "2024-03-22T14:58:48Z")

</div>

I’m trying to code a simple b-spline (and then NURBS) on the GPU using kernels to speed up a ton of simultaneous evaluations. I’m having trouble because Cox-deBoor is recursive and throws an error when I try to use this…

---

## [VSCodeServer.VSCodeLogger freeze when changed logged and ripristinated](https://discourse.julialang.org/t/vscodeserver-vscodelogger-freeze-when-changed-logged-and-ripristinated/111870)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 9\
**Last updated:** [March 22, 2024, 1:53pm UTC](https://discourse.julialang.org/t/vscodeserver-vscodelogger-freeze-when-changed-logged-and-ripristinated/111870 "2024-03-22T13:53:10Z")

</div>

In a console this work as expected: julia\> using Logging julia\> old\_logger = Logging.current\_logger() ConsoleLogger(IOBuffer(data=UInt8\[...\], readable=false, writable=false, seekable=false, append=false, size=0, maxsiz…

---

## [Solving a Fokker-Plank PDE](https://discourse.julialang.org/t/solving-a-fokker-plank-pde/109101)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 5\
**Last updated:** [March 22, 2024, 1:45pm UTC](https://discourse.julialang.org/t/solving-a-fokker-plank-pde/109101 "2024-03-22T13:45:46Z")

</div>

I have a Wiener process in 3 dimensions, with time-independent drift and diffusion, ie x \\in \\mathbb{R}^3, dx\_t = \\mu(x) dt + \\sigma(x) dw where \\sigma(x) is diagonal. I also have an initial distribution/density p\_0(x…

---

## [Book: The Elements of Differentiable Programming](https://discourse.julialang.org/t/book-the-elements-of-differentiable-programming/111968)

<div class="topic-metadata">

**Author:** [@christianhvejsel](https://discourse.julialang.org/u/christianhvejsel)\
**Replies:** 1\
**Last updated:** [March 22, 2024, 1:04pm UTC](https://discourse.julialang.org/t/book-the-elements-of-differentiable-programming/111968 "2024-03-22T13:04:25Z")

</div>

Not Julia related per se, but a lot people in this community make use of and contribute to Differentiable Programming, here is a freely available book on the topic: \[2403.14606\] The Elements of Differentiable Programming …

---

## [Looking for Developer for Backtesting App](https://discourse.julialang.org/t/looking-for-developer-for-backtesting-app/111773)

<div class="topic-metadata">

**Author:** [@Kruxigt](https://discourse.julialang.org/u/Kruxigt)\
**Replies:** 1\
**Last updated:** [March 22, 2024, 1:00pm UTC](https://discourse.julialang.org/t/looking-for-developer-for-backtesting-app/111773 "2024-03-22T13:00:29Z")

</div>

Hi! Last few years I have developed a backtesting framework in Julia. It’s now in its alpha stage. The core of the framework is closed source and the aim is to monetize it at a future stage. I’m now looking for a partn…

---

## [Bayesian Design of Experiment in Julia](https://discourse.julialang.org/t/bayesian-design-of-experiment-in-julia/106372)

<div class="topic-metadata">

**Author:** [@Dreycen\_Foiles](https://discourse.julialang.org/u/Dreycen_Foiles)\
**Replies:** 1\
**Last updated:** [March 22, 2024, 1:00pm UTC](https://discourse.julialang.org/t/bayesian-design-of-experiment-in-julia/106372 "2024-03-22T13:00:22Z")

</div>

I’m learning more about Bayesian experimental design and I’d be interested to know if anyone has experience doing something like this in Julia. Specifically, I’m looking for anything like this Pyro tutorial on an adapati…

---

## [How is division by a singular matrix performed in Julia even it is theoretically looks invalid?](https://discourse.julialang.org/t/how-is-division-by-a-singular-matrix-performed-in-julia-even-it-is-theoretically-looks-invalid/111961)

<div class="topic-metadata">

**Author:** [@shubh\_b](https://discourse.julialang.org/u/shubh_b)\
**Replies:** 3\
**Last updated:** [March 22, 2024, 12:43pm UTC](https://discourse.julialang.org/t/how-is-division-by-a-singular-matrix-performed-in-julia-even-it-is-theoretically-looks-invalid/111961 "2024-03-22T12:43:48Z")

</div>

These are two matrices A1 and A2: julia\> A1 = \[2 4 6; 5 8 3; 2 1 8; 7 6 3\] 4×3 Matrix{Int64}: 2 4 6 5 8 3 2 1 8 7 6 3 julia\> A2 = reshape(collect(11:19), (3, 3)) 3×3 Matrix{Int64}: 11 14 17 12 15 18 1…

---

## [Indexing of vector of arrays](https://discourse.julialang.org/t/indexing-of-vector-of-arrays/111958)

<div class="topic-metadata">

**Author:** [@TongTongYee](https://discourse.julialang.org/u/TongTongYee)\
**Replies:** 2\
**Last updated:** [March 22, 2024, 11:35am UTC](https://discourse.julialang.org/t/indexing-of-vector-of-arrays/111958 "2024-03-22T11:35:22Z")

</div>

Hi! How do I index vector vector of arrays? Of course, I can easily index them in some cases. But I don’t know how to index them in specific cases. For example, A = Vector{Vector{Matrix{Int64}}}() for i = 1:10 pu…

---

## [Help understand boxing](https://discourse.julialang.org/t/help-understand-boxing/111829)

<div class="topic-metadata">

**Author:** [@jlamb](https://discourse.julialang.org/u/jlamb)\
**Replies:** 17\
**Last updated:** [March 22, 2024, 10:21am UTC](https://discourse.julialang.org/t/help-understand-boxing/111829 "2024-03-22T10:21:37Z")

</div>

I am trying to resolve a type stability problem in a function where some local variables are showing up as the type Core.Box when I run @code\_warntype. I’ve narrowed the cause to where I use a comprehension or generator …

---

## [Julia community in Lund](https://discourse.julialang.org/t/julia-community-in-lund/111825)

<div class="topic-metadata">

**Author:** [@RahulManavalan](https://discourse.julialang.org/u/RahulManavalan)\
**Replies:** 4\
**Last updated:** [March 22, 2024, 9:36am UTC](https://discourse.julialang.org/t/julia-community-in-lund/111825 "2024-03-22T09:36:23Z")

</div>

Hi there I am new in Lund, Sweden; starting a PhD in Mathematical Statistics. I am curious if there exists a Julia community in Lund. Thanks in advance.

---

## [How to check the dimension names before eliminating some dimensions of a polyhedron?](https://discourse.julialang.org/t/how-to-check-the-dimension-names-before-eliminating-some-dimensions-of-a-polyhedron/111864)

<div class="topic-metadata">

**Author:** [@WuSiren](https://discourse.julialang.org/u/WuSiren)\
**Replies:** 3\
**Last updated:** [March 22, 2024, 8:12am UTC](https://discourse.julialang.org/t/how-to-check-the-dimension-names-before-eliminating-some-dimensions-of-a-polyhedron/111864 "2024-03-22T08:12:38Z")

</div>

using Polyhedra using JuMP using CDDLib poly = Model() @variable(poly, -1 \<= x \<= 1) @variable(poly, -1 \<= y \<= 1) P = polyhedron(poly, CDDLib.Library(:exact)) P\_y = eliminate(P, \[1\]) How can I make sure the index \[1\] …

---

## [\[ANN\] HampelOutliers.jl](https://discourse.julialang.org/t/ann-hampeloutliers-jl/111844)

<div class="topic-metadata">

**Author:** [@tobydriscoll](https://discourse.julialang.org/u/tobydriscoll)\
**Replies:** 4\
**Last updated:** [March 22, 2024, 5:28am UTC](https://discourse.julialang.org/t/ann-hampeloutliers-jl/111844 "2024-03-22T05:28:52Z")

</div>

I’ve written an implementation of the Hampel filter, which is a simple yet robust way to detect and replace outliers in univariate time series without smoothing legitimate jumps. It seems to offer more flexibility than a…

---

## [The julia expression of only one equation with two variables](https://discourse.julialang.org/t/the-julia-expression-of-only-one-equation-with-two-variables/111916)

<div class="topic-metadata">

**Author:** [@whucsu](https://discourse.julialang.org/u/whucsu)\
**Replies:** 2\
**Last updated:** [March 22, 2024, 4:14am UTC](https://discourse.julialang.org/t/the-julia-expression-of-only-one-equation-with-two-variables/111916 "2024-03-22T04:14:37Z")

</div>

I get one simulation with only one expression below: the P(t) and Q(t) with their derivative are the expected, all the R，C are known parameters. suppose the initial values are known ,but I don’t know how to define th…

---

## [A question about implementing double integral](https://discourse.julialang.org/t/a-question-about-implementing-double-integral/111914)

<div class="topic-metadata">

**Author:** [@swsw](https://discourse.julialang.org/u/swsw)\
**Replies:** 3\
**Last updated:** [March 22, 2024, 3:59am UTC](https://discourse.julialang.org/t/a-question-about-implementing-double-integral/111914 "2024-03-22T03:59:37Z")

</div>

I’d like to calculate the following double integral Integrate\[exp(-x)Integrate\[exp(-y),{y,0,x}\],{x,0,∞}\] I try the following in julia but i think i probably get it wrong. using QuadGK inner\_function(y) = exp(-y) …

---

## [Color conversion from terminal color to Colors.jl color](https://discourse.julialang.org/t/color-conversion-from-terminal-color-to-colors-jl-color/111889)

<div class="topic-metadata">

**Author:** [@Julia2001](https://discourse.julialang.org/u/Julia2001)\
**Replies:** 5\
**Last updated:** [March 22, 2024, 12:43am UTC](https://discourse.julialang.org/t/color-conversion-from-terminal-color-to-colors-jl-color/111889 "2024-03-22T00:43:05Z")

</div>

Hi, how can I convert the “one integer” terminal color consumed by printstyled() to a Colors.jl color usable within Luxor e.g.?

---

## [Annotations?](https://discourse.julialang.org/t/annotations/111949)

<div class="topic-metadata">

**Author:** [@mkitti](https://discourse.julialang.org/u/mkitti)\
**Replies:** 5\
**Last updated:** [March 21, 2024, 10:50pm UTC](https://discourse.julialang.org/t/annotations/111949 "2024-03-21T22:50:09Z")

</div>

What is the best way to annotate code that is both machine and human readable? My first thought would be parseable structured docstrings, but this seems too fragile. My second thought is macros that just emit the main e…

---

## [Derivative of user-defined function with variables as arguments in a PDE system](https://discourse.julialang.org/t/derivative-of-user-defined-function-with-variables-as-arguments-in-a-pde-system/111597)

<div class="topic-metadata">

**Author:** [@ehintz3](https://discourse.julialang.org/u/ehintz3)\
**Replies:** 1\
**Last updated:** [March 21, 2024, 10:47pm UTC](https://discourse.julialang.org/t/derivative-of-user-defined-function-with-variables-as-arguments-in-a-pde-system/111597 "2024-03-21T22:47:36Z")

</div>

Hello, I am new to Julia, so apologies if I am missing something basic, but I wasn’t getting any replies in New To Julia, so I’m not sure that was the right place for this post. I am trying to solve a system of PDEs usin…

---

## [Color gradient doesn't work](https://discourse.julialang.org/t/color-gradient-doesnt-work/111734)

<div class="topic-metadata">

**Author:** [@Enlil50](https://discourse.julialang.org/u/Enlil50)\
**Replies:** 6\
**Last updated:** [March 21, 2024, 9:39pm UTC](https://discourse.julialang.org/t/color-gradient-doesnt-work/111734 "2024-03-21T21:39:50Z")

</div>

I have a for loop. I want the color used to be evenly spaced from a color gradient such as inferno. I read various stuff online, but none works. for example trying the palette = : inferno in ax.plot of pythonplot gives: …

---

## [Documenter.jl: misplaced docstrings when using method signatures](https://discourse.julialang.org/t/documenter-jl-misplaced-docstrings-when-using-method-signatures/111951)

<div class="topic-metadata">

**Author:** [@matthias314](https://discourse.julialang.org/u/matthias314)\
**Replies:** 0\
**Last updated:** [March 21, 2024, 9:06pm UTC](https://discourse.julialang.org/t/documenter-jl-misplaced-docstrings-when-using-method-signatures/111951 "2024-03-21T21:06:36Z")

</div>

I get repeated docstrings in the output of Documenter.jl when I try to select methods by their signatures. Here is an example: Take this package module DocumenterTestPkg export f "vector doc" f(x::Vector{T}) where T =…

---

## [How to return a non-allocating view from a function?](https://discourse.julialang.org/t/how-to-return-a-non-allocating-view-from-a-function/111946)

<div class="topic-metadata">

**Author:** [@Ahmed\_Salih](https://discourse.julialang.org/u/Ahmed_Salih)\
**Replies:** 3\
**Last updated:** [March 21, 2024, 8:11pm UTC](https://discourse.julialang.org/t/how-to-return-a-non-allocating-view-from-a-function/111946 "2024-03-21T20:11:39Z")

</div>

Hello! I have a function in which I do a small amount of operations, before returning a view of the final result as such: # some operations .. ParticleRanges = view(ParticleSplitterLinearIndices, ParticleSplitt…

---

## [Memory allocation in recursion](https://discourse.julialang.org/t/memory-allocation-in-recursion/111929)

<div class="topic-metadata">

**Author:** [@PetrKryslUCSD](https://discourse.julialang.org/u/PetrKryslUCSD)\
**Replies:** 11\
**Last updated:** [March 21, 2024, 8:08pm UTC](https://discourse.julialang.org/t/memory-allocation-in-recursion/111929 "2024-03-21T20:08:27Z")

</div>

Does Julia allocate memory when recursing? I see a ton of memory allocation in a recursive quicksort…

---

## [Something wrong with my memory calculation](https://discourse.julialang.org/t/something-wrong-with-my-memory-calculation/111932)

<div class="topic-metadata">

**Author:** [@Clement\_Benoist](https://discourse.julialang.org/u/Clement_Benoist)\
**Replies:** 0\
**Last updated:** [March 21, 2024, 1:35pm UTC](https://discourse.julialang.org/t/something-wrong-with-my-memory-calculation/111932 "2024-03-21T13:35:41Z")

</div>

I’m coding a NeuralODE in Julia and tried to get a memory use data from the program. On python, when I change my adjoint method I get different but consistent (by method) results for my max RSS but on Julia I always get …

---

## [Warning: error closing screen](https://discourse.julialang.org/t/warning-error-closing-screen/111939)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 1\
**Last updated:** [March 21, 2024, 4:42pm UTC](https://discourse.julialang.org/t/warning-error-closing-screen/111939 "2024-03-21T16:42:33Z")

</div>

I can run my simulation fine from the REPL. You can find my code at: GitHub - aenarete/KiteControllers.jl: Controllers for airborne wind energy systems But if I run it from bash with the command: julia --project -t 7,…

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

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