# Specific Domains

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

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

**Page:** 102

---

## [Why isn’t Symbolic Regression used more?](https://discourse.julialang.org/t/why-isn-t-symbolic-regression-used-more/116215)

<div class="topic-metadata">

**Author:** [@Tarny\_GG\_Channie](https://discourse.julialang.org/u/Tarny_GG_Channie)\
**Replies:** 10\
**Last updated:** [June 27, 2024, 6:36pm UTC](https://discourse.julialang.org/t/why-isn-t-symbolic-regression-used-more/116215 "2024-06-27T18:36:07Z")

</div>

Thanks to its ability to dynamically compile expression and execute it efficiently, Julia positions itself very strongly in the symbolic regression landscape. Symbolic regression offers a way to discover formulas with f…

---

## [Bayesian parameter estimation of a ODE with Bernouilli multivariate likelihood](https://discourse.julialang.org/t/bayesian-parameter-estimation-of-a-ode-with-bernouilli-multivariate-likelihood/116303)

<div class="topic-metadata">

**Author:** [@Marta\_Pardo\_Araujo](https://discourse.julialang.org/u/Marta_Pardo_Araujo)\
**Replies:** 6\
**Last updated:** [June 27, 2024, 2:44pm UTC](https://discourse.julialang.org/t/bayesian-parameter-estimation-of-a-ode-with-bernouilli-multivariate-likelihood/116303 "2024-06-27T14:44:03Z")

</div>

Hi, I am trying to estimate the parameters of a ODE model with N equations. The model’s output are probabilities, ie, a vector of N probabilities each time step. Each coordinate of the vector defines the probability tha…

---

## [Tolerance in using JuMP optimizer](https://discourse.julialang.org/t/tolerance-in-using-jump-optimizer/116143)

<div class="topic-metadata">

**Author:** [@gustave\_Robichon](https://discourse.julialang.org/u/gustave_Robichon)\
**Replies:** 2\
**Last updated:** [June 25, 2024, 9:14am UTC](https://discourse.julialang.org/t/tolerance-in-using-jump-optimizer/116143 "2024-06-25T09:14:15Z")

</div>

Hello :slight\_smile: I am using JuMP, and a solver named Hypatia. I want to change the tol\_inconsistent (the tolerance to numerical imprecision) to 1e-12. I found on the web that I should use this : julia \> MOI.set(mo…

---

## [Using anonymous and or symbols to create dynamic variable names?](https://discourse.julialang.org/t/using-anonymous-and-or-symbols-to-create-dynamic-variable-names/83095)

<div class="topic-metadata">

**Author:** [@Ryan\_Vogt1](https://discourse.julialang.org/u/Ryan_Vogt1)\
**Replies:** 7\
**Last updated:** [June 27, 2024, 8:53am UTC](https://discourse.julialang.org/t/using-anonymous-and-or-symbols-to-create-dynamic-variable-names/83095 "2024-06-27T08:53:48Z")

</div>

Howdy, so I asked a question on stack overflow:julia - Define several JuMP variables names in loop based on iteration - Stack Overflow , and Oscar replied. I thought a bit more, and remember this board, and thought I’d s…

---

## [How to save algorithm iteration result when using JuMP](https://discourse.julialang.org/t/how-to-save-algorithm-iteration-result-when-using-jump/115684)

<div class="topic-metadata">

**Author:** [@Xuekui\_wang](https://discourse.julialang.org/u/Xuekui_wang)\
**Replies:** 6\
**Last updated:** [June 27, 2024, 12:35am UTC](https://discourse.julialang.org/t/how-to-save-algorithm-iteration-result-when-using-jump/115684 "2024-06-27T00:35:53Z")

</div>

When using Ipopt in JuMP，it will print these information on console, I want to know how to get these info, for example, save these data as a vector or dict. thanks. '' iter objective inf\_pr inf\_du lg(mu) ||d||…

---

## [Struggle at large scale problem of nlp using ipopt](https://discourse.julialang.org/t/struggle-at-large-scale-problem-of-nlp-using-ipopt/116240)

<div class="topic-metadata">

**Author:** [@Chi\_Jan](https://discourse.julialang.org/u/Chi_Jan)\
**Replies:** 10\
**Last updated:** [June 26, 2024, 10:47pm UTC](https://discourse.julialang.org/t/struggle-at-large-scale-problem-of-nlp-using-ipopt/116240 "2024-06-26T22:47:57Z")

</div>

hi there, I have a non-linear problem where the decision variable is a matrix of dimensions m \* n \* k. The constraints may also be somewhat complex. When I set m to a small subset of the dataset, e.g., 10, the program se…

---

## [User defined gradient in JuMP with Ipopt: gradient with user data](https://discourse.julialang.org/t/user-defined-gradient-in-jump-with-ipopt-gradient-with-user-data/116236)

<div class="topic-metadata">

**Author:** [@vict\_riag](https://discourse.julialang.org/u/vict_riag)\
**Replies:** 3\
**Last updated:** [June 26, 2024, 10:15pm UTC](https://discourse.julialang.org/t/user-defined-gradient-in-jump-with-ipopt-gradient-with-user-data/116236 "2024-06-26T22:15:56Z")

</div>

Hello, I have been trying to register a gradient in JuMP that includes user data. I think the issue is with how I have formulated the gradient and fed it into JuMP, but I am not sure how to use the @expression or @param…

---

## [Problem with ForwardDiff and a custom vector type in OrdinaryDiffEq](https://discourse.julialang.org/t/problem-with-forwarddiff-and-a-custom-vector-type-in-ordinarydiffeq/116260)

<div class="topic-metadata">

**Author:** [@JonasWickman](https://discourse.julialang.org/u/JonasWickman)\
**Replies:** 5\
**Last updated:** [June 26, 2024, 8:07pm UTC](https://discourse.julialang.org/t/problem-with-forwarddiff-and-a-custom-vector-type-in-ordinarydiffeq/116260 "2024-06-26T20:07:58Z")

</div>

I have an ODE system tracking the densities of a number of species. If certain criteria are fulfilled I add more species to the system and run the new slightly bigger ODE system, and repeat this process until certain oth…

---

## [Do a function like relu need a kernel ? When you need to write a GPU kernel rather than "just" using CuArray?](https://discourse.julialang.org/t/do-a-function-like-relu-need-a-kernel-when-you-need-to-write-a-gpu-kernel-rather-than-just-using-cuarray/116265)

<div class="topic-metadata">

**Author:** [@sylvaticus](https://discourse.julialang.org/u/sylvaticus)\
**Replies:** 3\
**Last updated:** [June 26, 2024, 4:54pm UTC](https://discourse.julialang.org/t/do-a-function-like-relu-need-a-kernel-when-you-need-to-write-a-gpu-kernel-rather-than-just-using-cuarray/116265 "2024-06-26T16:54:03Z")

</div>

The relu function works elementwise by returning the input that is non-negative, i.e. relu(x) = ifelse.(x .\> 0, x, 0). Does the fact that it works elementwise means I need to write a GPU kernel for it, or I can simply u…

---

## [GraphPPL RxInfer Model](https://discourse.julialang.org/t/graphppl-rxinfer-model/114255)

<div class="topic-metadata">

**Author:** [@ohmsweetohm1](https://discourse.julialang.org/u/ohmsweetohm1)\
**Replies:** 3\
**Last updated:** [June 26, 2024, 4:44pm UTC](https://discourse.julialang.org/t/graphppl-rxinfer-model/114255 "2024-06-26T16:44:04Z")

</div>

I have tried to extend the example from here to include an unknown parameter theta: @model function rotate\_ssm(y, x0, Q, P) x\_prior ~ MvNormalMeanCovariance(mean(x0), cov(x0)) x\_prev = x\_prior theta ~ Norma…

---

## [Does MixedModels assume a normal distribution on coefficient estimates?](https://discourse.julialang.org/t/does-mixedmodels-assume-a-normal-distribution-on-coefficient-estimates/116152)

<div class="topic-metadata">

**Author:** [@ergative](https://discourse.julialang.org/u/ergative)\
**Replies:** 1\
**Last updated:** [June 26, 2024, 3:44pm UTC](https://discourse.julialang.org/t/does-mixedmodels-assume-a-normal-distribution-on-coefficient-estimates/116152 "2024-06-26T15:44:10Z")

</div>

Hi–new to Julia here, sorry if this has been answered. I came over from lme4 in R-land, where you get p-values from lmers from the lmerTest package, which uses the Satterthwaite method for estimating degrees of freedom a…

---

## [Using Integrator in EventQueueABM model](https://discourse.julialang.org/t/using-integrator-in-eventqueueabm-model/116249)

<div class="topic-metadata">

**Author:** [@Sahil\_Khan](https://discourse.julialang.org/u/Sahil_Khan)\
**Replies:** 0\
**Last updated:** [June 26, 2024, 9:14am UTC](https://discourse.julialang.org/t/using-integrator-in-eventqueueabm-model/116249 "2024-06-26T09:14:30Z")

</div>

Hello everyone, I want to include some differential equation based model to my agents model. I want to integrate both of them. so at some points in simulation time, agent model get feed some info from differential model…

---

## [Writing AD rules for Rotations.jl](https://discourse.julialang.org/t/writing-ad-rules-for-rotations-jl/116222)

<div class="topic-metadata">

**Author:** [@sonofwatt](https://discourse.julialang.org/u/sonofwatt)\
**Replies:** 1\
**Last updated:** [June 26, 2024, 5:14am UTC](https://discourse.julialang.org/t/writing-ad-rules-for-rotations-jl/116222 "2024-06-26T05:14:46Z")

</div>

I’m trying to create a coordinate transformation, via LinearAlgebra and Rotations.jl (Home · Rotations.jl), and can be optimized via Zygote.jl, or the like. I’m running into a problem because some of the functions in Ro…

---

## [How is the parameter p\_steps evaluated in Morris method?](https://discourse.julialang.org/t/how-is-the-parameter-p-steps-evaluated-in-morris-method/116124)

<div class="topic-metadata">

**Author:** [@Nick1](https://discourse.julialang.org/u/Nick1)\
**Replies:** 3\
**Last updated:** [June 26, 2024, 5:03am UTC](https://discourse.julialang.org/t/how-is-the-parameter-p-steps-evaluated-in-morris-method/116124 "2024-06-26T05:03:24Z")

</div>

It is not clear from the tutorial how the p\_steps parameter is evaluated in Morris method. Usually, the parameters that are defined are the total number of trajectories and the number of trajectories with the highest spr…

---

## [2D interpolation of large data to calculate depth profiles](https://discourse.julialang.org/t/2d-interpolation-of-large-data-to-calculate-depth-profiles/116190)

<div class="topic-metadata">

**Author:** [@tamasgal](https://discourse.julialang.org/u/tamasgal)\
**Replies:** 32\
**Last updated:** [June 25, 2024, 10:14pm UTC](https://discourse.julialang.org/t/2d-interpolation-of-large-data-to-calculate-depth-profiles/116190 "2024-06-25T22:14:11Z")

</div>

I am working with a dataset which is essentially a depth sample map: x and y coordinates on an irregular grid and the corresponding depth values. I need to calculate tons of depth profiles for different lines defined by …

---

## [Injecting known primal solution (JuMP + Gurobi)](https://discourse.julialang.org/t/injecting-known-primal-solution-jump-gurobi/116213)

<div class="topic-metadata">

**Author:** [@HenriDeh](https://discourse.julialang.org/u/HenriDeh)\
**Replies:** 1\
**Last updated:** [June 25, 2024, 6:28pm UTC](https://discourse.julialang.org/t/injecting-known-primal-solution-jump-gurobi/116213 "2024-06-25T18:28:36Z")

</div>

Hi, Is there a way to provide a solver (specifically Gurobi) with a feasible primal solution to help it prune the tree ? The answer is probably somewhere in JuMP’s documentation but the keywords I used to search lead t…

---

## [Two y-axis - Issue with limits](https://discourse.julialang.org/t/two-y-axis-issue-with-limits/116205)

<div class="topic-metadata">

**Author:** [@NaSi](https://discourse.julialang.org/u/NaSi)\
**Replies:** 1\
**Last updated:** [June 25, 2024, 8:34pm UTC](https://discourse.julialang.org/t/two-y-axis-issue-with-limits/116205 "2024-06-25T20:34:30Z")

</div>

Hi :slight\_smile: I am using GLMakie and I would like to add a second yaxis. I have some troubles for defining the limits of both y-axis. The axis seem dependent (when I fix the limit of ax2, it changes the vizualisatio…

---

## [How to display moving objects in Makie.jl](https://discourse.julialang.org/t/how-to-display-moving-objects-in-makie-jl/54282)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 14\
**Last updated:** [June 25, 2024, 7:36pm UTC](https://discourse.julialang.org/t/how-to-display-moving-objects-in-makie-jl/54282 "2024-06-25T19:36:43Z")

</div>

Hello, I am now able to display a 3D scene in the way I want it. But how can I display a moving system? My main loop looks like this: for i = 0:4 # calculate a vector of 3D coordinates X = range(0,…

---

## [Incremental Technical Analysis indicators](https://discourse.julialang.org/t/incremental-technical-analysis-indicators/107844)

<div class="topic-metadata">

**Author:** [@FemtoTrader](https://discourse.julialang.org/u/FemtoTrader)\
**Replies:** 23\
**Last updated:** [June 25, 2024, 11:00am UTC](https://discourse.julialang.org/t/incremental-technical-analysis-indicators/107844 "2024-06-25T11:00:44Z")

</div>

Hi, Unlike existing libraries for technical analysis which typically have to work on the whole input vector in order to calculate new values of indicators I’m looking for a Julia library which implements technical analy…

---

## [StatsBase crosscor standard deviation does not have lags applied](https://discourse.julialang.org/t/statsbase-crosscor-standard-deviation-does-not-have-lags-applied/116182)

<div class="topic-metadata">

**Author:** [@BG-marine](https://discourse.julialang.org/u/BG-marine)\
**Replies:** 1\
**Last updated:** [June 25, 2024, 5:47am UTC](https://discourse.julialang.org/t/statsbase-crosscor-standard-deviation-does-not-have-lags-applied/116182 "2024-06-25T05:47:57Z")

</div>

In StatsBase.crosscor!() why is the standard deviation for normalisation calculated before the lag is applied to the input vectors, rather than after? This doesn’t appear to follow a standard correlation formula. Cross…

---

## [Slow cubature](https://discourse.julialang.org/t/slow-cubature/116120)

<div class="topic-metadata">

**Author:** [@e3c6](https://discourse.julialang.org/u/e3c6)\
**Replies:** 21\
**Last updated:** [June 25, 2024, 12:46am UTC](https://discourse.julialang.org/t/slow-cubature/116120 "2024-06-25T00:46:35Z")

</div>

I’m trying the perform a numerical integration over 4-dimensions. Here is a self-contained example: import HCubature using StaticArrays: SMatrix using StaticArrays: SVector using LinearAlgebra: Symmetric using LinearAlg…

---

## [Logical conditions](https://discourse.julialang.org/t/logical-conditions/115363)

<div class="topic-metadata">

**Author:** [@papar](https://discourse.julialang.org/u/papar)\
**Replies:** 14\
**Last updated:** [June 25, 2024, 12:40am UTC](https://discourse.julialang.org/t/logical-conditions/115363 "2024-06-25T00:40:18Z")

</div>

I’m working on an optimization problem in JuMP where I need to ensure that certain constraints are enforced based on a binary variable. Here are the variables and constraints I have defined in my model: Variable Definit…

---

## [Distributed computing or Multi-threading on Optimization](https://discourse.julialang.org/t/distributed-computing-or-multi-threading-on-optimization/116147)

<div class="topic-metadata">

**Author:** [@benjamin.nunez](https://discourse.julialang.org/u/benjamin.nunez)\
**Replies:** 1\
**Last updated:** [June 24, 2024, 7:05pm UTC](https://discourse.julialang.org/t/distributed-computing-or-multi-threading-on-optimization/116147 "2024-06-24T19:05:21Z")

</div>

Hello! I am currently working on an algorithm to solve stochastic optimization problems by decomposing them into scenarios (very similar to Progressive Hedging). I am not very knowledgeable in the area of task parallel…

---

## [Structural SImplify Unable to perform broadcast equation](https://discourse.julialang.org/t/structural-simplify-unable-to-perform-broadcast-equation/115967)

<div class="topic-metadata">

**Author:** [@HSnowden](https://discourse.julialang.org/u/HSnowden)\
**Replies:** 1\
**Last updated:** [June 24, 2024, 3:48pm UTC](https://discourse.julialang.org/t/structural-simplify-unable-to-perform-broadcast-equation/115967 "2024-06-24T15:48:06Z")

</div>

Hello, In a model I habe been building, I go to perform structural simplify upon it and it results in the following error MethodError: no method matching length(::SymbolicUtils.BasicSymbolic{Any}). The code to blame fo…

---

## [Moving Horizon Estimation](https://discourse.julialang.org/t/moving-horizon-estimation/115921)

<div class="topic-metadata">

**Author:** [@MarcoDellOro](https://discourse.julialang.org/u/MarcoDellOro)\
**Replies:** 8\
**Last updated:** [June 24, 2024, 11:39am UTC](https://discourse.julialang.org/t/moving-horizon-estimation/115921 "2024-06-24T11:39:08Z")

</div>

Hello everybody, I am having some issues with the Moving Horizon Estimation (MHE) code I am writing. To me, the task and the code seem quite easy, but I cannot obtain the desired results for some reason. I attach here …

---

## [A question about PSF deconvolution](https://discourse.julialang.org/t/a-question-about-psf-deconvolution/116108)

<div class="topic-metadata">

**Author:** [@Cyan](https://discourse.julialang.org/u/Cyan)\
**Replies:** 7\
**Last updated:** [June 24, 2024, 10:15am UTC](https://discourse.julialang.org/t/a-question-about-psf-deconvolution/116108 "2024-06-24T10:15:17Z")

</div>

Could I obtain the deconvolved image if I have an original image and its PSF? And if I have a deconvolved image and its PSF, how can I obtain the convolved image?

---

## [Media properties and component encapsulation in MTK](https://discourse.julialang.org/t/media-properties-and-component-encapsulation-in-mtk/83557)

<div class="topic-metadata">

**Author:** [@dodoplus](https://discourse.julialang.org/u/dodoplus)\
**Replies:** 6\
**Last updated:** [June 24, 2024, 9:36am UTC](https://discourse.julialang.org/t/media-properties-and-component-encapsulation-in-mtk/83557 "2024-06-24T09:36:29Z")

</div>

Hi, Modelica has a Media library that captures the properties of various materials. It implements the equations given in IF97 (a version of was apparently implemented in Julia). It provides e.g. calculation of the Cp, C…

---

## [\`sum\` syntax rule differs between Julia and \`JuMP\`](https://discourse.julialang.org/t/sum-syntax-rule-differs-between-julia-and-jump/116125)

<div class="topic-metadata">

**Author:** [@WuSiren](https://discourse.julialang.org/u/WuSiren)\
**Replies:** 2\
**Last updated:** [June 23, 2024, 11:59pm UTC](https://discourse.julialang.org/t/sum-syntax-rule-differs-between-julia-and-jump/116125 "2024-06-23T23:59:51Z")

</div>

a = \[1, 2, 3, 4\] sum(a\[n\] for n in 1:0) # ERROR: MethodError: reducing over an empty collection is not allowed; consider supplying \`init\` to the reducer using JuMP m = Model() @constraint(m, sum(a\[n\] for n in 1:0) \>= 0)…

---

## [How to ensure Makie's theme activation occurs when called in a module?](https://discourse.julialang.org/t/how-to-ensure-makies-theme-activation-occurs-when-called-in-a-module/116114)

<div class="topic-metadata">

**Author:** [@Datseris](https://discourse.julialang.org/u/Datseris)\
**Replies:** 2\
**Last updated:** [June 23, 2024, 6:27pm UTC](https://discourse.julialang.org/t/how-to-ensure-makies-theme-activation-occurs-when-called-in-a-module/116114 "2024-06-23T18:27:05Z")

</div>

I had some script with lots of convenience plotting functions for Makie, as well as a bunch of themes I like to use in my work. For reproducibility I decided to make it a Julia package so that it can be added and its git…

---

## [Convergence of gmres from IterativeSolvers in minimal examples](https://discourse.julialang.org/t/convergence-of-gmres-from-iterativesolvers-in-minimal-examples/115852)

<div class="topic-metadata">

**Author:** [@Pablo\_Marchant](https://discourse.julialang.org/u/Pablo_Marchant)\
**Replies:** 7\
**Last updated:** [June 23, 2024, 2:55pm UTC](https://discourse.julialang.org/t/convergence-of-gmres-from-iterativesolvers-in-minimal-examples/115852 "2024-06-23T14:55:07Z")

</div>

Hi all, I’m trying to see if I can apply gmres in a domain specific problem I have, where I need to solve a system Ax=b where A is a non-symmetric square block tridiagonal matrix. I have not managed to get this to work,…

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

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