# Specific Domains

**URL:** https://discourse.julialang.org/c/domain/10.md?page=242

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

**Page:** 243

---

## [Siamese network in Lux.jl / Re-using parts of network](https://discourse.julialang.org/t/siamese-network-in-lux-jl-re-using-parts-of-network/82013)

<div class="topic-metadata">

**Author:** [@svilupp](https://discourse.julialang.org/u/svilupp)\
**Replies:** 10\
**Last updated:** [October 24, 2022, 4:46pm UTC](https://discourse.julialang.org/t/siamese-network-in-lux-jl-re-using-parts-of-network/82013 "2022-10-24T16:46:38Z")

</div>

Hi everyone! Have you seen any examples of re-using parts of the networks in Lux.jl? I haven’t been able to find anything. I’ve been trying to put a toy example of Siamese network (ala Keras - Siamese Contrastive Loss,…

---

## [Lux and Flux GPU function definitions overlap](https://discourse.julialang.org/t/lux-and-flux-gpu-function-definitions-overlap/89149)

<div class="topic-metadata">

**Author:** [@erlebach](https://discourse.julialang.org/u/erlebach)\
**Replies:** 4\
**Last updated:** [October 24, 2022, 1:39pm UTC](https://discourse.julialang.org/t/lux-and-flux-gpu-function-definitions-overlap/89149 "2022-10-24T13:39:13Z")

</div>

I am running the NeuralODE example for GPUs. Here is the code: using DifferentialEquations, Flux, DiffEqFlux, SciMLSensitivity using Random rng = Random.default\_rng() model\_gpu = Chain(Dense(2, 50, tanh), Dense(50, 2)…

---

## [How to give names to the curves in 1D graph of 2D array?](https://discourse.julialang.org/t/how-to-give-names-to-the-curves-in-1d-graph-of-2d-array/89182)

<div class="topic-metadata">

**Author:** [@ryofurue](https://discourse.julialang.org/u/ryofurue)\
**Replies:** 1\
**Last updated:** [October 24, 2022, 11:25am UTC](https://discourse.julialang.org/t/how-to-give-names-to-the-curves-in-1d-graph-of-2d-array/89182 "2022-10-24T11:25:33Z")

</div>

In my M x N array, the columns represent y values and so plot(arr) will produce N curves, which is what I need. But my question is, how can one give a name to each curve? My first intuition was to give a vector of the n…

---

## [Awkward ODE problem](https://discourse.julialang.org/t/awkward-ode-problem/88435)

<div class="topic-metadata">

**Author:** [@josuagrw](https://discourse.julialang.org/u/josuagrw)\
**Replies:** 7\
**Last updated:** [October 24, 2022, 11:14am UTC](https://discourse.julialang.org/t/awkward-ode-problem/88435 "2022-10-24T11:14:02Z")

</div>

I have a differential equation of the kind \\dot q\_1 = f(q\_1, q\_2)\\\\ A(q\_1,q\_2)\\dot q\_2=B(q\_1,q\_2)\\dot q\_1 where q\_1, q\_2 are vectors, A is a symmetric, positive-definite real matrix, and B is a real (non-square) matri…

---

## [Analysis and Diagnostics for MCMC](https://discourse.julialang.org/t/analysis-and-diagnostics-for-mcmc/3104)

<div class="topic-metadata">

**Author:** [@jacobcvt12](https://discourse.julialang.org/u/jacobcvt12)\
**Replies:** 27\
**Last updated:** [October 24, 2022, 10:05am UTC](https://discourse.julialang.org/t/analysis-and-diagnostics-for-mcmc/3104 "2022-10-24T10:05:44Z")

</div>

Is there a package in Julia similar to the coda package in R? Most of the functions in coda are pretty straightforward, but it’s nice to have a bunch of MCMC diagnostics/analysis functions all grouped together. I see tha…

---

## [Question about quadratic cost function in JuMP](https://discourse.julialang.org/t/question-about-quadratic-cost-function-in-jump/89143)

<div class="topic-metadata">

**Author:** [@Steve\_Zhang](https://discourse.julialang.org/u/Steve_Zhang)\
**Replies:** 2\
**Last updated:** [October 24, 2022, 1:56am UTC](https://discourse.julialang.org/t/question-about-quadratic-cost-function-in-jump/89143 "2022-10-24T01:56:52Z")

</div>

I am new to Julia and transferred from Matlab. Now I am using JuMP to write an optimization problem with a quadratic cost function like: and my code for optimization problem is: function get\_cost\_func(sys, R, M) …

---

## [PATHsolver not working](https://discourse.julialang.org/t/pathsolver-not-working/89156)

<div class="topic-metadata">

**Author:** [@Muhammad\_Nadeem](https://discourse.julialang.org/u/Muhammad_Nadeem)\
**Replies:** 3\
**Last updated:** [October 23, 2022, 11:25pm UTC](https://discourse.julialang.org/t/pathsolver-not-working/89156 "2022-10-23T23:25:00Z")

</div>

Hello all, I am trying to run a simple mixed complementary problem using PATHsolver, but it is giving me this error “UndefVarError: solveMCP not defined” My code is: using PATHSolver M = \[0 0 -1 -1 ; 0 0 1 -2 ; …

---

## [Unsupported Argument in Optimization.solve](https://discourse.julialang.org/t/unsupported-argument-in-optimization-solve/89146)

<div class="topic-metadata">

**Author:** [@erlebach](https://discourse.julialang.org/u/erlebach)\
**Replies:** 1\
**Last updated:** [October 23, 2022, 5:29pm UTC](https://discourse.julialang.org/t/unsupported-argument-in-optimization-solve/89146 "2022-10-23T17:29:48Z")

</div>

I am running a demonstration of NeuralODE, fond at: https://diffeqflux.sciml.ai/stable/examples/neural\_ode/: result\_neuralode2 = Optimization.solve(optprob2, Optim.BFGS(initial\_stepnorm=0.01), callback = callb…

---

## [Error running highdim\_pde script from 2019 on 2022 SciML versions](https://discourse.julialang.org/t/error-running-highdim-pde-script-from-2019-on-2022-sciml-versions/89093)

<div class="topic-metadata">

**Author:** [@erlebach](https://discourse.julialang.org/u/erlebach)\
**Replies:** 12\
**Last updated:** [October 23, 2022, 3:33pm UTC](https://discourse.julialang.org/t/error-running-highdim-pde-script-from-2019-on-2022-sciml-versions/89093 "2022-10-23T15:33:50Z")

</div>

I am trying to run the “Universal Differential Equations” demos posted at the github repository: GitHub - ChrisRackauckas/universal\_differential\_equations: Repository for the Universal Differential Equations for Scienti…

---

## [How to obtain maximized ELBO value from Turing to use as approximation to model evidence?](https://discourse.julialang.org/t/how-to-obtain-maximized-elbo-value-from-turing-to-use-as-approximation-to-model-evidence/88787)

<div class="topic-metadata">

**Author:** [@neuro\_enthusiast](https://discourse.julialang.org/u/neuro_enthusiast)\
**Replies:** 5\
**Last updated:** [October 23, 2022, 11:02am UTC](https://discourse.julialang.org/t/how-to-obtain-maximized-elbo-value-from-turing-to-use-as-approximation-to-model-evidence/88787 "2022-10-23T11:02:33Z")

</div>

Hi all! I’m running variational inference in Turing. However, I can’t seem to find a way to extract the maximized ELBO from the Turing ADVI. If anyone can point me to relevant information that would be wonderful!

---

## [Connection timed out problem on SSHManager of Distributed.jl](https://discourse.julialang.org/t/connection-timed-out-problem-on-sshmanager-of-distributed-jl/89129)

<div class="topic-metadata">

**Author:** [@hsianktin](https://discourse.julialang.org/u/hsianktin)\
**Replies:** 0\
**Last updated:** [October 23, 2022, 5:51am UTC](https://discourse.julialang.org/t/connection-timed-out-problem-on-sshmanager-of-distributed-jl/89129 "2022-10-23T05:51:57Z")

</div>

Hi everyone, I met a problem in adding another Windows machine as a worker using SSHManager. Background information: I am able to direct ssh into the remote machine, and the machine is under the same local network wit…

---

## [Wasserstein distance for histograms](https://discourse.julialang.org/t/wasserstein-distance-for-histograms/89083)

<div class="topic-metadata">

**Author:** [@ignace](https://discourse.julialang.org/u/ignace)\
**Replies:** 2\
**Last updated:** [October 22, 2022, 11:38pm UTC](https://discourse.julialang.org/t/wasserstein-distance-for-histograms/89083 "2022-10-22T23:38:53Z")

</div>

Hello, How can one compute the Wasserstein distance between two histograms. To be clear, what I mean with a histogram is, e.g., the following h : using StatsBase x = randn(10^4) h = fit(Histogram, x) Looking forwar…

---

## [ABLines in Algebraofgraphics (x=y line, identity line, 1-to-1 line)](https://discourse.julialang.org/t/ablines-in-algebraofgraphics-x-y-line-identity-line-1-to-1-line/89106)

<div class="topic-metadata">

**Author:** [@mreichMPI-BGC](https://discourse.julialang.org/u/mreichMPI-BGC)\
**Replies:** 2\
**Last updated:** [October 22, 2022, 7:39pm UTC](https://discourse.julialang.org/t/ablines-in-algebraofgraphics-x-y-line-identity-line-1-to-1-line/89106 "2022-10-22T19:39:28Z")

</div>

Dear all, I am new to Julia and AoG - coming from the Tidyverse and ggplot… - and don’t (even) manage to add an x=y line to a scatter plot: using AlgebraOfGraphics, CairoMakie, Statistics x=rand(4) y = x .+ 0.2 .+ 0.2…

---

## [Obtaining the posterior distribution of parameters in Turing using the MCMC method](https://discourse.julialang.org/t/obtaining-the-posterior-distribution-of-parameters-in-turing-using-the-mcmc-method/89111)

<div class="topic-metadata">

**Author:** [@Aru](https://discourse.julialang.org/u/Aru)\
**Replies:** 3\
**Last updated:** [October 22, 2022, 6:35pm UTC](https://discourse.julialang.org/t/obtaining-the-posterior-distribution-of-parameters-in-turing-using-the-mcmc-method/89111 "2022-10-22T18:35:37Z")

</div>

Hello everyone. I am a student majoring in physics with no solid background in statistics. I recently tried to use Turing, but I found that there are some hard-to-understand problems that make it difficult for me to use …

---

## [JuMP constraint with limited dimension of an index](https://discourse.julialang.org/t/jump-constraint-with-limited-dimension-of-an-index/89100)

<div class="topic-metadata">

**Author:** [@MeryKay](https://discourse.julialang.org/u/MeryKay)\
**Replies:** 1\
**Last updated:** [October 22, 2022, 9:57am UTC](https://discourse.julialang.org/t/jump-constraint-with-limited-dimension-of-an-index/89100 "2022-10-22T09:57:15Z")

</div>

Hi everyboday! In one constraint, I want to limit my index to a specific intervall within the intervall which is used in all other cases. More specific this means: index p is normally defined in p = 1:P whereby P=10 N…

---

## [Optimizing noisy objective](https://discourse.julialang.org/t/optimizing-noisy-objective/25260)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 60\
**Last updated:** [October 22, 2022, 8:19am UTC](https://discourse.julialang.org/t/optimizing-noisy-objective/25260 "2022-10-22T08:19:10Z")

</div>

I have a f: \\mathbb{R}^n \\to \\mathbb{R}^m function, m \> n, and I am trying to minimize some norm of f(x). The norm itself is flexible, I would like to get something working. The problem is that f is stochastic as it com…

---

## [Makie workflow to access and change attributes of objects](https://discourse.julialang.org/t/makie-workflow-to-access-and-change-attributes-of-objects/89058)

<div class="topic-metadata">

**Author:** [@weltenbummler](https://discourse.julialang.org/u/weltenbummler)\
**Replies:** 4\
**Last updated:** [October 21, 2022, 8:05pm UTC](https://discourse.julialang.org/t/makie-workflow-to-access-and-change-attributes-of-objects/89058 "2022-10-21T20:05:55Z")

</div>

Hi, This post is about plotting with Makie. I quite often set out making a figure and then later add further plots or labels. My workflow below works but seems a bit cumbersome to me and I would like to have your input …

---

## [How many flops does it take to compute a square root?](https://discourse.julialang.org/t/how-many-flops-does-it-take-to-compute-a-square-root/89027)

<div class="topic-metadata">

**Author:** [@nvenkov1](https://discourse.julialang.org/u/nvenkov1)\
**Replies:** 6\
**Last updated:** [October 21, 2022, 5:20pm UTC](https://discourse.julialang.org/t/how-many-flops-does-it-take-to-compute-a-square-root/89027 "2022-10-21T17:20:27Z")

</div>

How many flops does it take to compute a square root?

---

## [MILP: formulation of sums with barriers](https://discourse.julialang.org/t/milp-formulation-of-sums-with-barriers/89067)

<div class="topic-metadata">

**Author:** [@MeryKay](https://discourse.julialang.org/u/MeryKay)\
**Replies:** 4\
**Last updated:** [October 21, 2022, 5:17pm UTC](https://discourse.julialang.org/t/milp-formulation-of-sums-with-barriers/89067 "2022-10-21T17:17:48Z")

</div>

I want to program a sum over an indexed variable, whereby the index has a barrier just for this constraint. sum(X\[t\] for t=1:T) is clear to me, but how can I integrate the barrier of t\<= R\_j as shown in the screenshot. …

---

## [ChainRulesCore and ForwardDiff](https://discourse.julialang.org/t/chainrulescore-and-forwarddiff/61705)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 15\
**Last updated:** [October 21, 2022, 4:14pm UTC](https://discourse.julialang.org/t/chainrulescore-and-forwarddiff/61705 "2022-10-21T16:14:20Z")

</div>

If I want to define derivatives for a function to work with ForwardDiff, is it sufficient to define a ChainRulesCore.frule, or is there something else? The rule does not seem to be used, instead the original function is…

---

## [Nonlinear semidefinite programming recommendation](https://discourse.julialang.org/t/nonlinear-semidefinite-programming-recommendation/75338)

<div class="topic-metadata">

**Author:** [@wujinq](https://discourse.julialang.org/u/wujinq)\
**Replies:** 17\
**Last updated:** [October 21, 2022, 2:38pm UTC](https://discourse.julialang.org/t/nonlinear-semidefinite-programming-recommendation/75338 "2022-10-21T14:38:10Z")

</div>

I want to do the following optimization: \\min E \\ \\ \\text{s.t.} \\ \\ x \\in \\mathbb{R}, \\ \\ \\text{$M$ is semidefinite}, where M\_{ij} is a polynominal of x and E. Is there a Julia package suitable for this problem?

---

## [CallWithMetadata has no field val error for simple PDE system in ModelingToolkit](https://discourse.julialang.org/t/callwithmetadata-has-no-field-val-error-for-simple-pde-system-in-modelingtoolkit/89055)

<div class="topic-metadata">

**Author:** [@hbooth](https://discourse.julialang.org/u/hbooth)\
**Replies:** 2\
**Last updated:** [October 21, 2022, 1:02pm UTC](https://discourse.julialang.org/t/callwithmetadata-has-no-field-val-error-for-simple-pde-system-in-modelingtoolkit/89055 "2022-10-21T13:02:26Z")

</div>

Hi using DifferentialEquations using ModelingToolkit, MethodOfLines @parameters x @variables t c(..) Dx = Differential(x) Dxx = Differential(x)^2 Dt = Differential(t) L = 1. D = 0.16 k = 1. j0 = 0.4 eq = Dt(c(t,x))…

---

## [Optimizing over the quantum entropy cone natively in JuMP](https://discourse.julialang.org/t/optimizing-over-the-quantum-entropy-cone-natively-in-jump/59833)

<div class="topic-metadata">

**Author:** [@ryancorywright](https://discourse.julialang.org/u/ryancorywright)\
**Replies:** 6\
**Last updated:** [October 21, 2022, 12:33pm UTC](https://discourse.julialang.org/t/optimizing-over-the-quantum-entropy-cone-natively-in-jump/59833 "2022-10-21T12:33:31Z")

</div>

Hi everyone, I was wondering if it was possible to optimize over the quantum entropy cone natively in JuMP. If you need a reminder, the cone is \\{(X\_1, X\_2, X\_3): X\_1 \\preceq X\_2^{\\frac{1}{2}} \\log(X\_2^{\\frac{-1}{2}}X…

---

## [Changing line color Theme in CairoMakie and GLMakie](https://discourse.julialang.org/t/changing-line-color-theme-in-cairomakie-and-glmakie/88990)

<div class="topic-metadata">

**Author:** [@garrek](https://discourse.julialang.org/u/garrek)\
**Replies:** 7\
**Last updated:** [October 21, 2022, 9:40am UTC](https://discourse.julialang.org/t/changing-line-color-theme-in-cairomakie-and-glmakie/88990 "2022-10-21T09:40:22Z")

</div>

I’m trying to define a custom plotting theme in Makie.jl. I want to remove the default color cycle for Lines and have black be the default color. I’m not sure if my syntax is correct, but it works in CairoMakie. But, whe…

---

## [How to integrate out latent variables/ intermediate values in Turing model](https://discourse.julialang.org/t/how-to-integrate-out-latent-variables-intermediate-values-in-turing-model/89003)

<div class="topic-metadata">

**Author:** [@neuro\_enthusiast](https://discourse.julialang.org/u/neuro_enthusiast)\
**Replies:** 4\
**Last updated:** [October 21, 2022, 9:27am UTC](https://discourse.julialang.org/t/how-to-integrate-out-latent-variables-intermediate-values-in-turing-model/89003 "2022-10-21T09:27:05Z")

</div>

Hi all, my generative model has the following form. P(y |\\theta) = \\int P(y|x) P(x|\\theta) dx I want to get the posterior of θ given fixed data y. But I unfortunately don’t have an analytical solution for the inte…

---

## [Matrix SDEs](https://discourse.julialang.org/t/matrix-sdes/89038)

<div class="topic-metadata">

**Author:** [@simoneb](https://discourse.julialang.org/u/simoneb)\
**Replies:** 2\
**Last updated:** [October 21, 2022, 6:20am UTC](https://discourse.julialang.org/t/matrix-sdes/89038 "2022-10-21T06:20:47Z")

</div>

Hi all, I am trying to fit the following matrix SDE using the following code: G = \[0 -1; 1 0\]; p = (A=G, B=G) function drift(du, u, p, t) ## drift function for the SDE du = t \* p.A; end # drift function diffusion…

---

## [Help with Arrow.jl and size of files](https://discourse.julialang.org/t/help-with-arrow-jl-and-size-of-files/89028)

<div class="topic-metadata">

**Author:** [@austin-putz](https://discourse.julialang.org/u/austin-putz)\
**Replies:** 23\
**Last updated:** [October 21, 2022, 2:40am UTC](https://discourse.julialang.org/t/help-with-arrow-jl-and-size-of-files/89028 "2022-10-21T02:40:20Z")

</div>

I was recently exploring the Arrow.jl package and arrow package in R. I had a project where I needed to write a function to convert a large file in Julia, then port it to R for another package we need to use there. I wr…

---

## [Reducing framework overhead of solving a simple JuMP model in a hot loop](https://discourse.julialang.org/t/reducing-framework-overhead-of-solving-a-simple-jump-model-in-a-hot-loop/88977)

<div class="topic-metadata">

**Author:** [@erwinkn](https://discourse.julialang.org/u/erwinkn)\
**Replies:** 4\
**Last updated:** [October 20, 2022, 7:34pm UTC](https://discourse.julialang.org/t/reducing-framework-overhead-of-solving-a-simple-jump-model-in-a-hot-loop/88977 "2022-10-20T19:34:57Z")

</div>

I have a nonlinear JuMP model that is used in a hot loop: it is run 4000-8000 times, each run with a different set of parameters. The reason for this loop is that I want to perform global sensitivity analysis on the para…

---

## [Parallel AMG](https://discourse.julialang.org/t/parallel-amg/88926)

<div class="topic-metadata">

**Author:** [@nvenkov1](https://discourse.julialang.org/u/nvenkov1)\
**Replies:** 8\
**Last updated:** [October 20, 2022, 4:28pm UTC](https://discourse.julialang.org/t/parallel-amg/88926 "2022-10-20T16:28:26Z")

</div>

Has anyone ever started working on a parallel implementation of algebraic multigrid solvers and preconditioners?

---

## [Strategies to optimize expensive multi-heatmap animation](https://discourse.julialang.org/t/strategies-to-optimize-expensive-multi-heatmap-animation/89019)

<div class="topic-metadata">

**Author:** [@aramirezreyes](https://discourse.julialang.org/u/aramirezreyes)\
**Replies:** 0\
**Last updated:** [October 20, 2022, 4:06pm UTC](https://discourse.julialang.org/t/strategies-to-optimize-expensive-multi-heatmap-animation/89019 "2022-10-20T16:06:08Z")

</div>

crossposting from Makie discord Hi all! Warning, no MWE for the moment due to the nature of the problem. This is a question about strategy. How to animate this multi-heatmap figure? I am using CairoMakie as I do this …

[Previous page](https://discourse.julialang.org/c/domain/10.md?page=241)

[Next page](https://discourse.julialang.org/c/domain/10.md?page=243)