# Specific Domains

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

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

**Page:** 409

---

## [Solving vector field PDE with DiffEqOperators.jl](https://discourse.julialang.org/t/solving-vector-field-pde-with-diffeqoperators-jl/50391)

<div class="topic-metadata">

**Author:** [@jarias9](https://discourse.julialang.org/u/jarias9)\
**Replies:** 1\
**Last updated:** [November 18, 2020, 7:19pm UTC](https://discourse.julialang.org/t/solving-vector-field-pde-with-diffeqoperators-jl/50391 "2020-11-18T19:19:52Z")

</div>

Hello all, I am wondering if there is a workaround to solving a vector field for lets say Maxwell’s equations such as \\frac{\\partial \\mathbf{E}}{\\partial t} = \\frac{1}{\\varepsilon\_0} \\times \\mathbf{H} \\frac{\\partial…

---

## [Query status of the solution when using Roots](https://discourse.julialang.org/t/query-status-of-the-solution-when-using-roots/50328)

<div class="topic-metadata">

**Author:** [@optimist](https://discourse.julialang.org/u/optimist)\
**Replies:** 4\
**Last updated:** [November 18, 2020, 7:08pm UTC](https://discourse.julialang.org/t/query-status-of-the-solution-when-using-roots/50328 "2020-11-18T19:08:51Z")

</div>

Edited to make it more specific I would like to query status of the solution when using the Roots package, before the function returns an answer. I can see in the source code of “find\_zero.jl” that it has some functiona…

---

## [Not sure how to apply a filter](https://discourse.julialang.org/t/not-sure-how-to-apply-a-filter/50265)

<div class="topic-metadata">

**Author:** [@klaff](https://discourse.julialang.org/u/klaff)\
**Replies:** 7\
**Last updated:** [November 18, 2020, 4:11pm UTC](https://discourse.julialang.org/t/not-sure-how-to-apply-a-filter/50265 "2020-11-18T16:11:53Z")

</div>

I’ve constructed a filter as follows: filter1 = let f1 = 6.31 f2 = 1258.9 Q1 = 0.71 f3 = 15.915 f4 = 15.915 Q2 = 0.64 SecondOrderSections(\[ Biquad(4π^2\*f2^2,0.0,0.0,2π\*f1/Q1, 4π^2\*f1^2), Biquad(0.0,0.0,1.0,…

---

## [How can I set a Specific chain name for a Structure？(BioStrcuture Package)](https://discourse.julialang.org/t/how-can-i-set-a-specific-chain-name-for-a-structure-biostrcuture-package/50211)

<div class="topic-metadata">

**Author:** [@Wang-Lin-boop](https://discourse.julialang.org/u/Wang-Lin-boop)\
**Replies:** 6\
**Last updated:** [November 18, 2020, 3:26pm UTC](https://discourse.julialang.org/t/how-can-i-set-a-specific-chain-name-for-a-structure-biostrcuture-package/50211 "2020-11-18T15:26:13Z")

</div>

Recently, I was working on a project dealing with a large number of coordinate files, which required me to convert each chain in some coordinate files with MMCIF/PDBXML/MMTF/PDB format into PDB format. Some of the chains…

---

## [New Gurobi.jl syntax trouble](https://discourse.julialang.org/t/new-gurobi-jl-syntax-trouble/50344)

<div class="topic-metadata">

**Author:** [@afazevedo](https://discourse.julialang.org/u/afazevedo)\
**Replies:** 2\
**Last updated:** [November 18, 2020, 1:04pm UTC](https://discourse.julialang.org/t/new-gurobi-jl-syntax-trouble/50344 "2020-11-18T13:04:16Z")

</div>

Hi all I’ve got some issues in adapting to the new Gurobi.jl syntax. For example, this is a callback that I made in the old syntax: function cb(cb\_data, cb\_where) if cb\_where == CB\_MIPNODE #first node …

---

## [Sum-Of-Squares SOC constraint in JuMP](https://discourse.julialang.org/t/sum-of-squares-soc-constraint-in-jump/50091)

<div class="topic-metadata">

**Author:** [@andrewrosemberg](https://discourse.julialang.org/u/andrewrosemberg)\
**Replies:** 5\
**Last updated:** [November 16, 2020, 8:52am UTC](https://discourse.julialang.org/t/sum-of-squares-soc-constraint-in-jump/50091 "2020-11-16T08:52:04Z")

</div>

I’m trying to add a constraint of type (x'Ax \<= b) in JuMP to be interpreted as SOC by the solver. A is Positive definite. Seems to be the same error as presented in this issue, but A+A' has positive eigenvalues: eigv…

---

## [Gpu\_rand inside @cuda](https://discourse.julialang.org/t/gpu-rand-inside-cuda/50322)

<div class="topic-metadata">

**Author:** [@dgatey](https://discourse.julialang.org/u/dgatey)\
**Replies:** 2\
**Last updated:** [November 18, 2020, 10:41am UTC](https://discourse.julialang.org/t/gpu-rand-inside-cuda/50322 "2020-11-18T10:41:14Z")

</div>

Hey Everyone, I’ve managed to use gpu\_rand via gpu\_call but was wondering if it is possible to use gpu\_rand via @cuda as I’m currently using a 2D block. The @cuda function call only accepts bit types so I’m struggling …

---

## [Arpack and MKL](https://discourse.julialang.org/t/arpack-and-mkl/50291)

<div class="topic-metadata">

**Author:** [@Jonas](https://discourse.julialang.org/u/Jonas)\
**Replies:** 1\
**Last updated:** [November 18, 2020, 9:07am UTC](https://discourse.julialang.org/t/arpack-and-mkl/50291 "2020-11-18T09:07:20Z")

</div>

Hello everyone, I wrote some code to calculate the first n eigenvalues and eigenvectors of a generalized eigenvalue problem over a two-dimensional mesh. For this I used the function eigs() of the package Arpack. Since t…

---

## [Julia DifferentialEquations: increase maxiter](https://discourse.julialang.org/t/julia-differentialequations-increase-maxiter/50318)

<div class="topic-metadata">

**Author:** [@Luigi\_Marongiu](https://discourse.julialang.org/u/Luigi_Marongiu)\
**Replies:** 6\
**Last updated:** [November 18, 2020, 7:57am UTC](https://discourse.julialang.org/t/julia-differentialequations-increase-maxiter/50318 "2020-11-18T07:57:54Z")

</div>

Hello, I have set an ODEProblem object, but the procedure is stopped (although I still get t a solution) because a higher maxiters is needed. How can I check what is the default value and what is the syntax to increase…

---

## [Julia DifferentialEquations: add delay into ODE](https://discourse.julialang.org/t/julia-differentialequations-add-delay-into-ode/50337)

<div class="topic-metadata">

**Author:** [@Luigi\_Marongiu](https://discourse.julialang.org/u/Luigi_Marongiu)\
**Replies:** 3\
**Last updated:** [November 18, 2020, 7:03am UTC](https://discourse.julialang.org/t/julia-differentialequations-add-delay-into-ode/50337 "2020-11-18T07:03:01Z")

</div>

Hello, I have a set of ODE to model the growth of bacteria nad phages: function multiInfect!(du, u, p, t) μ, ν, κ, φ, ψ, ω, β, Β, ρ, τ, η, Η = p #= du\[1\] = resident susceptible du\[2\] = resident infected…

---

## [Some structures cannot be differentiated automatically?](https://discourse.julialang.org/t/some-structures-cannot-be-differentiated-automatically/50358)

<div class="topic-metadata">

**Author:** [@yeruoforever](https://discourse.julialang.org/u/yeruoforever)\
**Replies:** 0\
**Last updated:** [November 18, 2020, 5:46am UTC](https://discourse.julialang.org/t/some-structures-cannot-be-differentiated-automatically/50358 "2020-11-18T05:46:43Z")

</div>

Dear All, I am a student learning neural network. I implemented a small graph neural network model using Flux and MetaGraphs, but there was a problem that puzzled me. I posted the problematic code below, hope my friends…

---

## [\[ANN\] PopGen.jl v0.4.0 and PopGenSims.jl v0.0.2](https://discourse.julialang.org/t/ann-popgen-jl-v0-4-0-and-popgensims-jl-v0-0-2/50316)

<div class="topic-metadata">

**Author:** [@pdimens](https://discourse.julialang.org/u/pdimens)\
**Replies:** 4\
**Last updated:** [November 17, 2020, 9:47pm UTC](https://discourse.julialang.org/t/ann-popgen-jl-v0-4-0-and-popgensims-jl-v0-0-2/50316 "2020-11-17T21:47:05Z")

</div>

Hello all, I’m thrilled to announce the next release of PopGen.jl v0.4.0 (although incorrectly versioned in the release, oops). Internal tweaks removing CategoricalArrays in favor of PooledArrays cleaning up and adding…

---

## [JuMP model Polyhedron](https://discourse.julialang.org/t/jump-model-polyhedron/50108)

<div class="topic-metadata">

**Author:** [@JohnZ](https://discourse.julialang.org/u/JohnZ)\
**Replies:** 11\
**Last updated:** [November 18, 2020, 1:33am UTC](https://discourse.julialang.org/t/jump-model-polyhedron/50108 "2020-11-18T01:33:34Z")

</div>

Hi, I am trying to plot the feasible set of a JuMP model as a 3D figure. I have used the syntax provided in Polyhedra.jl documentation to create a polyhedron: Optimization · Polyhedra A simple example is given below: …

---

## [How to get better performance from NUTS in Turing?](https://discourse.julialang.org/t/how-to-get-better-performance-from-nuts-in-turing/50349)

<div class="topic-metadata">

**Author:** [@mkarikom](https://discourse.julialang.org/u/mkarikom)\
**Replies:** 0\
**Last updated:** [November 18, 2020, 12:54am UTC](https://discourse.julialang.org/t/how-to-get-better-performance-from-nuts-in-turing/50349 "2020-11-18T00:54:32Z")

</div>

I’m wondering if anyone has further performance tips to speed this up? This MWE is an attempt to use NUTS on data that is similar in size to my actual use case, following the model introduced in this post, but no longer…

---

## [DBnomics - World Economic Database](https://discourse.julialang.org/t/dbnomics-world-economic-database/50242)

<div class="topic-metadata">

**Author:** [@tlorans](https://discourse.julialang.org/u/tlorans)\
**Replies:** 1\
**Last updated:** [November 18, 2020, 12:28am UTC](https://discourse.julialang.org/t/dbnomics-world-economic-database/50242 "2020-11-18T00:28:32Z")

</div>

Hello everybody, This is my first post in this forum. :slight\_smile: I would like to point out the existence of a Julia version of the DBnomics API. The DBnomics project is led by French academics and aims to concatena…

---

## [Substituting state and parameter values in jacobian of DifferentialEquations.jl system using ModelingToolkit](https://discourse.julialang.org/t/substituting-state-and-parameter-values-in-jacobian-of-differentialequations-jl-system-using-modelingtoolkit/50264)

<div class="topic-metadata">

**Author:** [@claudio20497](https://discourse.julialang.org/u/claudio20497)\
**Replies:** 2\
**Last updated:** [November 17, 2020, 11:48pm UTC](https://discourse.julialang.org/t/substituting-state-and-parameter-values-in-jacobian-of-differentialequations-jl-system-using-modelingtoolkit/50264 "2020-11-17T23:48:31Z")

</div>

So I would like to explicitly evaluate the R\_0 of an epidemiological model written in DifferentialEquations.jl using ModelingToolkit.jl. Here is a simple SIRD model: using DifferentialEquations, ModelingToolkit # Mode…

---

## [Float64 error when using ForwardDiff.gradient for nested functions](https://discourse.julialang.org/t/float64-error-when-using-forwarddiff-gradient-for-nested-functions/50343)

<div class="topic-metadata">

**Author:** [@yips](https://discourse.julialang.org/u/yips)\
**Replies:** 0\
**Last updated:** [November 17, 2020, 11:22pm UTC](https://discourse.julialang.org/t/float64-error-when-using-forwarddiff-gradient-for-nested-functions/50343 "2020-11-17T23:22:54Z")

</div>

I am trying to use ForwardDiff.gradient find the gradient of the following function: function loss\_test(x::Array{T,2},SecreteFuncs,pInf::Array{T,2},cmat::Array{T,2}) where {T \<: Real} cmat\_sol = Findcss(x,SecreteFuncs…

---

## [Optimization of Curves](https://discourse.julialang.org/t/optimization-of-curves/50312)

<div class="topic-metadata">

**Author:** [@brett\_knoss](https://discourse.julialang.org/u/brett_knoss)\
**Replies:** 8\
**Last updated:** [November 17, 2020, 8:50pm UTC](https://discourse.julialang.org/t/optimization-of-curves/50312 "2020-11-17T20:50:55Z")

</div>

What is a good option for finding the x value of a tangent point of a curve (such as matching price and value), or finding the minimum or maximum of a curve.

---

## [Scatter3d and pyplot backend in Plots.jl](https://discourse.julialang.org/t/scatter3d-and-pyplot-backend-in-plots-jl/50238)

<div class="topic-metadata">

**Author:** [@keithjlee](https://discourse.julialang.org/u/keithjlee)\
**Replies:** 4\
**Last updated:** [November 17, 2020, 3:43pm UTC](https://discourse.julialang.org/t/scatter3d-and-pyplot-backend-in-plots-jl/50238 "2020-11-17T15:43:46Z")

</div>

Hello! I’m having a very strange issue trying to visualize a 3D scatterplot using the pyplot backend in Plots.jl. The Z-axis position of my inputted coordinates do not seem to register, and only stays on the 0 value, whi…

---

## [DifferentialEquations - Derivatives in ODE function/ nesting AD](https://discourse.julialang.org/t/differentialequations-derivatives-in-ode-function-nesting-ad/50100)

<div class="topic-metadata">

**Author:** [@Andreas\_Schlaginhauf](https://discourse.julialang.org/u/Andreas_Schlaginhauf)\
**Replies:** 6\
**Last updated:** [November 17, 2020, 2:57pm UTC](https://discourse.julialang.org/t/differentialequations-derivatives-in-ode-function-nesting-ad/50100 "2020-11-17T14:57:21Z")

</div>

If the ODE function of a neural ODE includes the gradient f’(x) of a neural network f(x), we need to calculate second order gradients of f(x) in the backward pass. For me this led to all sorts of problems when the AD met…

---

## [Colorbar cutoff in Plots.jl with pyplot](https://discourse.julialang.org/t/colorbar-cutoff-in-plots-jl-with-pyplot/48178)

<div class="topic-metadata">

**Author:** [@gideonsimpson](https://discourse.julialang.org/u/gideonsimpson)\
**Replies:** 2\
**Last updated:** [November 17, 2020, 12:47pm UTC](https://discourse.julialang.org/t/colorbar-cutoff-in-plots-jl-with-pyplot/48178 "2020-11-17T12:47:06Z")

</div>

I’ve got the following issue (see attached figure), where, with my font selections, the colorbar tick labels are cutoff. How do I expand the figure to include them or shrink the width of the colorbar?

---

## [\`norm\` function error on CuDeviceArray](https://discourse.julialang.org/t/norm-function-error-on-cudevicearray/50273)

<div class="topic-metadata">

**Author:** [@HaoxuanGuo](https://discourse.julialang.org/u/HaoxuanGuo)\
**Replies:** 1\
**Last updated:** [November 17, 2020, 7:12am UTC](https://discourse.julialang.org/t/norm-function-error-on-cudevicearray/50273 "2020-11-17T07:12:21Z")

</div>

This code will get an error about norm function using CUDA using Adapt using LinearAlgebra struct TestType{T\<:AbstractArray} a::T b::T c::Float32 end Adapt.adapt\_structure(to, x::TestType) = TestType(adapt…

---

## [In-place orthogonal/unitary transforms?](https://discourse.julialang.org/t/in-place-orthogonal-unitary-transforms/50111)

<div class="topic-metadata">

**Author:** [@pkairys](https://discourse.julialang.org/u/pkairys)\
**Replies:** 12\
**Last updated:** [November 17, 2020, 10:46am UTC](https://discourse.julialang.org/t/in-place-orthogonal-unitary-transforms/50111 "2020-11-17T10:46:43Z")

</div>

Say I want to perform a unitary/orthogonal transformation of an operator: U^\\dagger A U Is there a way to do this without allocating a new array? I see that there is a method for lmul!(A,B) that will do AB and place t…

---

## [Random walks inside a model](https://discourse.julialang.org/t/random-walks-inside-a-model/49994)

<div class="topic-metadata">

**Author:** [@Camilo1704](https://discourse.julialang.org/u/Camilo1704)\
**Replies:** 4\
**Last updated:** [November 17, 2020, 8:49am UTC](https://discourse.julialang.org/t/random-walks-inside-a-model/49994 "2020-11-17T08:49:02Z")

</div>

Hi guys, I’m wondering if there is a way to include a variable which is a random walk in a model. Particularly, if inside my model I have the lines of code shown below, is there a walk around to do it more efficiently? …

---

## [Are there any way to copy a Union Array to GPU with CUDA.jl?](https://discourse.julialang.org/t/are-there-any-way-to-copy-a-union-array-to-gpu-with-cuda-jl/50272)

<div class="topic-metadata">

**Author:** [@HaoxuanGuo](https://discourse.julialang.org/u/HaoxuanGuo)\
**Replies:** 1\
**Last updated:** [November 17, 2020, 7:13am UTC](https://discourse.julialang.org/t/are-there-any-way-to-copy-a-union-array-to-gpu-with-cuda-jl/50272 "2020-11-17T07:13:41Z")

</div>

Look at this example abstract type AbstractTest end struct ATest \<: AbstractTest end struct BTest \<: AbstractTest end struct CTest{T \<: AbstractTest} a::T end using CUDA CuArray(\[CTest(ATest()), CTest(BTest())\]) …

---

## [Are there any way to copy a Dict or a custom datatype which contains an Array to GPU by CUDA?](https://discourse.julialang.org/t/are-there-any-way-to-copy-a-dict-or-a-custom-datatype-which-contains-an-array-to-gpu-by-cuda/50217)

<div class="topic-metadata">

**Author:** [@HaoxuanGuo](https://discourse.julialang.org/u/HaoxuanGuo)\
**Replies:** 6\
**Last updated:** [November 17, 2020, 7:13am UTC](https://discourse.julialang.org/t/are-there-any-way-to-copy-a-dict-or-a-custom-datatype-which-contains-an-array-to-gpu-by-cuda/50217 "2020-11-17T07:13:01Z")

</div>

I want to use the LinearAlgebra library to simplify my programming. So, I warp some Arrays in an struct. But I find I can not copy it to GPU because it is a non-bitstype. Similarly, I want to look up some tables which i…

---

## [Really small plot picture](https://discourse.julialang.org/t/really-small-plot-picture/18247)

<div class="topic-metadata">

**Author:** [@muriloasouza](https://discourse.julialang.org/u/muriloasouza)\
**Replies:** 8\
**Last updated:** [November 17, 2020, 5:48am UTC](https://discourse.julialang.org/t/really-small-plot-picture/18247 "2020-11-17T05:48:42Z")

</div>

I’m trying to plot a groupedbar graphic using Statsplots, but i’m getting a really small image, like this: Is there anyway i can plot another groupedbar graphic in this same image? Like above or below that one?

---

## [Request for comments: filenames with metadata](https://discourse.julialang.org/t/request-for-comments-filenames-with-metadata/50266)

<div class="topic-metadata">

**Author:** [@alhirzel](https://discourse.julialang.org/u/alhirzel)\
**Replies:** 1\
**Last updated:** [November 17, 2020, 3:05am UTC](https://discourse.julialang.org/t/request-for-comments-filenames-with-metadata/50266 "2020-11-17T03:05:12Z")

</div>

When managing some bills using a toy framework I’m developing in Julia, I am finding it’s useful to embed some metadata in the filename. I wanted to share an example here for comments: does anyone know of a library in an…

---

## [Makie 3D mesh plot is not as good as before](https://discourse.julialang.org/t/makie-3d-mesh-plot-is-not-as-good-as-before/39808)

<div class="topic-metadata">

**Author:** [@yusri-dh](https://discourse.julialang.org/u/yusri-dh)\
**Replies:** 6\
**Last updated:** [November 16, 2020, 11:16pm UTC](https://discourse.julialang.org/t/makie-3d-mesh-plot-is-not-as-good-as-before/39808 "2020-11-16T23:16:24Z")

</div>

In the Makie v0.10.0, plotting 3D mesh is not as good as before. Here is minimal working example using Makie v0.10.0, shell\> wget https://github.com/libigl/libigl-tutorial-data/raw/master/bumpy-cube.obj using Makie, Fi…

---

## [Sankey diagrams](https://discourse.julialang.org/t/sankey-diagrams/8668)

<div class="topic-metadata">

**Author:** [@Alvaro\_Jose\_Gonzalez](https://discourse.julialang.org/u/Alvaro_Jose_Gonzalez)\
**Replies:** 3\
**Last updated:** [November 16, 2020, 9:48pm UTC](https://discourse.julialang.org/t/sankey-diagrams/8668 "2020-11-16T21:48:30Z")

</div>

Hello everybody, Does anyone know how to create a Sankey diagram in Julia? Thanks! Alvaro

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

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