# Specific Domains

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

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

**Page:** 82

---

## [Plot the confidence interval for a model fit](https://discourse.julialang.org/t/plot-the-confidence-interval-for-a-model-fit/37767)

<div class="topic-metadata">

**Author:** [@ffevotte](https://discourse.julialang.org/u/ffevotte)\
**Replies:** 20\
**Last updated:** [September 24, 2021, 12:43pm UTC](https://discourse.julialang.org/t/plot-the-confidence-interval-for-a-model-fit/37767 "2021-09-24T12:43:53Z")

</div>

I’m no data scientist, so it’s likely that my question does not make much sense… But here goes: I recently tried to do some data analysis that involved fitting models. I’m able to fit the model, compute represent a pred…

---

## [Avoid storing intermediate results from the forward pass by default](https://discourse.julialang.org/t/avoid-storing-intermediate-results-from-the-forward-pass-by-default/119694)

<div class="topic-metadata">

**Author:** [@Jakub\_Mitura](https://discourse.julialang.org/u/Jakub_Mitura)\
**Replies:** 14\
**Last updated:** [September 29, 2024, 4:10pm UTC](https://discourse.julialang.org/t/avoid-storing-intermediate-results-from-the-forward-pass-by-default/119694 "2024-09-29T16:10:01Z")

</div>

Hello, Is it possible in Zygote to change the default behavior for computing derivatives during backpropagation to use gradient checkpointing? I have a memory-constrained problem with a Lux.jl model that uses Zygote fo…

---

## [SymbolicRegression -- two simple questions](https://discourse.julialang.org/t/symbolicregression-two-simple-questions/119908)

<div class="topic-metadata">

**Author:** [@BLI](https://discourse.julialang.org/u/BLI)\
**Replies:** 2\
**Last updated:** [September 29, 2024, 8:09am UTC](https://discourse.julialang.org/t/symbolicregression-two-simple-questions/119908 "2024-09-29T08:09:24Z")

</div>

I figured out the answers.

---

## [Having a hard time with Makie Recipes](https://discourse.julialang.org/t/having-a-hard-time-with-makie-recipes/119978)

<div class="topic-metadata">

**Author:** [@co1emi11er](https://discourse.julialang.org/u/co1emi11er)\
**Replies:** 7\
**Last updated:** [September 28, 2024, 1:45pm UTC](https://discourse.julialang.org/t/having-a-hard-time-with-makie-recipes/119978 "2024-09-28T13:45:27Z")

</div>

Current Struggle I am currently building an extension for both Makie.jl and Plots.jl. While I have managed to get it working in Plots.jl, I have not had the same pleasant experience with Makie. I keep coming back to it a…

---

## [Do we need separate package for OwensQ function (power calk)?](https://discourse.julialang.org/t/do-we-need-separate-package-for-owensq-function-power-calk/120346)

<div class="topic-metadata">

**Author:** [@PharmCat](https://discourse.julialang.org/u/PharmCat)\
**Replies:** 0\
**Last updated:** [September 28, 2024, 1:06pm UTC](https://discourse.julialang.org/t/do-we-need-separate-package-for-owensq-function-power-calk/120346 "2024-09-28T13:06:40Z")

</div>

Do we need separate package for OwensQ function (power calk) ? So this is the question. I’m using OwensQ function for sample size estimation for BE study in ClinicalTrialUtilities.jl , but this package was made when I …

---

## [Trying to serialize runtime generated function and MTK ODESystem](https://discourse.julialang.org/t/trying-to-serialize-runtime-generated-function-and-mtk-odesystem/119961)

<div class="topic-metadata">

**Author:** [@Bart\_van\_de\_Lint](https://discourse.julialang.org/u/Bart_van_de_Lint)\
**Replies:** 5\
**Last updated:** [September 28, 2024, 11:52am UTC](https://discourse.julialang.org/t/trying-to-serialize-runtime-generated-function-and-mtk-odesystem/119961 "2024-09-28T11:52:37Z")

</div>

I have a system with around 1500 equations. It takes a long time to run the following function, so I want to serialize and deserialize f\_ip, dvs, psym and io\_sys so I can use them between sessions. (\_, f\_ip), dvs, psym,…

---

## [Discrete Fourier Transform (DFT) matrix and inverse](https://discourse.julialang.org/t/discrete-fourier-transform-dft-matrix-and-inverse/119952)

<div class="topic-metadata">

**Author:** [@Nikos\_Gianniotis](https://discourse.julialang.org/u/Nikos_Gianniotis)\
**Replies:** 18\
**Last updated:** [September 28, 2024, 7:02am UTC](https://discourse.julialang.org/t/discrete-fourier-transform-dft-matrix-and-inverse/119952 "2024-09-28T07:02:40Z")

</div>

Hello everyone, I want to calculate the explicit DFT matrix (and its inverse). I have been looking at related Julia packages to see if this functionality is somewhere available. While I could calculate it myself, I am a…

---

## [Nested and different AD methods altogether: How to add AD calculations inside my loss function when using neural differential equations?](https://discourse.julialang.org/t/nested-and-different-ad-methods-altogether-how-to-add-ad-calculations-inside-my-loss-function-when-using-neural-differential-equations/108985)

<div class="topic-metadata">

**Author:** [@facusapienza](https://discourse.julialang.org/u/facusapienza)\
**Replies:** 9\
**Last updated:** [September 28, 2024, 12:16am UTC](https://discourse.julialang.org/t/nested-and-different-ad-methods-altogether-how-to-add-ad-calculations-inside-my-loss-function-when-using-neural-differential-equations/108985 "2024-09-28T00:16:46Z")

</div>

Hi all, I am implementing regularization penalties inside Universal Differential Equations (also applicable to Physics-Informed neural networks) where I need to differentiate a (loss) function that includes in its calcu…

---

## [Using Juniper with Convex.jl](https://discourse.julialang.org/t/using-juniper-with-convex-jl/119974)

<div class="topic-metadata">

**Author:** [@CaG21](https://discourse.julialang.org/u/CaG21)\
**Replies:** 2\
**Last updated:** [September 27, 2024, 9:45pm UTC](https://discourse.julialang.org/t/using-juniper-with-convex-jl/119974 "2024-09-27T21:45:23Z")

</div>

I am trying to solve a MISOCP, modeled using Convex.jl, and I thought that trying to solve it with Juniper would be great since it is one of the non-commercial solvers that can handle this problem. However, when I call s…

---

## [Efficiently check if points are contained in polygons](https://discourse.julialang.org/t/efficiently-check-if-points-are-contained-in-polygons/74415)

<div class="topic-metadata">

**Author:** [@ahmoreira](https://discourse.julialang.org/u/ahmoreira)\
**Replies:** 33\
**Last updated:** [September 27, 2024, 9:15pm UTC](https://discourse.julialang.org/t/efficiently-check-if-points-are-contained-in-polygons/74415 "2024-09-27T21:15:34Z")

</div>

Hi all, I am new to Julia and I am trying to check if points (plant occurrence points) fall within polygons from a shapefile. The code I wrote works, but I would like to make it more efficient. Currently, the code takes …

---

## [Best Practice for Type Declarations in CUDA Kernels](https://discourse.julialang.org/t/best-practice-for-type-declarations-in-cuda-kernels/119968)

<div class="topic-metadata">

**Author:** [@jamesmulhern](https://discourse.julialang.org/u/jamesmulhern)\
**Replies:** 3\
**Last updated:** [September 27, 2024, 4:40pm UTC](https://discourse.julialang.org/t/best-practice-for-type-declarations-in-cuda-kernels/119968 "2024-09-27T16:40:10Z")

</div>

I am writing a number of CUDA kernels for my optimization project. In the rest of my code I am using parametric type declarations so that the code can be run using different types. An example of an existing function is…

---

## [Search for eigenvalues under a specific upper boundary](https://discourse.julialang.org/t/search-for-eigenvalues-under-a-specific-upper-boundary/119913)

<div class="topic-metadata">

**Author:** [@tianyi\_zhang](https://discourse.julialang.org/u/tianyi_zhang)\
**Replies:** 6\
**Last updated:** [September 27, 2024, 2:58pm UTC](https://discourse.julialang.org/t/search-for-eigenvalues-under-a-specific-upper-boundary/119913 "2024-09-27T14:58:27Z")

</div>

Hello, I am trying to solve an eigenvalue problem in which eigenvalues lesser than a upper boundary are solved. The linear operator is hermitian or real symmetric, and the eigenvalues are guaranteed to be positive, so I…

---

## [EAGO unbounded variables](https://discourse.julialang.org/t/eago-unbounded-variables/119865)

<div class="topic-metadata">

**Author:** [@kiki\_Van\_Elsen](https://discourse.julialang.org/u/kiki_Van_Elsen)\
**Replies:** 3\
**Last updated:** [September 27, 2024, 7:52am UTC](https://discourse.julialang.org/t/eago-unbounded-variables/119865 "2024-09-27T07:52:39Z")

</div>

Hello! I want to run the following code: function bivariate(solver, type) try #choose which solver you want to use m = Model(solver.Optimizer) #function we want to Optimize f(x1,x2)…

---

## [Building a custom "variable in set" constraint](https://discourse.julialang.org/t/building-a-custom-variable-in-set-constraint/119927)

<div class="topic-metadata">

**Author:** [@janneshb](https://discourse.julialang.org/u/janneshb)\
**Replies:** 9\
**Last updated:** [September 26, 2024, 11:02pm UTC](https://discourse.julialang.org/t/building-a-custom-variable-in-set-constraint/119927 "2024-09-26T23:02:14Z")

</div>

Hi, I implemented my own type of sets (let’s call it MySetType) and now I want to override Julia’s in functionality, similar to what is already supported for SecondOrderCones. So, I would like to make this work with Ju…

---

## [Turing models, vectors, and type stability?](https://discourse.julialang.org/t/turing-models-vectors-and-type-stability/119928)

<div class="topic-metadata">

**Author:** [@rkube](https://discourse.julialang.org/u/rkube)\
**Replies:** 0\
**Last updated:** [September 26, 2024, 7:04pm UTC](https://discourse.julialang.org/t/turing-models-vectors-and-type-stability/119928 "2024-09-26T19:04:06Z")

</div>

Hi, I’m puzzled by a type instability issue in the Turing model below. When the model defines the vector p2 before iterating, the model is type stable. When using the \[...\] notation, the model is unstable. using Linear…

---

## [How to approach my unconstrained optimization problem?](https://discourse.julialang.org/t/how-to-approach-my-unconstrained-optimization-problem/119919)

<div class="topic-metadata">

**Author:** [@nsajko](https://discourse.julialang.org/u/nsajko)\
**Replies:** 4\
**Last updated:** [September 26, 2024, 6:36pm UTC](https://discourse.julialang.org/t/how-to-approach-my-unconstrained-optimization-problem/119919 "2024-09-26T18:36:32Z")

</div>

I’ve got an objective function to maximize. The domain is {\\mathbb{R}}^n, where n is about thirty. The codomain is a bounded subset of the nonnegative reals. The objective value is zero for most of the domain. The object…

---

## [Non-intuitive difference of copyto! between CuArray and Array](https://discourse.julialang.org/t/non-intuitive-difference-of-copyto-between-cuarray-and-array/119904)

<div class="topic-metadata">

**Author:** [@0samuraiE](https://discourse.julialang.org/u/0samuraiE)\
**Replies:** 3\
**Last updated:** [September 26, 2024, 2:43pm UTC](https://discourse.julialang.org/t/non-intuitive-difference-of-copyto-between-cuarray-and-array/119904 "2024-09-26T14:43:01Z")

</div>

Hi all. I trapped an error when using copyto! to CuArray. Error is ERROR: CUDA error: invalid argument (code 1, ERROR\_INVALID\_VALUE) Stacktrace: \[1\] throw\_api\_error(res::CUDA.cudaError\_enum) @ CUDA ~/.julia/packa…

---

## [Confusion with variable indices with JuMP/CPLEX](https://discourse.julialang.org/t/confusion-with-variable-indices-with-jump-cplex/119920)

<div class="topic-metadata">

**Author:** [@mike\_k](https://discourse.julialang.org/u/mike_k)\
**Replies:** 1\
**Last updated:** [September 26, 2024, 2:35pm UTC](https://discourse.julialang.org/t/confusion-with-variable-indices-with-jump-cplex/119920 "2024-09-26T14:35:15Z")

</div>

I try to set branching priorities on variables using CPXcopyorder(env, lp, cnt, indices, priority, direction). Particlularly I want to get the indices of the varialbes of the JuMP model in the CPLEX model. I tried: inde…

---

## [Differentiable Regularizers for CUDA, CPUs?](https://discourse.julialang.org/t/differentiable-regularizers-for-cuda-cpus/119915)

<div class="topic-metadata">

**Author:** [@roflmaostc](https://discourse.julialang.org/u/roflmaostc)\
**Replies:** 0\
**Last updated:** [September 26, 2024, 12:12pm UTC](https://discourse.julialang.org/t/differentiable-regularizers-for-cuda-cpus/119915 "2024-09-26T12:12:55Z")

</div>

Hi, is anyone aware of a package which implements Total Variation (TV) or other kinds of regularizer efficiently and differentiable? A naive, differentiable way of TV would be: julia\> arr = rand(10, 10) 10×10 Matrix{F…

---

## [Parameters in the neural network not updating after training](https://discourse.julialang.org/t/parameters-in-the-neural-network-not-updating-after-training/119914)

<div class="topic-metadata">

**Author:** [@Ashima\_Kalathingal](https://discourse.julialang.org/u/Ashima_Kalathingal)\
**Replies:** 0\
**Last updated:** [September 26, 2024, 12:04pm UTC](https://discourse.julialang.org/t/parameters-in-the-neural-network-not-updating-after-training/119914 "2024-09-26T12:04:53Z")

</div>

I am implementing a Neural ODE in julia. When I complete the training using ADAM optimizer, the loss function decreases as expected. But after training, the parameters reset to the original initial value. So I am not abl…

---

## [How to write a new gradient descent algorithm and benchmark against existing methods](https://discourse.julialang.org/t/how-to-write-a-new-gradient-descent-algorithm-and-benchmark-against-existing-methods/119869)

<div class="topic-metadata">

**Author:** [@cdawg](https://discourse.julialang.org/u/cdawg)\
**Replies:** 4\
**Last updated:** [September 26, 2024, 11:06am UTC](https://discourse.julialang.org/t/how-to-write-a-new-gradient-descent-algorithm-and-benchmark-against-existing-methods/119869 "2024-09-26T11:06:07Z")

</div>

I have an idea for a new gradient descent algorithm. I just want to see if it does better or worse than a standard momentum methods. How would I go about making use of the existing benchmarking suites CUTEst.jl or Optim…

---

## [Bonito slider emit event if mousedrag stops](https://discourse.julialang.org/t/bonito-slider-emit-event-if-mousedrag-stops/119897)

<div class="topic-metadata">

**Author:** [@baxmittens](https://discourse.julialang.org/u/baxmittens)\
**Replies:** 0\
**Last updated:** [September 26, 2024, 8:09am UTC](https://discourse.julialang.org/t/bonito-slider-emit-event-if-mousedrag-stops/119897 "2024-09-26T08:09:32Z")

</div>

Hi there, I want something like this: using Bonito using WGLMakie import Bonito.TailwindDashboard as D app = App() do session::Session sliders = \[Bonito.Slider(-1:.1:1, value=0.0) for i = 1:6\] vals = map!(Observable…

---

## [How to add a manually created legend inside the axis?](https://discourse.julialang.org/t/how-to-add-a-manually-created-legend-inside-the-axis/119893)

<div class="topic-metadata">

**Author:** [@WuSiren](https://discourse.julialang.org/u/WuSiren)\
**Replies:** 5\
**Last updated:** [September 26, 2024, 5:56am UTC](https://discourse.julialang.org/t/how-to-add-a-manually-created-legend-inside-the-axis/119893 "2024-09-26T05:56:27Z")

</div>

using CairoMakie f = Figure() Axis(f\[1, 1\]) elem\_1 = \[LineElement(color = :red, linestyle = nothing), MarkerElement(color = :blue, marker = 'x', markersize = 15, strokecolor = :black)\] elem\_2 = \[P…

---

## [Help finding multiple roots to nonlinear equation - IntervalRootFinding not working](https://discourse.julialang.org/t/help-finding-multiple-roots-to-nonlinear-equation-intervalrootfinding-not-working/119575)

<div class="topic-metadata">

**Author:** [@lpz18](https://discourse.julialang.org/u/lpz18)\
**Replies:** 13\
**Last updated:** [September 25, 2024, 8:42pm UTC](https://discourse.julialang.org/t/help-finding-multiple-roots-to-nonlinear-equation-intervalrootfinding-not-working/119575 "2024-09-25T20:42:44Z")

</div>

I have a set of nonlinear equations (which happen to be steady-state DEs) that I am trying to solve that have multiple real roots and 3 variables (x, y, & z in the code below) In Mathematica I can get the multiple roots…

---

## [Basinattractionfor4D](https://discourse.julialang.org/t/basinattractionfor4d/119793)

<div class="topic-metadata">

**Author:** [@Sarbast\_Hussein\_Mika](https://discourse.julialang.org/u/Sarbast_Hussein_Mika)\
**Replies:** 3\
**Last updated:** [September 25, 2024, 7:59pm UTC](https://discourse.julialang.org/t/basinattractionfor4d/119793 "2024-09-25T19:59:28Z")

</div>

Hello everyone Is there any code for Basin attraction for 4D nonliner differential system with line equilibrium? To classify regions Chaotic Periodic Line equilibrium Hope someone could help me in this regards

---

## [Good to copy GPU arrays back to CPU?](https://discourse.julialang.org/t/good-to-copy-gpu-arrays-back-to-cpu/119778)

<div class="topic-metadata">

**Author:** [@huiyuxie](https://discourse.julialang.org/u/huiyuxie)\
**Replies:** 6\
**Last updated:** [September 25, 2024, 7:51pm UTC](https://discourse.julialang.org/t/good-to-copy-gpu-arrays-back-to-cpu/119778 "2024-09-25T19:51:27Z")

</div>

I was wondering if it’s a good practice to always copy GPU arrays back to the CPU before a program ends. I’m running a large program that performs simulations with a large amount of data on the GPU. I know that some res…

---

## [Mismatched ForwardDiff and FiniteDifference gradients with ImplicitDifferentiation](https://discourse.julialang.org/t/mismatched-forwarddiff-and-finitedifference-gradients-with-implicitdifferentiation/119756)

<div class="topic-metadata">

**Author:** [@benjaminfaber](https://discourse.julialang.org/u/benjaminfaber)\
**Replies:** 12\
**Last updated:** [September 25, 2024, 7:31pm UTC](https://discourse.julialang.org/t/mismatched-forwarddiff-and-finitedifference-gradients-with-implicitdifferentiation/119756 "2024-09-25T19:31:05Z")

</div>

@gdalle I am attempting to use ImplicitDifferentiation with ForwardDiff to compute the gradient of an optimization problem. I’ve tried to figure out how to reduce the problem to a minimum working example, but I have not…

---

## [Change color of 'walls' in Axis3 in Makie](https://discourse.julialang.org/t/change-color-of-walls-in-axis3-in-makie/119882)

<div class="topic-metadata">

**Author:** [@JonasWickman](https://discourse.julialang.org/u/JonasWickman)\
**Replies:** 2\
**Last updated:** [September 25, 2024, 6:47pm UTC](https://discourse.julialang.org/t/change-color-of-walls-in-axis3-in-makie/119882 "2024-09-25T18:47:16Z")

</div>

Using Makie, I want to group some subfigures together with colored boxes as described here: How to draw boxes around subfigures | Makie. However, some of the subfigures are 3d plots, and the ‘walls’ and ‘floor’ of the Ax…

---

## [Switching stiff solvers when one is unstable](https://discourse.julialang.org/t/switching-stiff-solvers-when-one-is-unstable/119609)

<div class="topic-metadata">

**Author:** [@pulk\_jain](https://discourse.julialang.org/u/pulk_jain)\
**Replies:** 15\
**Last updated:** [September 25, 2024, 6:03pm UTC](https://discourse.julialang.org/t/switching-stiff-solvers-when-one-is-unstable/119609 "2024-09-25T18:03:59Z")

</div>

I am solving an ode (~40 variables) for many different initial conditions. Mostly TRBDF2 gives fastest results. But sometimes it shows that the system is unstable. In such cases I switched to Rodas5P or RadauIIA5 which g…

---

## [Solve Large Scale Underdetermined Linear Equation with per Element Equality Constraint](https://discourse.julialang.org/t/solve-large-scale-underdetermined-linear-equation-with-per-element-equality-constraint/119655)

<div class="topic-metadata">

**Author:** [@RoyiAvital](https://discourse.julialang.org/u/RoyiAvital)\
**Replies:** 16\
**Last updated:** [September 25, 2024, 5:10pm UTC](https://discourse.julialang.org/t/solve-large-scale-underdetermined-linear-equation-with-per-element-equality-constraint/119655 "2024-09-25T17:10:56Z")

</div>

I have the following linear system to solve: \\boldsymbol{A} \\boldsymbol{x} = \\boldsymbol{0}, \\quad \\text{subject to} \\; {x}\_{i} = {v}\_{i} \\; \\forall i \\in \\mathcal{V} Where \\boldsymbol{A} \\in \\mathbb{R}^{m \\times n} w…

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

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