# Latest

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

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

**Page:** 159

---

## [Matrix multiply breakthrough, AlphaTensor (could also do for other algorithms): "AlphaTensor discovered algorithms that are more efficient than the state of the art for many matrix sizes."](https://discourse.julialang.org/t/matrix-multiply-breakthrough-alphatensor-could-also-do-for-other-algorithms-alphatensor-discovered-algorithms-that-are-more-efficient-than-the-state-of-the-art-for-many-matrix-sizes/88455)

<div class="topic-metadata">

**Author:** [@Palli](https://discourse.julialang.org/u/Palli)\
**Replies:** 6\
**Last updated:** [May 30, 2025, 9:52am UTC](https://discourse.julialang.org/t/matrix-multiply-breakthrough-alphatensor-could-also-do-for-other-algorithms-alphatensor-discovered-algorithms-that-are-more-efficient-than-the-state-of-the-art-for-many-matrix-sizes/88455 "2025-05-30T09:52:09Z")

</div>

@Elrod I believe you and others (OpenBLAS, MKL etc) use the Strassen algorithm, which has now been practically improved on. TL;DR, multiplying 10-20% faster is considered a breakthrough, now achieved: we adapted Alpha…

---

## [CUDA suddenly crashes with check-bounds=no, used to work fine](https://discourse.julialang.org/t/cuda-suddenly-crashes-with-check-bounds-no-used-to-work-fine/129245)

<div class="topic-metadata">

**Author:** [@0samuraiE](https://discourse.julialang.org/u/0samuraiE)\
**Replies:** 1\
**Last updated:** [May 30, 2025, 9:43am UTC](https://discourse.julialang.org/t/cuda-suddenly-crashes-with-check-bounds-no-used-to-work-fine/129245 "2025-05-30T09:43:52Z")

</div>

Hi all, I have some code that used to work perfectly fine with check-bounds=no, but recently it started crashing with a bounds error. If I run it with check-bounds=yes, it works without any issues. I noticed that the s…

---

## [Makie datashader color corruption](https://discourse.julialang.org/t/makie-datashader-color-corruption/129460)

<div class="topic-metadata">

**Author:** [@cmarcotte](https://discourse.julialang.org/u/cmarcotte)\
**Replies:** 6\
**Last updated:** [May 30, 2025, 9:39am UTC](https://discourse.julialang.org/t/makie-datashader-color-corruption/129460 "2025-05-30T09:39:35Z")

</div>

Hello, I am running into a bizarre coloration issue with Makie (GLMakie) with the datashader – this leads to a bright blue color on-screen, while the saved figure is red (as it is in the example). Additionally, the on-…

---

## [Turing with both Poisson and Gaussian noise](https://discourse.julialang.org/t/turing-with-both-poisson-and-gaussian-noise/129369)

<div class="topic-metadata">

**Author:** [@Cyan](https://discourse.julialang.org/u/Cyan)\
**Replies:** 3\
**Last updated:** [May 30, 2025, 7:54am UTC](https://discourse.julialang.org/t/turing-with-both-poisson-and-gaussian-noise/129369 "2025-05-30T07:54:18Z")

</div>

I have a set of data with both Poisson and Gaussian noise, which distribution I should choose in likelihood? My dataset looks like obs = data + Poisson + Gaussian PS My data is a series of images, it contains Gaussia…

---

## [What should be the acceptable tolerance value between the lower bound and the estimated upper bound for the results of SDDP.jl to be reliable?](https://discourse.julialang.org/t/what-should-be-the-acceptable-tolerance-value-between-the-lower-bound-and-the-estimated-upper-bound-for-the-results-of-sddp-jl-to-be-reliable/129450)

<div class="topic-metadata">

**Author:** [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)\
**Replies:** 4\
**Last updated:** [May 30, 2025, 7:40am UTC](https://discourse.julialang.org/t/what-should-be-the-acceptable-tolerance-value-between-the-lower-bound-and-the-estimated-upper-bound-for-the-results-of-sddp-jl-to-be-reliable/129450 "2025-05-30T07:40:04Z")

</div>

I am dealing with a multi-stage mixed-integer stochastic global supply chain optimization problem. In this regard, I want to determine a tolerance value for the gap between the lower bound and the estimated upper bound …

---

## [ExtendableFEM - Penalty Method - Total Pressure Inlet Condition](https://discourse.julialang.org/t/extendablefem-penalty-method-total-pressure-inlet-condition/129448)

<div class="topic-metadata">

**Author:** [@ziolai](https://discourse.julialang.org/u/ziolai)\
**Replies:** 2\
**Last updated:** [May 30, 2025, 7:30am UTC](https://discourse.julialang.org/t/extendablefem-penalty-method-total-pressure-inlet-condition/129448 "2025-05-30T07:30:38Z")

</div>

Greetings. Does it make sense to attempt extending ExtendableFEM.jl/examples/Example235\_StokesIteratedPenalty.jl at master · WIAS-PDELib/ExtendableFEM.jl · GitHub to implement a total pressure inlet boundary condition? …

---

## [How can the Julia language write a structure to a file and read that structure from the file?](https://discourse.julialang.org/t/how-can-the-julia-language-write-a-structure-to-a-file-and-read-that-structure-from-the-file/129436)

<div class="topic-metadata">

**Author:** [@dk-zen](https://discourse.julialang.org/u/dk-zen)\
**Replies:** 5\
**Last updated:** [May 30, 2025, 4:57am UTC](https://discourse.julialang.org/t/how-can-the-julia-language-write-a-structure-to-a-file-and-read-that-structure-from-the-file/129436 "2025-05-30T04:57:50Z")

</div>

How can the Julia language write a structure to a file and read that structure from the file? I tried it myself. Although I could read it correctly, when there are dozens of keywords in the structure, the reading process…

---

## [Where is the "rate\_a" in PowerModels.jl from](https://discourse.julialang.org/t/where-is-the-rate-a-in-powermodels-jl-from/129424)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 8\
**Last updated:** [May 29, 2025, 11:52pm UTC](https://discourse.julialang.org/t/where-is-the-rate-a-in-powermodels-jl-from/129424 "2025-05-29T23:52:24Z")

</div>

The data provided by PowerModels.jl is inconsistent with my local raw Matpower file. The latter is all zero (seems to mean that there is no limit), but the former has a vector of numbers that I can’t tell where they com…

---

## [Get an eigvec by eigval, using eigvecs, only allow SymTridiagonal?](https://discourse.julialang.org/t/get-an-eigvec-by-eigval-using-eigvecs-only-allow-symtridiagonal/127408)

<div class="topic-metadata">

**Author:** [@WalterMadelim](https://discourse.julialang.org/u/WalterMadelim)\
**Replies:** 12\
**Last updated:** [May 29, 2025, 10:40pm UTC](https://discourse.julialang.org/t/get-an-eigvec-by-eigval-using-eigvecs-only-allow-symtridiagonal/127408 "2025-05-29T22:40:55Z")

</div>

Hi, there I wonder why the eigvecs method does not work for a A::Symmetric{Int64, Matrix{Int64}}. using LinearAlgebra A = SymTridiagonal(\[1.; 2.; 1.\], \[2.; 3.\]) val3 = eigvals(A)\[3\] vec3 = reshape(eigvecs(A, \[val3\]), (…

---

## [How to point Julia to a python executable](https://discourse.julialang.org/t/how-to-point-julia-to-a-python-executable/129388)

<div class="topic-metadata">

**Author:** [@keynescoefen](https://discourse.julialang.org/u/keynescoefen)\
**Replies:** 9\
**Last updated:** [May 29, 2025, 10:30pm UTC](https://discourse.julialang.org/t/how-to-point-julia-to-a-python-executable/129388 "2025-05-29T22:30:34Z")

</div>

Hello, When trying to precompile PythonCall I received the following error IOError: could not spawn \`/Users/jerome/.julia/environments/v1.11/.CondaPkg/.pixi/envs/default/bin/python --version\`: no such file or directory …

---

## [DifferentialEquations.jl Mass Matrix DAE Explicit](https://discourse.julialang.org/t/differentialequations-jl-mass-matrix-dae-explicit/129447)

<div class="topic-metadata">

**Author:** [@ziolai](https://discourse.julialang.org/u/ziolai)\
**Replies:** 3\
**Last updated:** [May 29, 2025, 8:07pm UTC](https://discourse.julialang.org/t/differentialequations-jl-mass-matrix-dae-explicit/129447 "2025-05-29T20:07:25Z")

</div>

Are any of the methods listed at Mass Matrix and Fully Implicit DAE Solvers · DifferentialEquations.jl explicit in time, and therefore not require a Jacobian? Reason for asking is that solving time-dependent Stokes flow…

---

## [Type instability with Type{T} selectors in keyword arguments](https://discourse.julialang.org/t/type-instability-with-type-t-selectors-in-keyword-arguments/129439)

<div class="topic-metadata">

**Author:** [@marcsgil](https://discourse.julialang.org/u/marcsgil)\
**Replies:** 13\
**Last updated:** [May 29, 2025, 7:47pm UTC](https://discourse.julialang.org/t/type-instability-with-type-t-selectors-in-keyword-arguments/129439 "2025-05-29T19:47:02Z")

</div>

I’m trying to write a function that outputs arrays that are filled according to a certain computation. I want the user to be able to select the type of the array according to an argument. Consider these two example func…

---

## [How to solve infinite horizon problem under Makovian policy graph in SDDP.jl?](https://discourse.julialang.org/t/how-to-solve-infinite-horizon-problem-under-makovian-policy-graph-in-sddp-jl/129420)

<div class="topic-metadata">

**Author:** [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)\
**Replies:** 16\
**Last updated:** [May 29, 2025, 7:28pm UTC](https://discourse.julialang.org/t/how-to-solve-infinite-horizon-problem-under-makovian-policy-graph-in-sddp-jl/129420 "2025-05-29T19:28:08Z")

</div>

I am dealing with a multi-stage mixed-integer stochastic global supply chain optimization problem and want to solve infinite-horizon problem instances under a Markovian policy graph. I am using the Markovian policy graph…

---

## [Passing a Cache to optimization function](https://discourse.julialang.org/t/passing-a-cache-to-optimization-function/129321)

<div class="topic-metadata">

**Author:** [@miguelborrero](https://discourse.julialang.org/u/miguelborrero)\
**Replies:** 21\
**Last updated:** [May 29, 2025, 6:24pm UTC](https://discourse.julialang.org/t/passing-a-cache-to-optimization-function/129321 "2025-05-29T18:24:57Z")

</div>

Hi! In my original problem the computation of the objective value function is expensive and during the computation I can collaterally compute the constraints of my problem. To make my problem more efficient I tried to i…

---

## [Incorrect manifest after adding package as dependency that has forked package as dependency](https://discourse.julialang.org/t/incorrect-manifest-after-adding-package-as-dependency-that-has-forked-package-as-dependency/129435)

<div class="topic-metadata">

**Author:** [@KSepetanc](https://discourse.julialang.org/u/KSepetanc)\
**Replies:** 6\
**Last updated:** [May 29, 2025, 5:41pm UTC](https://discourse.julialang.org/t/incorrect-manifest-after-adding-package-as-dependency-that-has-forked-package-as-dependency/129435 "2025-05-29T17:41:43Z")

</div>

What I am doing is creating a custom package. In that package I add forked package as dependency and resulting manifest is correct (so my package has both manifest.toml and project.toml). However, once I add my package t…

---

## [Multiply vector elementwise by \[1, v, v, ...\]](https://discourse.julialang.org/t/multiply-vector-elementwise-by-1-v-v/129442)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 2\
**Last updated:** [May 29, 2025, 3:28pm UTC](https://discourse.julialang.org/t/multiply-vector-elementwise-by-1-v-v/129442 "2025-05-29T15:28:22Z")

</div>

I have a calculation where I need to multiple various vectors elementwise by \[1, v, v, ...\]. The inputs may be any kind of vector and I want code to adapt accordingly, ie the input x is an SVector, the output should als…

---

## [Universal Differential Equations](https://discourse.julialang.org/t/universal-differential-equations/121363)

<div class="topic-metadata">

**Author:** [@Piyush](https://discourse.julialang.org/u/Piyush)\
**Replies:** 4\
**Last updated:** [May 29, 2025, 3:13pm UTC](https://discourse.julialang.org/t/universal-differential-equations/121363 "2025-05-29T15:13:43Z")

</div>

Hi, I have been working on developing hybrid models specifically Universal Differential Equations (UDE) using Julia recently. I wanted to know if I can make use of \`EnsembleProblem’ along with hybrid model (mechanistic …

---

## [Metaheuristics.jl ECA() requires redefinition between solves, making it incompatible with Optimization.jl Ensemble Problems](https://discourse.julialang.org/t/metaheuristics-jl-eca-requires-redefinition-between-solves-making-it-incompatible-with-optimization-jl-ensemble-problems/129286)

<div class="topic-metadata">

**Author:** [@johnb](https://discourse.julialang.org/u/johnb)\
**Replies:** 1\
**Last updated:** [May 29, 2025, 2:55pm UTC](https://discourse.julialang.org/t/metaheuristics-jl-eca-requires-redefinition-between-solves-making-it-incompatible-with-optimization-jl-ensemble-problems/129286 "2025-05-29T14:55:54Z")

</div>

I’m performing an optimization where the number of parameters used is varied across the ensemble. (e.g. prob 1 uses 1 parameter, prob 2 uses 2 parameters) Unfortunately, the algorithm that seems to work best with my spec…

---

## [BMS Current Control using Callbacks in ModelingToolKit](https://discourse.julialang.org/t/bms-current-control-using-callbacks-in-modelingtoolkit/129336)

<div class="topic-metadata">

**Author:** [@madtorresga](https://discourse.julialang.org/u/madtorresga)\
**Replies:** 2\
**Last updated:** [May 29, 2025, 1:38pm UTC](https://discourse.julialang.org/t/bms-current-control-using-callbacks-in-modelingtoolkit/129336 "2025-05-29T13:38:17Z")

</div>

I want to implement a battery management system that simulates the internal behaviour of a cell (ECM), and controls the incoming current based on a voltage limit (here 3.8V). I’m trying to do this using ModelingToolkit.j…

---

## [Genie Builder 1.0 Now in VS Code: Build Genie Web Apps around your Julia code 10x faster](https://discourse.julialang.org/t/genie-builder-1-0-now-in-vs-code-build-genie-web-apps-around-your-julia-code-10x-faster/114481)

<div class="topic-metadata">

**Author:** [@essenciary](https://discourse.julialang.org/u/essenciary)\
**Replies:** 3\
**Last updated:** [May 29, 2025, 11:07am UTC](https://discourse.julialang.org/t/genie-builder-1-0-now-in-vs-code-build-genie-web-apps-around-your-julia-code-10x-faster/114481 "2025-05-29T11:07:49Z")

</div>

Hi everybody, Super happy to share that Genie Builder 1.0 is now available on the VS Code Marketplace! Genie Builder is our new and powerful tool that makes it easy and quick to build web applications around your Julia …

---

## [Facing high loss values when training the SciML model](https://discourse.julialang.org/t/facing-high-loss-values-when-training-the-sciml-model/129273)

<div class="topic-metadata">

**Author:** [@Rohan\_Kumar\_123](https://discourse.julialang.org/u/Rohan_Kumar_123)\
**Replies:** 1\
**Last updated:** [May 29, 2025, 10:59am UTC](https://discourse.julialang.org/t/facing-high-loss-values-when-training-the-sciml-model/129273 "2025-05-29T10:59:47Z")

</div>

Greetings, i am currently working on a problem to predict the state vectors of a satellite using a dataset which is generated by a propagator. I have made a Universal Differential Equation(UDE) to predict the state vecto…

---

## [Test a Julia Package as a Git Hook](https://discourse.julialang.org/t/test-a-julia-package-as-a-git-hook/39140)

<div class="topic-metadata">

**Author:** [@logankilpatrick](https://discourse.julialang.org/u/logankilpatrick)\
**Replies:** 6\
**Last updated:** [May 29, 2025, 9:30am UTC](https://discourse.julialang.org/t/test-a-julia-package-as-a-git-hook/39140 "2025-05-29T09:30:12Z")

</div>

Hey all, I want to try to set my git configuration such that when I commit and before I push, my git hooks fire and test my certain package. Is anyone doing this with Julia and can share how it’s done?

---

## [Read from file a (Convex.jl) problem](https://discourse.julialang.org/t/read-from-file-a-convex-jl-problem/129402)

<div class="topic-metadata">

**Author:** [@CaG21](https://discourse.julialang.org/u/CaG21)\
**Replies:** 4\
**Last updated:** [May 29, 2025, 6:20am UTC](https://discourse.julialang.org/t/read-from-file-a-convex-jl-problem/129402 "2025-05-29T06:20:38Z")

</div>

I saved a convex problem, not solved yet, created by using Convex.jl, in a MOF/JSON file, with the write\_to\_file function. Now, I want to read that problem and start the opimization using the function Convex.solve!. How…

---

## [Launching external programs from Julia](https://discourse.julialang.org/t/launching-external-programs-from-julia/129419)

<div class="topic-metadata">

**Author:** [@mpeters2](https://discourse.julialang.org/u/mpeters2)\
**Replies:** 13\
**Last updated:** [May 29, 2025, 2:19am UTC](https://discourse.julialang.org/t/launching-external-programs-from-julia/129419 "2025-05-29T02:19:54Z")

</div>

I’ve been trying to figure out how to launch an external program from Julia, and I am stumped. In Python, I can do this: os.system("exec python -m PyQt5.uic.pyuic evaluate.ui -o ui\_Evaluate.py -x") I’ve tried most of …

---

## [SciMLOperators "\*" function for TransposedOperator](https://discourse.julialang.org/t/scimloperators-function-for-transposedoperator/129422)

<div class="topic-metadata">

**Author:** [@fdekerme](https://discourse.julialang.org/u/fdekerme)\
**Replies:** 2\
**Last updated:** [May 29, 2025, 1:55am UTC](https://discourse.julialang.org/t/scimloperators-function-for-transposedoperator/129422 "2025-05-29T01:55:57Z")

</div>

Hello :grinning:, I want to use SciMLOperators.jl with LinearSolve.jl to solve a linear system like \\left( A^T V^{-1} A + \\lambda (W - I)^T (W - I) \\right) x = A^T V^{-1} \\mu I created a FunctionOperator for A, V\_in…

---

## [DBInterface prepared statement with named parameters not working](https://discourse.julialang.org/t/dbinterface-prepared-statement-with-named-parameters-not-working/129287)

<div class="topic-metadata">

**Author:** [@slwu89](https://discourse.julialang.org/u/slwu89)\
**Replies:** 6\
**Last updated:** [May 29, 2025, 1:19am UTC](https://discourse.julialang.org/t/dbinterface-prepared-statement-with-named-parameters-not-working/129287 "2025-05-29T01:19:52Z")

</div>

Based on reading the docs at GitHub - JuliaDatabases/DBInterface.jl: Database interface definitions for Julia, I thought that I could use named parameters in a prepared statement to update data, which did not work. Howev…

---

## [Remake problem with modeltoolkitized ODE function](https://discourse.julialang.org/t/remake-problem-with-modeltoolkitized-ode-function/118339)

<div class="topic-metadata">

**Author:** [@Vitor\_Patricio\_Canta](https://discourse.julialang.org/u/Vitor_Patricio_Canta)\
**Replies:** 1\
**Last updated:** [May 29, 2025, 12:08am UTC](https://discourse.julialang.org/t/remake-problem-with-modeltoolkitized-ode-function/118339 "2025-05-29T00:08:40Z")

</div>

Hi, I have a model that follows more or less the reccomendations of the “Solving Large Stiff Equations” tutorial that I was trying to optimize the parameters for. However I have encountered a problem with the remake app…

---

## [How to resolve a nonlinear equations system in GPU (parallel)](https://discourse.julialang.org/t/how-to-resolve-a-nonlinear-equations-system-in-gpu-parallel/118226)

<div class="topic-metadata">

**Author:** [@Juan\_jose](https://discourse.julialang.org/u/Juan_jose)\
**Replies:** 4\
**Last updated:** [May 29, 2025, 12:05am UTC](https://discourse.julialang.org/t/how-to-resolve-a-nonlinear-equations-system-in-gpu-parallel/118226 "2025-05-29T00:05:20Z")

</div>

Hello, I’m new using Julia and I’m trying to solve a large system of nonlinear equations (\>768 equations, \>200 unknowns). I was reading the documentation for NonlinearSolve.jl, and although it mentioned that for solvers …

---

## [Solving ODE on GPU from Python with DifferentialEquations.jl](https://discourse.julialang.org/t/solving-ode-on-gpu-from-python-with-differentialequations-jl/109999)

<div class="topic-metadata">

**Author:** [@SantiagoMorandi](https://discourse.julialang.org/u/SantiagoMorandi)\
**Replies:** 11\
**Last updated:** [May 29, 2025, 12:02am UTC](https://discourse.julialang.org/t/solving-ode-on-gpu-from-python-with-differentialequations-jl/109999 "2025-05-29T00:02:13Z")

</div>

Hi everyone, I am new in Julia. My goal is to solve stiff ODE systems (N\>1000) representing chemical reaction networks up to steady-state conditions. Since the Python Scipy implementation is too slow for the purpose, I…

---

## [Calling Python function JIT compiled with JAX from Julia without overhead](https://discourse.julialang.org/t/calling-python-function-jit-compiled-with-jax-from-julia-without-overhead/123552)

<div class="topic-metadata">

**Author:** [@sefffal](https://discourse.julialang.org/u/sefffal)\
**Replies:** 25\
**Last updated:** [May 28, 2025, 11:40pm UTC](https://discourse.julialang.org/t/calling-python-function-jit-compiled-with-jax-from-julia-without-overhead/123552 "2025-05-28T23:40:44Z")

</div>

It could be very useful to call JIT compiled JAX functions from Python libraries from Julia, and in turn, to be able for Julia libraries to call user-provided JAX functions. This is possible today via PythonCall, e.g.: …

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

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