# Specific Domains

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

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

**Page:** 144

---

## [Differential Equations running out of GPU memory](https://discourse.julialang.org/t/differential-equations-running-out-of-gpu-memory/101048)

<div class="topic-metadata">

**Author:** [@Tommy\_Fischer](https://discourse.julialang.org/u/Tommy_Fischer)\
**Replies:** 16\
**Last updated:** [January 7, 2024, 2:53pm UTC](https://discourse.julialang.org/t/differential-equations-running-out-of-gpu-memory/101048 "2024-01-07T14:53:11Z")

</div>

Hi all, I posted earlier about minimising memory usage on gpu modelling a system here. Since then I’ve tidied up my code a bit and have managed to reduce memory usage, but it still runs out of memory when I try to scan o…

---

## [In DataDrivenDiffEq, how to get SINDy results as an array of numeric coefficients?](https://discourse.julialang.org/t/in-datadrivendiffeq-how-to-get-sindy-results-as-an-array-of-numeric-coefficients/108449)

<div class="topic-metadata">

**Author:** [@kmd](https://discourse.julialang.org/u/kmd)\
**Replies:** 1\
**Last updated:** [January 7, 2024, 1:56pm UTC](https://discourse.julialang.org/t/in-datadrivendiffeq-how-to-get-sindy-results-as-an-array-of-numeric-coefficients/108449 "2024-01-07T13:56:08Z")

</div>

Hi, I would like to obtain the result of a SINDy model fitting as an array of coefficients I can convert to Float64. I can do this in PySINDy with the function \`\`model.coefficients()‘’, which will return the weights in…

---

## [How to get type of a column using Tables.Schema from Tables.jl?](https://discourse.julialang.org/t/how-to-get-type-of-a-column-using-tables-schema-from-tables-jl/108432)

<div class="topic-metadata">

**Author:** [@FemtoTrader](https://discourse.julialang.org/u/FemtoTrader)\
**Replies:** 7\
**Last updated:** [January 6, 2024, 11:22pm UTC](https://discourse.julialang.org/t/how-to-get-type-of-a-column-using-tables-schema-from-tables-jl/108432 "2024-01-06T23:22:24Z")

</div>

Hello, I have a TSFrame like so julia\> table 50×5 TSFrame with Date Index Index Open High Low Close Volume Date Float64 Float64 Float64 Float64 Float64 ───────────────────────────────…

---

## [Is there round(..., digits=2) in CUDA?](https://discourse.julialang.org/t/is-there-round-digits-2-in-cuda/108450)

<div class="topic-metadata">

**Author:** [@BMval](https://discourse.julialang.org/u/BMval)\
**Replies:** 2\
**Last updated:** [January 7, 2024, 7:51am UTC](https://discourse.julialang.org/t/is-there-round-digits-2-in-cuda/108450 "2024-01-07T07:51:14Z")

</div>

I would like to use the function round(a, digits=2) on the GPU, but it seems that this is not possible. Is there another function for this purpose? Thank you in advance

---

## [Value.(x) with direct mode](https://discourse.julialang.org/t/value-x-with-direct-mode/108428)

<div class="topic-metadata">

**Author:** [@SingingKim](https://discourse.julialang.org/u/SingingKim)\
**Replies:** 2\
**Last updated:** [January 7, 2024, 4:50am UTC](https://discourse.julialang.org/t/value-x-with-direct-mode/108428 "2024-01-07T04:50:57Z")

</div>

I’m trying to reconstruct the optimization problem from the LP file that was generated from some commercial tool that has limited flexibility. I found a few relevant postings here and learned that I have to use the dire…

---

## [ReversedDiff/Zygote with SDE DifferentialEquations fails to compute gradients after a certain number of parameters](https://discourse.julialang.org/t/reverseddiff-zygote-with-sde-differentialequations-fails-to-compute-gradients-after-a-certain-number-of-parameters/108451)

<div class="topic-metadata">

**Author:** [@MatthieuDarcy](https://discourse.julialang.org/u/MatthieuDarcy)\
**Replies:** 2\
**Last updated:** [January 7, 2024, 3:37am UTC](https://discourse.julialang.org/t/reverseddiff-zygote-with-sde-differentialequations-fails-to-compute-gradients-after-a-certain-number-of-parameters/108451 "2024-01-07T03:37:10Z")

</div>

Hello, I’m trying to solve a parameter estimation problem for an SDE. The parameters are a small, one-layer neural network (technically, it’s a random feature model). My parameters are contained in one matrix of size (4…

---

## [Gurobi: best bound 0.0e+00, gap 100.0000%](https://discourse.julialang.org/t/gurobi-best-bound-0-0e-00-gap-100-0000/108191)

<div class="topic-metadata">

**Author:** [@quais](https://discourse.julialang.org/u/quais)\
**Replies:** 7\
**Last updated:** [January 6, 2024, 9:13pm UTC](https://discourse.julialang.org/t/gurobi-best-bound-0-0e-00-gap-100-0000/108191 "2024-01-06T21:13:28Z")

</div>

Hi, I am currently working on a cost minimization model. However, I have the problem when solving with Gurobi (v10) that my best bounds is stuck at 0.0 all the time. As a result, my gap remains constant at 100%. Addition…

---

## [How to improve performance in a function that repeatedly defines and multiplies matrices](https://discourse.julialang.org/t/how-to-improve-performance-in-a-function-that-repeatedly-defines-and-multiplies-matrices/108153)

<div class="topic-metadata">

**Author:** [@Uranium238](https://discourse.julialang.org/u/Uranium238)\
**Replies:** 64\
**Last updated:** [January 6, 2024, 2:22pm UTC](https://discourse.julialang.org/t/how-to-improve-performance-in-a-function-that-repeatedly-defines-and-multiplies-matrices/108153 "2024-01-06T14:22:15Z")

</div>

I am trying to write a code that involves running 3 nested loops, and calling a computationally expensive function within the loop. The expense mainly comes from defining , multiplying and inverting large \[~100x100\] matr…

---

## [Checkpointing and adjoints for custom kernels](https://discourse.julialang.org/t/checkpointing-and-adjoints-for-custom-kernels/108281)

<div class="topic-metadata">

**Author:** [@smartalecH](https://discourse.julialang.org/u/smartalecH)\
**Replies:** 2\
**Last updated:** [January 6, 2024, 12:32am UTC](https://discourse.julialang.org/t/checkpointing-and-adjoints-for-custom-kernels/108281 "2024-01-06T00:32:52Z")

</div>

I’ve written a PDE finite-difference solver that uses KernelAbstractions.jl to update my stencil each time step. In my previous life, when I’d write this stuff in C/C++, I’d roll my own adjoints, and then hook them up t…

---

## [Plotting points and boxplots on the same graph in Gadfly](https://discourse.julialang.org/t/plotting-points-and-boxplots-on-the-same-graph-in-gadfly/108413)

<div class="topic-metadata">

**Author:** [@bottaw](https://discourse.julialang.org/u/bottaw)\
**Replies:** 0\
**Last updated:** [January 5, 2024, 10:48pm UTC](https://discourse.julialang.org/t/plotting-points-and-boxplots-on-the-same-graph-in-gadfly/108413 "2024-01-05T22:48:23Z")

</div>

I am trying to make a gadfly graph that includes a series of data points on it with a minimalist boxplot of those points right next to it for statistical comparison. Something like this graph here: I have come pretty…

---

## [A lagged table-join for finding the last ten rows that meets a criterion, for each row](https://discourse.julialang.org/t/a-lagged-table-join-for-finding-the-last-ten-rows-that-meets-a-criterion-for-each-row/108352)

<div class="topic-metadata">

**Author:** [@Cong](https://discourse.julialang.org/u/Cong)\
**Replies:** 6\
**Last updated:** [January 4, 2024, 10:54pm UTC](https://discourse.julialang.org/t/a-lagged-table-join-for-finding-the-last-ten-rows-that-meets-a-criterion-for-each-row/108352 "2024-01-04T22:54:45Z")

</div>

I have a dataframe like below. I am trying to add 6 columns to it: Out1Bac1, Out1Bac2, Out1Bac3, Out2Bac1, Out2Bac2, Out2Bac3. For each row, find the last row whose abschoice equals the CP1 of the current row, and get t…

---

## [ABM - Scheduler.bytype just for 1 agent function](https://discourse.julialang.org/t/abm-scheduler-bytype-just-for-1-agent-function/107611)

<div class="topic-metadata">

**Author:** [@elisadonati23](https://discourse.julialang.org/u/elisadonati23)\
**Replies:** 2\
**Last updated:** [January 4, 2024, 5:25pm UTC](https://discourse.julialang.org/t/abm-scheduler-bytype-just-for-1-agent-function/107611 "2024-01-04T17:25:25Z")

</div>

Dear all, I am running an ABM model with NoSpaceAgent representing sardines populations over time with three lifestages identified by the field type::Symbol and with 2 sexes (sex::Symbol). Adults can spawn in a period …

---

## [ModelingToolkit model with structural\_simplify error: How to initialize variables correctly?](https://discourse.julialang.org/t/modelingtoolkit-model-with-structural-simplify-error-how-to-initialize-variables-correctly/108110)

<div class="topic-metadata">

**Author:** [@juanmvenegas](https://discourse.julialang.org/u/juanmvenegas)\
**Replies:** 5\
**Last updated:** [January 4, 2024, 2:49pm UTC](https://discourse.julialang.org/t/modelingtoolkit-model-with-structural-simplify-error-how-to-initialize-variables-correctly/108110 "2024-01-04T14:49:41Z")

</div>

Hello everybody, I am currently developing a ModelingToolkit.jl model to simulate axial flow in a tubular pipe with a chemical reaction, including axial diffusion. This system is currently an initial value problem, and …

---

## [JuMP: Good first issue](https://discourse.julialang.org/t/jump-good-first-issue/108313)

<div class="topic-metadata">

**Author:** [@Paulnkk](https://discourse.julialang.org/u/Paulnkk)\
**Replies:** 3\
**Last updated:** [January 4, 2024, 12:05pm UTC](https://discourse.julialang.org/t/jump-good-first-issue/108313 "2024-01-04T12:05:05Z")

</div>

Hey, I am new in the jump community and wanted to learn more about it due to various applications for my phd research. I have some developer experience in cvxpy and happy to learn more! Here I wanted to ask if there are …

---

## [Calculate Equilibrium points for continuous system DynamicalSystems.jl](https://discourse.julialang.org/t/calculate-equilibrium-points-for-continuous-system-dynamicalsystems-jl/85247)

<div class="topic-metadata">

**Author:** [@Sergey\_Novak](https://discourse.julialang.org/u/Sergey_Novak)\
**Replies:** 9\
**Last updated:** [January 4, 2024, 4:50am UTC](https://discourse.julialang.org/t/calculate-equilibrium-points-for-continuous-system-dynamicalsystems-jl/85247 "2024-01-04T04:50:49Z")

</div>

I need calculate all equilibrium points for dynamical system Hindmarsh Rose (HR). Code for system: function HR!(du, u, p, t) function sigma(x) return 1.0 / ( 1.0 + exp( -10.0 \* ( x - ( - 0.25 ) ) )…

---

## [Lux Loss Not Decreasing](https://discourse.julialang.org/t/lux-loss-not-decreasing/108277)

<div class="topic-metadata">

**Author:** [@Dale\_James\_Black](https://discourse.julialang.org/u/Dale_James_Black)\
**Replies:** 1\
**Last updated:** [January 4, 2024, 12:18am UTC](https://discourse.julialang.org/t/lux-loss-not-decreasing/108277 "2024-01-04T00:18:42Z")

</div>

I am currently developing a comprehensive tutorial for my lab and the broader Julia community on training a Lux model for image segmentation. This tutorial is contained within a Pluto notebook that encompasses the entire…

---

## [Beta-Binomial Bayesian Model - larger samples result in more uncertainty?](https://discourse.julialang.org/t/beta-binomial-bayesian-model-larger-samples-result-in-more-uncertainty/108314)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 1\
**Last updated:** [January 3, 2024, 6:43pm UTC](https://discourse.julialang.org/t/beta-binomial-bayesian-model-larger-samples-result-in-more-uncertainty/108314 "2024-01-03T18:43:23Z")

</div>

I have a question about the Beta-Binomial Bayesian model. Let’s say I’d like to implement some quality control using this method. Based on a previous review of 10 widgets, I start with the prior belief that 90% of widget…

---

## [Incorrect output when using Base.Threads](https://discourse.julialang.org/t/incorrect-output-when-using-base-threads/108290)

<div class="topic-metadata">

**Author:** [@Uranium238](https://discourse.julialang.org/u/Uranium238)\
**Replies:** 3\
**Last updated:** [January 3, 2024, 4:14pm UTC](https://discourse.julialang.org/t/incorrect-output-when-using-base-threads/108290 "2024-01-03T16:14:46Z")

</div>

I am trying to write a code which has a computationally expensive triply nested loop. I am trying to parallelize the top lop using the @threads macro from Base.Threads. The loops involve calling to two functions GFC and G…

---

## [Is there a package that supports gamut calculations for the xyY color space?](https://discourse.julialang.org/t/is-there-a-package-that-supports-gamut-calculations-for-the-xyy-color-space/108128)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 4\
**Last updated:** [January 3, 2024, 1:52pm UTC](https://discourse.julialang.org/t/is-there-a-package-that-supports-gamut-calculations-for-the-xyy-color-space/108128 "2024-01-03T13:52:48Z")

</div>

Is there a Julia package that allows checking if a color in XYZ or xyY space is valid, and if not, allow finding a nearest point (eg by bisecting the ray from the white point)?

---

## [Failed test of CUDA package](https://discourse.julialang.org/t/failed-test-of-cuda-package/107378)

<div class="topic-metadata">

**Author:** [@nihil](https://discourse.julialang.org/u/nihil)\
**Replies:** 6\
**Last updated:** [January 3, 2024, 8:27am UTC](https://discourse.julialang.org/t/failed-test-of-cuda-package/107378 "2024-01-03T08:27:19Z")

</div>

There are a fair few failed tests from ‘test CUDA’ in the package manager. Is this normal, or is this from a miss-match of driver and julia pagage? CUDA.versioninfo() CUDA runtime 12.3, artifact installation CUDA dri…

---

## [How to get location of sequence in GenomicAnnotations](https://discourse.julialang.org/t/how-to-get-location-of-sequence-in-genomicannotations/106731)

<div class="topic-metadata">

**Author:** [@michikazjp](https://discourse.julialang.org/u/michikazjp)\
**Replies:** 2\
**Last updated:** [January 3, 2024, 3:02am UTC](https://discourse.julialang.org/t/how-to-get-location-of-sequence-in-genomicannotations/106731 "2024-01-03T03:02:54Z")

</div>

Hello I am newcomer in BioJulia and getting interested in it. Now I am starting to use GenomicAnnotations package. For example, I open a genbank file of E. coli K12 as using GenomicAnnotations K12 = readgbk(“K12.gbff”…

---

## [Piecewise linear regression with automatic knot selection](https://discourse.julialang.org/t/piecewise-linear-regression-with-automatic-knot-selection/93984)

<div class="topic-metadata">

**Author:** [@meghapatnaik](https://discourse.julialang.org/u/meghapatnaik)\
**Replies:** 9\
**Last updated:** [January 2, 2024, 8:42pm UTC](https://discourse.julialang.org/t/piecewise-linear-regression-with-automatic-knot-selection/93984 "2024-01-02T20:42:04Z")

</div>

Hi everyone, I am new to Julia, mostly work with R. I am looking for the equivalent of the “segmented” package in R, to fit a piecewise linear (constant in my special case) regression where the knot location is automate…

---

## [Safety Download files](https://discourse.julialang.org/t/safety-download-files/108233)

<div class="topic-metadata">

**Author:** [@CeterisPartybus](https://discourse.julialang.org/u/CeterisPartybus)\
**Replies:** 3\
**Last updated:** [January 2, 2024, 4:19pm UTC](https://discourse.julialang.org/t/safety-download-files/108233 "2024-01-02T16:19:41Z")

</div>

This might be a stupid question as I am not a network/security guy. I am thinking about implementing a feature in a package I develop that downloads a file from an external server. Should safety be a concern here? Spec…

---

## [Access variables of another subsystem in a callback](https://discourse.julialang.org/t/access-variables-of-another-subsystem-in-a-callback/108248)

<div class="topic-metadata">

**Author:** [@HKruenaegel](https://discourse.julialang.org/u/HKruenaegel)\
**Replies:** 1\
**Last updated:** [January 2, 2024, 8:23am UTC](https://discourse.julialang.org/t/access-variables-of-another-subsystem-in-a-callback/108248 "2024-01-02T08:23:13Z")

</div>

Is it possible to read a variable of another subsystem in a callback? For example if I want tu use decay2.u in a callback defined in decay1 in the following example code: using ModelingToolkit function affect!(integ,u,p…

---

## [Reduce(logaddexp,x) has arbitrarily better accuracy than logsumexp(x)](https://discourse.julialang.org/t/reduce-logaddexp-x-has-arbitrarily-better-accuracy-than-logsumexp-x/108229)

<div class="topic-metadata">

**Author:** [@cdawg](https://discourse.julialang.org/u/cdawg)\
**Replies:** 11\
**Last updated:** [January 2, 2024, 4:01am UTC](https://discourse.julialang.org/t/reduce-logaddexp-x-has-arbitrarily-better-accuracy-than-logsumexp-x/108229 "2024-01-02T04:01:31Z")

</div>

It occurred to me this morning that logsumexp is associative and commutative and can be calculated with a simple reduce. We could write logsumexp(x) = reduce(logaddexp,x) Isn’t that pretty! So how does that compare t…

---

## [PreallocationTools.jl with nested ForwardDiff and sparsity pattern detection errors](https://discourse.julialang.org/t/preallocationtools-jl-with-nested-forwarddiff-and-sparsity-pattern-detection-errors/107897)

<div class="topic-metadata">

**Author:** [@asterycs](https://discourse.julialang.org/u/asterycs)\
**Replies:** 6\
**Last updated:** [January 1, 2024, 9:05pm UTC](https://discourse.julialang.org/t/preallocationtools-jl-with-nested-forwarddiff-and-sparsity-pattern-detection-errors/107897 "2024-01-01T21:05:57Z")

</div>

Hello! I’m trying to combine sparsity pattern detection with a nested ForwardDiff call together with PreallocationTools. A contrived, minimal example is below: using ForwardDiff using PreallocationTools using Symbolics…

---

## [How to tweak dodge bar plot's label position of AOG](https://discourse.julialang.org/t/how-to-tweak-dodge-bar-plots-label-position-of-aog/108221)

<div class="topic-metadata">

**Author:** [@math4mad](https://discourse.julialang.org/u/math4mad)\
**Replies:** 0\
**Last updated:** [January 1, 2024, 1:26am UTC](https://discourse.julialang.org/t/how-to-tweak-dodge-bar-plots-label-position-of-aog/108221 "2024-01-01T01:26:33Z")

</div>

I want to clone this chart : Views about AI’s impact on society in the next 20 years, 2021 here is my code using Tidier, DataFrames, CSV, AlgebraOfGraphics, Pipe countrys = \["Japan", "China", "Germany", "France"\] c…

---

## [SIGVIEW - interface with DAQ via JULIA](https://discourse.julialang.org/t/sigview-interface-with-daq-via-julia/108202)

<div class="topic-metadata">

**Author:** [@Jawad\_Malik](https://discourse.julialang.org/u/Jawad_Malik)\
**Replies:** 8\
**Last updated:** [December 31, 2023, 5:46pm UTC](https://discourse.julialang.org/t/sigview-interface-with-daq-via-julia/108202 "2023-12-31T17:46:51Z")

</div>

I am working on a Condition Based Monitoring System for a small Marine Vessel. Need a FFT Software at the end of the line to keep a check on the Diesel Engines. I have four sets of Piezo + DAQ to pick up vibrations; 24…

---

## [Advice on Command Line Usage on HPC with Julia](https://discourse.julialang.org/t/advice-on-command-line-usage-on-hpc-with-julia/108176)

<div class="topic-metadata">

**Author:** [@TI36XPro](https://discourse.julialang.org/u/TI36XPro)\
**Replies:** 6\
**Last updated:** [December 31, 2023, 3:50am UTC](https://discourse.julialang.org/t/advice-on-command-line-usage-on-hpc-with-julia/108176 "2023-12-31T03:50:50Z")

</div>

Hello, I have been working and re-working the method I use to generate and submit runs on my HPC and I’d appreciate any advice or comments on the process (I’m sure there is plenty of room for improvement). Apologies if …

---

## [Convex.jl, possible to have a product as an operation?](https://discourse.julialang.org/t/convex-jl-possible-to-have-a-product-as-an-operation/108101)

<div class="topic-metadata">

**Author:** [@1m1](https://discourse.julialang.org/u/1m1)\
**Replies:** 8\
**Last updated:** [December 30, 2023, 5:21pm UTC](https://discourse.julialang.org/t/convex-jl-possible-to-have-a-product-as-an-operation/108101 "2023-12-30T17:21:17Z")

</div>

My constraints have a product. It is a convex set. I cannot find an operation to simply have the product of the components of a vector (or matrix/tensor in general). Is that not possible?

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

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