# Specific Domains

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

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

**Page:** 222

---

## [Different behaviors on Latin Hypercube Sampling from two packages](https://discourse.julialang.org/t/different-behaviors-on-latin-hypercube-sampling-from-two-packages/93076)

<div class="topic-metadata">

**Author:** [@MJulia](https://discourse.julialang.org/u/MJulia)\
**Replies:** 8\
**Last updated:** [January 23, 2023, 11:51am UTC](https://discourse.julialang.org/t/different-behaviors-on-latin-hypercube-sampling-from-two-packages/93076 "2023-01-23T11:51:22Z")

</div>

I tried QuasiMonteCarlo.jl and LatinHypercubeSampling.jl. #for QuasiMonteCarlo.jl s1 = QuasiMonteCarlo.sample(25, \[0.0,0.0\], \[1.0,1.0\], LatinHypercubeSample()) 2×25 Matrix{Float64}: 0.4 0.56 0.96 1.0 0.68 0.16…

---

## [Frameworks, or libraries to use for feasibility problems?](https://discourse.julialang.org/t/frameworks-or-libraries-to-use-for-feasibility-problems/93114)

<div class="topic-metadata">

**Author:** [@rakshith95](https://discourse.julialang.org/u/rakshith95)\
**Replies:** 4\
**Last updated:** [January 23, 2023, 10:53am UTC](https://discourse.julialang.org/t/frameworks-or-libraries-to-use-for-feasibility-problems/93114 "2023-01-23T10:53:49Z")

</div>

Hello, I have a problem where I need to find the closest (min norm) point to my data which is feasible under some (non convex) polynomial constraints. What would be the best framework/ libraries to use for a problem lik…

---

## [Ignore, please (uninformed question)](https://discourse.julialang.org/t/ignore-please-uninformed-question/93390)

<div class="topic-metadata">

**Author:** [@mcreel](https://discourse.julialang.org/u/mcreel)\
**Replies:** 1\
**Last updated:** [January 23, 2023, 10:32am UTC](https://discourse.julialang.org/t/ignore-please-uninformed-question/93390 "2023-01-23T10:32:28Z")

</div>

In the course of preparing materials to explain boxplots to students, I am seeing something odd with the whiskers in a boxplot. The code ## example box X = randn(100) M = median(X) Q1 = quantile(X,0.25) Q3 = quantile(X,…

---

## [KNITRO.jl. Could not load library .dylib](https://discourse.julialang.org/t/knitro-jl-could-not-load-library-dylib/93332)

<div class="topic-metadata">

**Author:** [@nkashaev](https://discourse.julialang.org/u/nkashaev)\
**Replies:** 9\
**Last updated:** [January 22, 2023, 8:51pm UTC](https://discourse.julialang.org/t/knitro-jl-could-not-load-library-dylib/93332 "2023-01-22T20:51:11Z")

</div>

Hi, When testing KNITRO, I got the following error: could not load library “” dlopen(.dylib, 0x0001): tried: ‘/Users/XXX/knitro-12.3.0-MacOS-64/lib/.dylib’ (no such file), ‘.dylib’ (no such file), ‘/System/Volumes/Pr…

---

## [Functions embedded in Julia](https://discourse.julialang.org/t/functions-embedded-in-julia/93368)

<div class="topic-metadata">

**Author:** [@Andrea\_Del\_Bravo](https://discourse.julialang.org/u/Andrea_Del_Bravo)\
**Replies:** 2\
**Last updated:** [January 22, 2023, 6:17pm UTC](https://discourse.julialang.org/t/functions-embedded-in-julia/93368 "2023-01-22T18:17:24Z")

</div>

I installed Julia on a laptop running ferdora 36 operating system. I tried to capture some xeta function values but the answer was that zeta(x) was not recognize a a valid command. I read on the Julia manual that zeta is…

---

## [Understanding Lighting in Makie](https://discourse.julialang.org/t/understanding-lighting-in-makie/93356)

<div class="topic-metadata">

**Author:** [@Whyborn](https://discourse.julialang.org/u/Whyborn)\
**Replies:** 2\
**Last updated:** [January 22, 2023, 2:37pm UTC](https://discourse.julialang.org/t/understanding-lighting-in-makie/93356 "2023-01-22T14:37:56Z")

</div>

Hi, I’m having some problems understanding how the lighting works in Makie. I’m trying to draw a surface displaying a cone with a blunt leading edge, and then draw some flow fields extending from the surface of the cone…

---

## [Matrix Product NL expression](https://discourse.julialang.org/t/matrix-product-nl-expression/93345)

<div class="topic-metadata">

**Author:** [@apateonas](https://discourse.julialang.org/u/apateonas)\
**Replies:** 1\
**Last updated:** [January 22, 2023, 9:35am UTC](https://discourse.julialang.org/t/matrix-product-nl-expression/93345 "2023-01-22T09:35:32Z")

</div>

Hello, I am new to JuMP. I want to create the following non-linear expression: KdagK = @NLexpression(model, K' \* K) Here, K is a matrix of complex variables. What is the best way to do this? Thank you!

---

## [Frobenius Norm objective with JuMP](https://discourse.julialang.org/t/frobenius-norm-objective-with-jump/93338)

<div class="topic-metadata">

**Author:** [@apateonas](https://discourse.julialang.org/u/apateonas)\
**Replies:** 2\
**Last updated:** [January 22, 2023, 1:14am UTC](https://discourse.julialang.org/t/frobenius-norm-objective-with-jump/93338 "2023-01-22T01:14:47Z")

</div>

I am looking to use the frobenius norm of a matrix as an objective. Here is the setup: nkrauss = 4 model = Model(Ipopt.Optimizer) # Krauss operator variables nKelems = nkrauss \* dims^2 Kelems = \[@variable(model, set =…

---

## [JuMP solution is not correct: PiecewiseLinearOpt](https://discourse.julialang.org/t/jump-solution-is-not-correct-piecewiselinearopt/93321)

<div class="topic-metadata">

**Author:** [@salvaguerrero](https://discourse.julialang.org/u/salvaguerrero)\
**Replies:** 1\
**Last updated:** [January 22, 2023, 12:16am UTC](https://discourse.julialang.org/t/jump-solution-is-not-correct-piecewiselinearopt/93321 "2023-01-22T00:16:47Z")

</div>

Good morning, I’m working on a MILP model where a non-linear function is linearized using PiecewiseLinearOpt and Gurobi. The model is solved and the objective value reached is correct, but the constraints are not met. …

---

## [PyMNE for julia](https://discourse.julialang.org/t/pymne-for-julia/93315)

<div class="topic-metadata">

**Author:** [@skaarup91](https://discourse.julialang.org/u/skaarup91)\
**Replies:** 3\
**Last updated:** [January 21, 2023, 9:17pm UTC](https://discourse.julialang.org/t/pymne-for-julia/93315 "2023-01-21T21:17:31Z")

</div>

Hi. I am trying to import raw EEG data in the .cnt format to julia. I have downloaded PyMNE (the Python MNE package for julia) and I am trying with the PyMNE.io.read\_raw\_cnt function. I run: PyMNE.io.read\_raw\_cnt(“C:\\Us…

---

## [Trust-region constrained optimization](https://discourse.julialang.org/t/trust-region-constrained-optimization/93077)

<div class="topic-metadata">

**Author:** [@shce](https://discourse.julialang.org/u/shce)\
**Replies:** 32\
**Last updated:** [January 21, 2023, 8:14pm UTC](https://discourse.julialang.org/t/trust-region-constrained-optimization/93077 "2023-01-21T20:14:33Z")

</div>

Hello, I was wondering if there exists a Julia package that implements a trust-region constrained optimization method. The packages I have seen so far are for unconstrained or bound-constrained problems. I am looking fo…

---

## [Julia wrapper for SU2 (The Open-Source CFD Code) (Stanford University Unstructured)](https://discourse.julialang.org/t/julia-wrapper-for-su2-the-open-source-cfd-code-stanford-university-unstructured/93318)

<div class="topic-metadata">

**Author:** [@Mr.Benz](https://discourse.julialang.org/u/Mr.Benz)\
**Replies:** 0\
**Last updated:** [January 21, 2023, 3:40pm UTC](https://discourse.julialang.org/t/julia-wrapper-for-su2-the-open-source-cfd-code-stanford-university-unstructured/93318 "2023-01-21T15:40:20Z")

</div>

Python wrapper is available for SU2 (The Open-Source CFD Code). Please provide Julia wrapper for SU2. SU2 is a suite of open-source software tools written in C++ for the numerical solution of partial differential equat…

---

## [Rank 1 Update for LU factors](https://discourse.julialang.org/t/rank-1-update-for-lu-factors/93298)

<div class="topic-metadata">

**Author:** [@maxchendt](https://discourse.julialang.org/u/maxchendt)\
**Replies:** 0\
**Last updated:** [January 21, 2023, 8:08am UTC](https://discourse.julialang.org/t/rank-1-update-for-lu-factors/93298 "2023-01-21T08:08:56Z")

</div>

Is there a Julia implementation? as that from the MATLAB Central File Exchange:

---

## [Zygote error in backprop through NN](https://discourse.julialang.org/t/zygote-error-in-backprop-through-nn/92920)

<div class="topic-metadata">

**Author:** [@LucasMSpereira](https://discourse.julialang.org/u/LucasMSpereira)\
**Replies:** 4\
**Last updated:** [January 21, 2023, 1:54am UTC](https://discourse.julialang.org/t/zygote-error-in-backprop-through-nn/92920 "2023-01-21T01:54:27Z")

</div>

Following my last post about implementing WGAN-GP, now I’m running into another problem with gradients. Changing the discriminator’s BatchNorm() to Metalhead’s ChannelLayerNorm() solves the issue with the inner gradient.…

---

## [My model can't be transfered to gpu. Adapt.jl problem?](https://discourse.julialang.org/t/my-model-cant-be-transfered-to-gpu-adapt-jl-problem/92585)

<div class="topic-metadata">

**Author:** [@Euhan](https://discourse.julialang.org/u/Euhan)\
**Replies:** 12\
**Last updated:** [January 21, 2023, 1:45am UTC](https://discourse.julialang.org/t/my-model-cant-be-transfered-to-gpu-adapt-jl-problem/92585 "2023-01-21T01:45:00Z")

</div>

I have a model I’ve been working with and when I’m trying to expand it to work with arbitrary numbers of channels (some of which can be switched off) transfer to gpu keeps choking on the line: sparams = ntuple(i-\>F.para…

---

## [Plotting only certain trajectories of ensemble solution](https://discourse.julialang.org/t/plotting-only-certain-trajectories-of-ensemble-solution/93288)

<div class="topic-metadata">

**Author:** [@Torkel](https://discourse.julialang.org/u/Torkel)\
**Replies:** 1\
**Last updated:** [January 20, 2023, 11:56pm UTC](https://discourse.julialang.org/t/plotting-only-certain-trajectories-of-ensemble-solution/93288 "2023-01-20T23:56:48Z")

</div>

If I soleve an ensemble problem: using DifferentialEquations using Plots prob = ODEProblem((u,p,t)-\>1.01u,0.5,(0.0,1.0)) function prob\_func(prob,i,repeat) remake(prob,u0=rand()\*prob.u0) end ensemble\_prob = Ensemble…

---

## [Incorrect Infeasibility Status While Using MIPVerify](https://discourse.julialang.org/t/incorrect-infeasibility-status-while-using-mipverify/93033)

<div class="topic-metadata">

**Author:** [@Ahmed\_Anwar](https://discourse.julialang.org/u/Ahmed_Anwar)\
**Replies:** 2\
**Last updated:** [January 20, 2023, 12:26pm UTC](https://discourse.julialang.org/t/incorrect-infeasibility-status-while-using-mipverify/93033 "2023-01-20T12:26:12Z")

</div>

Hello Everyone, I am trying to use MIPVerify library for neural network verification against adversarial attacks. My goal is to replace the ReLU function by a trainable mask that depends on the input image. Therefore fo…

---

## [Fitting parameters together with initial conditions of an ODE to real data](https://discourse.julialang.org/t/fitting-parameters-together-with-initial-conditions-of-an-ode-to-real-data/93255)

<div class="topic-metadata">

**Author:** [@moeinkh88](https://discourse.julialang.org/u/moeinkh88)\
**Replies:** 1\
**Last updated:** [January 20, 2023, 11:10am UTC](https://discourse.julialang.org/t/fitting-parameters-together-with-initial-conditions-of-an-ode-to-real-data/93255 "2023-01-20T11:10:06Z")

</div>

I want to fit some parameters and some initial conditions of an ODE system to data. I am using Turing based on this tutorial. Example code: @model function fitprob(data) # Prior distributions. σ ~ InverseGamma…

---

## [Cuda makes Julia freeze](https://discourse.julialang.org/t/cuda-makes-julia-freeze/92040)

<div class="topic-metadata">

**Author:** [@J.J.sakurai](https://discourse.julialang.org/u/J.J.sakurai)\
**Replies:** 8\
**Last updated:** [January 20, 2023, 5:01am UTC](https://discourse.julialang.org/t/cuda-makes-julia-freeze/92040 "2023-01-20T05:01:33Z")

</div>

I want to use julia on GPU. So I installed CUDA by add CUDA and wrote the following code. julia\> using CUDA julia\> y = CUDA.fill(2.0, 3) I have been waiting more than an hour but still no response. It seems to freeze…

---

## [Using DataFrames \`combine\` is there a way to programmatically pass multiple functions to apply to the same same column?](https://discourse.julialang.org/t/using-dataframes-combine-is-there-a-way-to-programmatically-pass-multiple-functions-to-apply-to-the-same-same-column/93238)

<div class="topic-metadata">

**Author:** [@alex-s-gardner](https://discourse.julialang.org/u/alex-s-gardner)\
**Replies:** 9\
**Last updated:** [January 20, 2023, 3:52am UTC](https://discourse.julialang.org/t/using-dataframes-combine-is-there-a-way-to-programmatically-pass-multiple-functions-to-apply-to-the-same-same-column/93238 "2023-01-20T03:52:48Z")

</div>

Is there a way to pass a Vector{Function} like this: methods = \[mean, max\] to DataFrames combine? I was thinking something like this? df = combine(groupdf, nrow, 3 =\> methods) of df = combine(groupdf, nrow, 3 =\> met…

---

## [Almost but not quite univariate trend filtering / compressive sensing](https://discourse.julialang.org/t/almost-but-not-quite-univariate-trend-filtering-compressive-sensing/93237)

<div class="topic-metadata">

**Author:** [@sparrowhawk](https://discourse.julialang.org/u/sparrowhawk)\
**Replies:** 3\
**Last updated:** [January 20, 2023, 12:49am UTC](https://discourse.julialang.org/t/almost-but-not-quite-univariate-trend-filtering-compressive-sensing/93237 "2023-01-20T00:49:21Z")

</div>

I need to minimize for x a convex (I think it is convex from plots I have made, but I don’t understand L1 very well) subproblem of the form ||y - Ax||² + λ²|Dx|\_1 where D is a first difference operator for piecewise co…

---

## [New: FastParzenWindows.jl](https://discourse.julialang.org/t/new-fastparzenwindows-jl/93231)

<div class="topic-metadata">

**Author:** [@Nikos\_Gianniotis](https://discourse.julialang.org/u/Nikos_Gianniotis)\
**Replies:** 2\
**Last updated:** [January 19, 2023, 9:27pm UTC](https://discourse.julialang.org/t/new-fastparzenwindows-jl/93231 "2023-01-19T21:27:37Z")

</div>

Hello everyone, I have finally managed to sufficiently polish a package so that it can hopefully be useful to others. I have now registered it with the general registry. The package can be found at: The package imple…

---

## [Error in example: Parameter Estimation on Highly Stiff Systems](https://discourse.julialang.org/t/error-in-example-parameter-estimation-on-highly-stiff-systems/93229)

<div class="topic-metadata">

**Author:** [@kaido975](https://discourse.julialang.org/u/kaido975)\
**Replies:** 0\
**Last updated:** [January 19, 2023, 7:48pm UTC](https://discourse.julialang.org/t/error-in-example-parameter-estimation-on-highly-stiff-systems/93229 "2023-01-19T19:48:41Z")

</div>

The example provided at Parameter Estimation on Highly Stiff Systems · SciMLSensitivity.jl, does not seem to be working with Optimization.jl but works fine with DiffEqFlux.sciml\_train. I am getting the following error wi…

---

## [Flux + Automatic Differentiation](https://discourse.julialang.org/t/flux-automatic-differentiation/93224)

<div class="topic-metadata">

**Author:** [@parf](https://discourse.julialang.org/u/parf)\
**Replies:** 3\
**Last updated:** [January 19, 2023, 6:04pm UTC](https://discourse.julialang.org/t/flux-automatic-differentiation/93224 "2023-01-19T18:04:45Z")

</div>

Hi! I am trying to construct PINN to solve 1D Burgers equation in Flux.jl without using NeuralPDE.jl. The velocity field u(t,x) is defined by the neural net and I have to calculate the gradients with respect to t and x i…

---

## [Plotly: Add Slider to Vector of Heatmaps](https://discourse.julialang.org/t/plotly-add-slider-to-vector-of-heatmaps/93211)

<div class="topic-metadata">

**Author:** [@phiG](https://discourse.julialang.org/u/phiG)\
**Replies:** 2\
**Last updated:** [January 19, 2023, 5:53pm UTC](https://discourse.julialang.org/t/plotly-add-slider-to-vector-of-heatmaps/93211 "2023-01-19T17:53:55Z")

</div>

Hello everyone, this is my first post here so I hope I don’t ask something too trivial. I have a vector of matrices, each of same size, and I am able to plot a single one of them as a heatmap with PlotlyJS by accessing…

---

## [Finding periodic orbits with BifurcationKit.jl](https://discourse.julialang.org/t/finding-periodic-orbits-with-bifurcationkit-jl/93085)

<div class="topic-metadata">

**Author:** [@rusandris](https://discourse.julialang.org/u/rusandris)\
**Replies:** 5\
**Last updated:** [January 19, 2023, 5:35pm UTC](https://discourse.julialang.org/t/finding-periodic-orbits-with-bifurcationkit-jl/93085 "2023-01-19T17:35:02Z")

</div>

Hi everyone! Currently I’m working on a problem where I have non-standard Hopf-bifurcation points (Hopf-Hopf or degenerate Hopf) which get an “nd” (not documented) mark on my diagrams. Since these are not identified as …

---

## [Constraint declaration](https://discourse.julialang.org/t/constraint-declaration/92683)

<div class="topic-metadata">

**Author:** [@anonymus12](https://discourse.julialang.org/u/anonymus12)\
**Replies:** 8\
**Last updated:** [January 19, 2023, 4:51pm UTC](https://discourse.julialang.org/t/constraint-declaration/92683 "2023-01-19T16:51:33Z")

</div>

i want to build an MEXCLP using Julia/Jump in VS Code. Everything is working fine with the declaration except for this particular constraint. I tried this declaration: @constraint(model, c1, sum(x\[j\] for j in W\[i\]) \>=…

---

## [How to profile Zygote gradients](https://discourse.julialang.org/t/how-to-profile-zygote-gradients/90398)

<div class="topic-metadata">

**Author:** [@renatobellotti](https://discourse.julialang.org/u/renatobellotti)\
**Replies:** 2\
**Last updated:** [January 19, 2023, 3:43pm UTC](https://discourse.julialang.org/t/how-to-profile-zygote-gradients/90398 "2023-01-19T15:43:46Z")

</div>

Hi, I’m currently running an expensive optimisation with a code that uses Zygote for automatic differentiation. The gradient computation seems to dominate the computational cost. Is it possible to profile which parts of…

---

## [Converting entire custom structure from cpu to gpu and viceversa](https://discourse.julialang.org/t/converting-entire-custom-structure-from-cpu-to-gpu-and-viceversa/93203)

<div class="topic-metadata">

**Author:** [@carlogalli](https://discourse.julialang.org/u/carlogalli)\
**Replies:** 2\
**Last updated:** [January 19, 2023, 2:40pm UTC](https://discourse.julialang.org/t/converting-entire-custom-structure-from-cpu-to-gpu-and-viceversa/93203 "2023-01-19T14:40:56Z")

</div>

Hi! I am looking for a good way to convert a custom structure of Arrays to CuArrays and viceversa. I suspect the answer is related to them and to this topic, but I am not well versed in the usage of convert and adapt so…

---

## [Failure to enforce boundary conditions with NeuralPDE.jl (for a system with mixed partial derivatives)](https://discourse.julialang.org/t/failure-to-enforce-boundary-conditions-with-neuralpde-jl-for-a-system-with-mixed-partial-derivatives/93178)

<div class="topic-metadata">

**Author:** [@gianmariomanca](https://discourse.julialang.org/u/gianmariomanca)\
**Replies:** 0\
**Last updated:** [January 18, 2023, 6:37pm UTC](https://discourse.julialang.org/t/failure-to-enforce-boundary-conditions-with-neuralpde-jl-for-a-system-with-mixed-partial-derivatives/93178 "2023-01-18T18:37:19Z")

</div>

I’m trying to use NeuralPDE.jl to solve a simple system of PDEs and it seems to completely ignore the boundary conditions. Can someone suggest a modification of the following code that will shows the correct behavior? u…

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

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