# Specific Domains

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

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

**Page:** 133

---

## [Reduce memory usage for ensemble problems](https://discourse.julialang.org/t/reduce-memory-usage-for-ensemble-problems/90591)

<div class="topic-metadata">

**Author:** [@henry2004y](https://discourse.julialang.org/u/henry2004y)\
**Replies:** 4\
**Last updated:** [February 22, 2024, 4:24am UTC](https://discourse.julialang.org/t/reduce-memory-usage-for-ensemble-problems/90591 "2024-02-22T04:24:32Z")

</div>

Hi, I have a test particle tracer which is built on top of OrdinaryDiffEq.jl. This model can have 3D numerical fields as input parameters and solve ODEs for each set of initial conditions. However, when testing on the …

---

## [Hide stroke in 3D scatter using Makie.jl](https://discourse.julialang.org/t/hide-stroke-in-3d-scatter-using-makie-jl/106383)

<div class="topic-metadata">

**Author:** [@juliohm](https://discourse.julialang.org/u/juliohm)\
**Replies:** 6\
**Last updated:** [February 21, 2024, 10:31pm UTC](https://discourse.julialang.org/t/hide-stroke-in-3d-scatter-using-makie-jl/106383 "2024-02-21T22:31:42Z")

</div>

Compare the two scatter plots produced by GLMakie in 2D and 3D: scatter(rand(1000), rand(1000)) scatter(rand(1000), rand(1000), rand(1000)) Is there an option to hide the white stroke around the dots in 3D? When …

---

## [Error when using LawsonEuler method on SplitODE](https://discourse.julialang.org/t/error-when-using-lawsoneuler-method-on-splitode/110495)

<div class="topic-metadata">

**Author:** [@ian](https://discourse.julialang.org/u/ian)\
**Replies:** 1\
**Last updated:** [February 21, 2024, 9:52pm UTC](https://discourse.julialang.org/t/error-when-using-lawsoneuler-method-on-splitode/110495 "2024-02-21T21:52:32Z")

</div>

Hi all. I have a Split/Semiinear ODE I want to solve efficiently. The docs suggest the LawsonEuler algorithm. I tried to implement but am running into an error I don’t understand. Here is a MWE using DifferentialEquatio…

---

## [Plotting ODEProblem solutions with Makie](https://discourse.julialang.org/t/plotting-odeproblem-solutions-with-makie/58096)

<div class="topic-metadata">

**Author:** [@airpmb](https://discourse.julialang.org/u/airpmb)\
**Replies:** 6\
**Last updated:** [February 21, 2024, 8:39pm UTC](https://discourse.julialang.org/t/plotting-odeproblem-solutions-with-makie/58096 "2024-02-21T20:39:26Z")

</div>

Whatever (plot recipes?) allows for the solution output from solving an ODEProblem in ModelingToolkit.jl (or maybe just OrdinaryDiffEq.jl?) to be just plotted using Plots.jl is wonderful and magical. I love how it’s poss…

---

## [Lazy constraints & issues with presolve](https://discourse.julialang.org/t/lazy-constraints-issues-with-presolve/110449)

<div class="topic-metadata">

**Author:** [@salvaguerrero](https://discourse.julialang.org/u/salvaguerrero)\
**Replies:** 3\
**Last updated:** [February 21, 2024, 7:32pm UTC](https://discourse.julialang.org/t/lazy-constraints-issues-with-presolve/110449 "2024-02-21T19:32:55Z")

</div>

Good morning community, I’m trying to use lazy constraints, but I’m finding some issues. I would like to solve a problem iterative (benders decomposition) where in each iteration an additional constraint is added as la…

---

## [Finite difference Laplacian with five-point stencil](https://discourse.julialang.org/t/finite-difference-laplacian-with-five-point-stencil/25014)

<div class="topic-metadata">

**Author:** [@Vortico](https://discourse.julialang.org/u/Vortico)\
**Replies:** 8\
**Last updated:** [February 21, 2024, 7:23pm UTC](https://discourse.julialang.org/t/finite-difference-laplacian-with-five-point-stencil/25014 "2024-02-21T19:23:31Z")

</div>

This is my first day learning Julia, and I’m not sure of an elegant way to compute the finite difference Laplacian on an NxM Float32 array. \\Delta f(x,y) \\approx f(x-1,y) + f(x+1,y) + f(x,y-1) + f(x,y+1) - 4f(x,y) Base…

---

## [CairoMakie get 3d plots to respect z order](https://discourse.julialang.org/t/cairomakie-get-3d-plots-to-respect-z-order/110428)

<div class="topic-metadata">

**Author:** [@p\_f](https://discourse.julialang.org/u/p_f)\
**Replies:** 5\
**Last updated:** [February 21, 2024, 6:09pm UTC](https://discourse.julialang.org/t/cairomakie-get-3d-plots-to-respect-z-order/110428 "2024-02-21T18:09:59Z")

</div>

I’m trying to export figure as a PDF using CairoMakie. Some of my subplots are 3d axes. I’m plotting several surfaces that wrap around each other, so some parts of surface A are behind surface B and some parts are in fro…

---

## [Help to find package for DSA (Discrete Stochastic Arithmetic) implementation](https://discourse.julialang.org/t/help-to-find-package-for-dsa-discrete-stochastic-arithmetic-implementation/110391)

<div class="topic-metadata">

**Author:** [@hungpham3112](https://discourse.julialang.org/u/hungpham3112)\
**Replies:** 5\
**Last updated:** [February 21, 2024, 4:30pm UTC](https://discourse.julialang.org/t/help-to-find-package-for-dsa-discrete-stochastic-arithmetic-implementation/110391 "2024-02-21T16:30:44Z")

</div>

Hi, I’m working in a project where we use DSA method (based-on CESTAC) to validate results of numerical software. Currently, I’m using Cadna (written in C) to solve this problem. The advantage of Cadna is it has example …

---

## [Subplots CairoMakie](https://discourse.julialang.org/t/subplots-cairomakie/110504)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 4\
**Last updated:** [February 21, 2024, 9:43am UTC](https://discourse.julialang.org/t/subplots-cairomakie/110504 "2024-02-21T09:43:43Z")

</div>

Hello, i want plot N subplots and off for all subplots ticks on axis. For example i do it in python with matplotlib. How can i do this with CairoMakie? Thank you for your helps

---

## [How to hand-craft an edge into a Pluto's a computational graph?](https://discourse.julialang.org/t/how-to-hand-craft-an-edge-into-a-plutos-a-computational-graph/108619)

<div class="topic-metadata">

**Author:** [@vvbond](https://discourse.julialang.org/u/vvbond)\
**Replies:** 11\
**Last updated:** [February 21, 2024, 9:00am UTC](https://discourse.julialang.org/t/how-to-hand-craft-an-edge-into-a-plutos-a-computational-graph/108619 "2024-02-21T09:00:33Z")

</div>

Consider the following contrived example of a Pluto notebook: First of all, I have a hard time even to visualise the corresponding computational graph. Below is my attempt to indicate by the dashed arrow that the com…

---

## [Shared legend for subplots in Plotly](https://discourse.julialang.org/t/shared-legend-for-subplots-in-plotly/110497)

<div class="topic-metadata">

**Author:** [@PeX](https://discourse.julialang.org/u/PeX)\
**Replies:** 0\
**Last updated:** [February 21, 2024, 4:58am UTC](https://discourse.julialang.org/t/shared-legend-for-subplots-in-plotly/110497 "2024-02-21T04:58:12Z")

</div>

I’ve seen variations of this question, but nothing fits my case exactly so I’m trying.: I want to create three subplots in Plotly that share the same legend, I use the following code: p = \[p\_1; p\_2; p\_3\] relayout!(p, h…

---

## [Trouble writing custom rule in Enzyme: \`AssertionError: !(overwritten\[end\])\`](https://discourse.julialang.org/t/trouble-writing-custom-rule-in-enzyme-assertionerror-overwritten-end/108939)

<div class="topic-metadata">

**Author:** [@Gianluca\_Fuwa](https://discourse.julialang.org/u/Gianluca_Fuwa)\
**Replies:** 3\
**Last updated:** [February 20, 2024, 11:25pm UTC](https://discourse.julialang.org/t/trouble-writing-custom-rule-in-enzyme-assertionerror-overwritten-end/108939 "2024-02-20T23:25:57Z")

</div>

Hi, I am currently trying to use Enzyme to get the derivative of a function that takes a multidimensional array of 3x3 matrices and outputs a scalar. The caveat is that these matrices are elements of the Lie-group of sp…

---

## [LoadError: UndefVarError: \`libglpk\` not defined](https://discourse.julialang.org/t/loaderror-undefvarerror-libglpk-not-defined/110480)

<div class="topic-metadata">

**Author:** [@Carl\_Baier](https://discourse.julialang.org/u/Carl_Baier)\
**Replies:** 7\
**Last updated:** [February 20, 2024, 11:16pm UTC](https://discourse.julialang.org/t/loaderror-undefvarerror-libglpk-not-defined/110480 "2024-02-20T23:16:52Z")

</div>

Hello, i have this code: using libglpk using GLPK using JuMP using Gurobi using DataFrames using BlockDecomposition, Coluna; I\_list = \[1, 2, 3\] T\_list = \[1, 2, 3, 4, 5, 6, 7\] K\_list = \[1, 2, 3\] alpha=0.01 M = 10000 Ma…

---

## [Matrix for storage sol from DifferentialEquations](https://discourse.julialang.org/t/matrix-for-storage-sol-from-differentialequations/110464)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 2\
**Last updated:** [February 20, 2024, 5:12pm UTC](https://discourse.julialang.org/t/matrix-for-storage-sol-from-differentialequations/110464 "2024-02-20T17:12:07Z")

</div>

Hello, i solve system with different initial condition and i want to save all sol. I’m thinking about creating a matrix or vector and storing solutions there for various initial conditions. How can this be done? Thank yo…

---

## [Autodiff Enzyme reverse mode sometimes fails on on repeated evaluation](https://discourse.julialang.org/t/autodiff-enzyme-reverse-mode-sometimes-fails-on-on-repeated-evaluation/110017)

<div class="topic-metadata">

**Author:** [@dchang10](https://discourse.julialang.org/u/dchang10)\
**Replies:** 7\
**Last updated:** [February 20, 2024, 3:27pm UTC](https://discourse.julialang.org/t/autodiff-enzyme-reverse-mode-sometimes-fails-on-on-repeated-evaluation/110017 "2024-02-20T15:27:04Z")

</div>

I’m developing a package which has a function defined the following way function emission\_radius(pix::Krang.AbstractPixel, θs::T, isindir, n) where {T} α, β = screen\_coordinate(pix) θo = inclination(pix) met…

---

## [Recording of Makie hist changes output](https://discourse.julialang.org/t/recording-of-makie-hist-changes-output/110440)

<div class="topic-metadata">

**Author:** [@DrBurner](https://discourse.julialang.org/u/DrBurner)\
**Replies:** 2\
**Last updated:** [February 20, 2024, 1:11pm UTC](https://discourse.julialang.org/t/recording-of-makie-hist-changes-output/110440 "2024-02-20T13:11:26Z")

</div>

I’d like to make a short video of the build up of an histogram with a little bit of white space between the bins . Unfortunately the record block seems to ignore the width keyword. using GLMakie samples=rand(1:10,100) …

---

## [Access JuMP variables queried as string](https://discourse.julialang.org/t/access-jump-variables-queried-as-string/110393)

<div class="topic-metadata">

**Author:** [@math\_opt](https://discourse.julialang.org/u/math_opt)\
**Replies:** 10\
**Last updated:** [February 20, 2024, 12:16pm UTC](https://discourse.julialang.org/t/access-jump-variables-queried-as-string/110393 "2024-02-20T12:16:06Z")

</div>

If a user wants to know the variable values for a JuMP model that is not exposed to the user directly, what are the options? For example: using JuMP using HiGHS #some arbitrary JuMP model model = Model(HiGHS.Optimizer)…

---

## [How to extract convergence status in Optim.jl?](https://discourse.julialang.org/t/how-to-extract-convergence-status-in-optim-jl/110430)

<div class="topic-metadata">

**Author:** [@Wooheon](https://discourse.julialang.org/u/Wooheon)\
**Replies:** 2\
**Last updated:** [February 20, 2024, 12:41am UTC](https://discourse.julialang.org/t/how-to-extract-convergence-status-in-optim-jl/110430 "2024-02-20T00:41:02Z")

</div>

I ran several optimization problems with different models, but some didn’t seem to calculate the result but spit out some random results. So, I would like to know how to extract the status of the ‘optimized’ results, but…

---

## [How to pass a constraint as Lazy with Gurobi?](https://discourse.julialang.org/t/how-to-pass-a-constraint-as-lazy-with-gurobi/110294)

<div class="topic-metadata">

**Author:** [@K5Julien](https://discourse.julialang.org/u/K5Julien)\
**Replies:** 5\
**Last updated:** [February 19, 2024, 5:57pm UTC](https://discourse.julialang.org/t/how-to-pass-a-constraint-as-lazy-with-gurobi/110294 "2024-02-19T17:57:50Z")

</div>

I want to pass a constraint as Lazy in the beginning of my JuMP model. Note that I do not want it to be a lazy constraint called in a Callback. How am I supposed to do it? The corresponding Gurobi parameter is: It …

---

## [Makie plotting: GLib-GIO-WARNING](https://discourse.julialang.org/t/makie-plotting-glib-gio-warning/64478)

<div class="topic-metadata">

**Author:** [@JosephPollacco](https://discourse.julialang.org/u/JosephPollacco)\
**Replies:** 2\
**Last updated:** [February 19, 2024, 4:40pm UTC](https://discourse.julialang.org/t/makie-plotting-glib-gio-warning/64478 "2024-02-19T16:40:13Z")

</div>

Hallo :grinning\_face\_with\_smiling\_eyes:, I am using the excellent Makie plotting package. When I am compiling the code I am getting the following Warnings, GLib-GIO-WARNING \*which I suspect is caused by the package Cair…

---

## [Methods to solve system of linear equations with equiality constraint](https://discourse.julialang.org/t/methods-to-solve-system-of-linear-equations-with-equiality-constraint/110403)

<div class="topic-metadata">

**Author:** [@Gregstrq](https://discourse.julialang.org/u/Gregstrq)\
**Replies:** 4\
**Last updated:** [February 19, 2024, 3:00pm UTC](https://discourse.julialang.org/t/methods-to-solve-system-of-linear-equations-with-equiality-constraint/110403 "2024-02-19T15:00:02Z")

</div>

Hello, everyone. I want to solve a system of linear equations Ax = b with some linear equality constraints. A is square but non-symmetric. The equation corresponds to a linearized stationary Boltzman equation, and the…

---

## [Client code for SSE (Server Sent Events) with HTTP.jl](https://discourse.julialang.org/t/client-code-for-sse-server-sent-events-with-http-jl/110386)

<div class="topic-metadata">

**Author:** [@screw\_dog](https://discourse.julialang.org/u/screw_dog)\
**Replies:** 1\
**Last updated:** [February 19, 2024, 2:08pm UTC](https://discourse.julialang.org/t/client-code-for-sse-server-sent-events-with-http-jl/110386 "2024-02-19T14:08:07Z")

</div>

I’m attempting to receive SSE from a Mastodon instance, following the API described in the official docs. Firstly, I can’t get HTTP.open to work. Attempting: HTTP.open("get", "https://aus.social/api/v1/streaming/public…

---

## [Can't run a Resnet code](https://discourse.julialang.org/t/cant-run-a-resnet-code/110303)

<div class="topic-metadata">

**Author:** [@delcencen](https://discourse.julialang.org/u/delcencen)\
**Replies:** 3\
**Last updated:** [February 19, 2024, 1:12pm UTC](https://discourse.julialang.org/t/cant-run-a-resnet-code/110303 "2024-02-19T13:12:10Z")

</div>

Hello, in this code i ave the 2 following errors : -ERROR: LoadError: UndefVarError: \`Resnet\` not defined -ERROR: LoadError: UndefVarError: \`get\_classes\` not defined the code of my book : using Flux using M…

---

## [Rolling (time) average of a timeseries of images](https://discourse.julialang.org/t/rolling-time-average-of-a-timeseries-of-images/110352)

<div class="topic-metadata">

**Author:** [@JonathanH](https://discourse.julialang.org/u/JonathanH)\
**Replies:** 3\
**Last updated:** [February 19, 2024, 9:59am UTC](https://discourse.julialang.org/t/rolling-time-average-of-a-timeseries-of-images/110352 "2024-02-19T09:59:36Z")

</div>

Hi all, I am new to Julia, and this may be a very straightforward question. I have a timeseries of images, and I am trying to calculate a pixel-wise rolling time average. The obvious choice for doing this would be the m…

---

## [Do not update neural network weights with a value of 0](https://discourse.julialang.org/t/do-not-update-neural-network-weights-with-a-value-of-0/110335)

<div class="topic-metadata">

**Author:** [@kuanchiun](https://discourse.julialang.org/u/kuanchiun)\
**Replies:** 5\
**Last updated:** [February 19, 2024, 2:42am UTC](https://discourse.julialang.org/t/do-not-update-neural-network-weights-with-a-value-of-0/110335 "2024-02-19T02:42:38Z")

</div>

Hello, I have just been exposed to the field of deep learning for a short time. Recently I was trying to create the 3-dimensional Neural Network, by giving the 3D Array as a pretrained weights. Each 2D matrix in the wei…

---

## [Read\_parquet() not defined?](https://discourse.julialang.org/t/read-parquet-not-defined/110380)

<div class="topic-metadata">

**Author:** [@austin-putz](https://discourse.julialang.org/u/austin-putz)\
**Replies:** 9\
**Last updated:** [February 19, 2024, 1:02am UTC](https://discourse.julialang.org/t/read-parquet-not-defined/110380 "2024-02-19T01:02:26Z")

</div>

Hello, I feel kind of silly here. No idea what I did. I couldn’t seem to load a dataset with Arrow from R (feather). So I jumped to Parquet to test. I used the arrow package in R, have a .parquet file on my hard disk and…

---

## [Labeling a categorical color dimension in 3d scatter in Makie](https://discourse.julialang.org/t/labeling-a-categorical-color-dimension-in-3d-scatter-in-makie/110072)

<div class="topic-metadata">

**Author:** [@tp2750](https://discourse.julialang.org/u/tp2750)\
**Replies:** 1\
**Last updated:** [February 18, 2024, 6:43pm UTC](https://discourse.julialang.org/t/labeling-a-categorical-color-dimension-in-3d-scatter-in-makie/110072 "2024-02-18T18:43:14Z")

</div>

I’m trying to do a 3D Iris plot in Makie, but I can not figure out how to label the color dimension. Here is what I do: using DataFrames using WGLMakie using RDatasets using CategoricalArrays iris = dataset("datasets",…

---

## [Solving linear system with sum of kronecker products](https://discourse.julialang.org/t/solving-linear-system-with-sum-of-kronecker-products/109253)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 21\
**Last updated:** [February 18, 2024, 11:35am UTC](https://discourse.julialang.org/t/solving-linear-system-with-sum-of-kronecker-products/109253 "2024-02-18T11:35:39Z")

</div>

Suppose I have ((A \\otimes B) + (C \\otimes D))x = v where A and C and B and D have the same size, respectively. Suppose also that I particularly want to avoid instantiating the Kronecker products since they are large. …

---

## [DifferentialEquations.jl -- change in \`remake\`?](https://discourse.julialang.org/t/differentialequations-jl-change-in-remake/110283)

<div class="topic-metadata">

**Author:** [@BLI](https://discourse.julialang.org/u/BLI)\
**Replies:** 8\
**Last updated:** [February 18, 2024, 10:48am UTC](https://discourse.julialang.org/t/differentialequations-jl-change-in-remake/110283 "2024-02-18T10:48:07Z")

</div>

Has there been a change in the remake function recently? The following code worked one month ago: prob\_rm = remake(prob, u0=\[\], p=p\_modify) sol = solve(prob\_rm, QBDF(), reltol=1e-5,tstops=t\_u\_switch) Now, I get an err…

---

## [Optimizing a CUDA-based Loss Function for Neural Network Training](https://discourse.julialang.org/t/optimizing-a-cuda-based-loss-function-for-neural-network-training/110365)

<div class="topic-metadata">

**Author:** [@josemanuel22](https://discourse.julialang.org/u/josemanuel22)\
**Replies:** 0\
**Last updated:** [February 18, 2024, 8:54am UTC](https://discourse.julialang.org/t/optimizing-a-cuda-based-loss-function-for-neural-network-training/110365 "2024-02-18T08:54:36Z")

</div>

This post is hidden. You deleted this post 8 hours ago. I have been working on the following loss function to train neural networks based on a statistical invariant. I have tried to port the function to GPU/CUDA, and I …

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

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