# Specific Domains

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

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

**Page:** 530

---

## [Comparison of automatic differentiation tools from 2016 still accurate?](https://discourse.julialang.org/t/comparison-of-automatic-differentiation-tools-from-2016-still-accurate/8341)

<div class="topic-metadata">

**Author:** [@anriseth](https://discourse.julialang.org/u/anriseth)\
**Replies:** 41\
**Last updated:** [August 16, 2018, 6:18pm UTC](https://discourse.julialang.org/t/comparison-of-automatic-differentiation-tools-from-2016-still-accurate/8341 "2018-08-16T18:18:50Z")

</div>

I was made aware of a comprehensive comparison of AD tools at a talk I attended yesterday. The best reference I could find online was these slides: http://www.autodiff.org/ad16/Oral/Srajer\_Benchmark.pdf Does anyone know…

---

## [Flux for deep learning in particle physics](https://discourse.julialang.org/t/flux-for-deep-learning-in-particle-physics/13258)

<div class="topic-metadata">

**Author:** [@johnh](https://discourse.julialang.org/u/johnh)\
**Replies:** 6\
**Last updated:** [August 16, 2018, 4:25pm UTC](https://discourse.julialang.org/t/flux-for-deep-learning-in-particle-physics/13258 "2018-08-16T16:25:15Z")

</div>

There is a CodaLab code competition for applying deep learning to particle physics https://competitions.codalab.org/competitions/19818#learn\_the\_details As I have a background in the field (see me? See a CERN PhD whic…

---

## [User defined functions with vector constraints](https://discourse.julialang.org/t/user-defined-functions-with-vector-constraints/13574)

<div class="topic-metadata">

**Author:** [@Libbum](https://discourse.julialang.org/u/Libbum)\
**Replies:** 0\
**Last updated:** [August 16, 2018, 3:25pm UTC](https://discourse.julialang.org/t/user-defined-functions-with-vector-constraints/13574 "2018-08-16T15:25:37Z")

</div>

I’ve been thrown in the deep end with a problem that looks to be possible with current Julia tooling, but being new to the language and this particular subset of the mathematics I’m having a hard time. In a MWE format: …

---

## [Choosing order in QuadGK](https://discourse.julialang.org/t/choosing-order-in-quadgk/13534)

<div class="topic-metadata">

**Author:** [@sverek](https://discourse.julialang.org/u/sverek)\
**Replies:** 3\
**Last updated:** [August 16, 2018, 10:10am UTC](https://discourse.julialang.org/t/choosing-order-in-quadgk/13534 "2018-08-16T10:10:01Z")

</div>

Anybody has a rough idea how I should set the order of QuadGK for a standard mass matrix element calculation (Lagrange element) with order p^2 and BigFloat precision N i read this in documentation: It is advisable to i…

---

## [Invalid variable type present in SDP Optimal Power Flow formulation](https://discourse.julialang.org/t/invalid-variable-type-present-in-sdp-optimal-power-flow-formulation/13536)

<div class="topic-metadata">

**Author:** [@Alvaro\_Jose\_Gonzalez](https://discourse.julialang.org/u/Alvaro_Jose_Gonzalez)\
**Replies:** 3\
**Last updated:** [August 16, 2018, 2:21am UTC](https://discourse.julialang.org/t/invalid-variable-type-present-in-sdp-optimal-power-flow-formulation/13536 "2018-08-16T02:21:49Z")

</div>

I am trying to recreate the SDP formulation of the OPF given by Lavaei and Low 2013. After defining the whole model, I receive a warning of an invalid variable. Am I missing any type declaration besides the one for the …

---

## [How to index a \`CatagoricalArray\` then make a new array with the same levels](https://discourse.julialang.org/t/how-to-index-a-catagoricalarray-then-make-a-new-array-with-the-same-levels/12874)

<div class="topic-metadata">

**Author:** [@pdeffebach](https://discourse.julialang.org/u/pdeffebach)\
**Replies:** 4\
**Last updated:** [August 15, 2018, 5:51pm UTC](https://discourse.julialang.org/t/how-to-index-a-catagoricalarray-then-make-a-new-array-with-the-same-levels/12874 "2018-08-15T17:51:59Z")

</div>

I have been looking at the source code for CategoricalArray but I can’t figure this out. vec = CategoricalArray(\[1,2,1,2\]) t = vec\[1\] # We know what levels t.pool # CategoricalArrays.CategoricalPool{Int64,UInt32}(\[1,2\]…

---

## [Implementing parallel computing in finite element problems](https://discourse.julialang.org/t/implementing-parallel-computing-in-finite-element-problems/13473)

<div class="topic-metadata">

**Author:** [@thehalfspace](https://discourse.julialang.org/u/thehalfspace)\
**Replies:** 6\
**Last updated:** [August 15, 2018, 3:00pm UTC](https://discourse.julialang.org/t/implementing-parallel-computing-in-finite-element-problems/13473 "2018-08-15T15:00:13Z")

</div>

Hi all, I am very new to parallel programming, and relatively new to Julia. I have a spectral element code (higher order finite element) for numerical simulation of earthquakes, that I wish to parallelize. I have read …

---

## [JuMP/Ipopt failing on simple problem](https://discourse.julialang.org/t/jump-ipopt-failing-on-simple-problem/13209)

<div class="topic-metadata">

**Author:** [@Arrigo\_Benedetti](https://discourse.julialang.org/u/Arrigo_Benedetti)\
**Replies:** 4\
**Last updated:** [August 14, 2018, 10:44pm UTC](https://discourse.julialang.org/t/jump-ipopt-failing-on-simple-problem/13209 "2018-08-14T22:44:38Z")

</div>

In order to debug the issue that I mentioned in the thread https://discourse.julialang.org/t/error-in-jump-nlexpression/13147 I decided to write a simple program to solve a constrained nonlinear optimization problem in o…

---

## [Finding Roots with Roots.jl on ODE solutions](https://discourse.julialang.org/t/finding-roots-with-roots-jl-on-ode-solutions/13463)

<div class="topic-metadata">

**Author:** [@OvidiusCicero](https://discourse.julialang.org/u/OvidiusCicero)\
**Replies:** 6\
**Last updated:** [August 14, 2018, 10:08pm UTC](https://discourse.julialang.org/t/finding-roots-with-roots-jl-on-ode-solutions/13463 "2018-08-14T22:08:03Z")

</div>

I can find roots of an ODE solution nicely in this working, simplified example using DifferentialEquations using Plots using Roots prob = ODEProblem( function(u::Float64,p,t) return -u end, 1.0, (0.0,1.…

---

## [Share an already defined variable to all processes](https://discourse.julialang.org/t/share-an-already-defined-variable-to-all-processes/13421)

<div class="topic-metadata">

**Author:** [@e3c6](https://discourse.julialang.org/u/e3c6)\
**Replies:** 3\
**Last updated:** [August 14, 2018, 5:20pm UTC](https://discourse.julialang.org/t/share-an-already-defined-variable-to-all-processes/13421 "2018-08-14T17:20:55Z")

</div>

If I have a variable already defined in one process, how can I share it with all the other processes? @everywhere foo=1 does this for a new definition of foo. But suppose I already assigned foo to a value in the current…

---

## [Solving a system of differential equations with DifferentailEqations.jl](https://discourse.julialang.org/t/solving-a-system-of-differential-equations-with-differentaileqations-jl/13406)

<div class="topic-metadata">

**Author:** [@OvidiusCicero](https://discourse.julialang.org/u/OvidiusCicero)\
**Replies:** 2\
**Last updated:** [August 13, 2018, 9:27pm UTC](https://discourse.julialang.org/t/solving-a-system-of-differential-equations-with-differentaileqations-jl/13406 "2018-08-13T21:27:08Z")

</div>

Hello there when I try to solve a system of ordinary differential equations with the following code: using DifferentialEquations function dgl(du,u,p,t) du\[1\] = u\[3\] du\[2\] = u\[4\] du\[3\] = -sqrt(u\[4\]^2+u\[5\]^2)\*u\[3\] d…

---

## [Way to make radial gradient in 2-d PGFPlot?](https://discourse.julialang.org/t/way-to-make-radial-gradient-in-2-d-pgfplot/13271)

<div class="topic-metadata">

**Author:** [@djsegal](https://discourse.julialang.org/u/djsegal)\
**Replies:** 2\
**Last updated:** [August 13, 2018, 1:46pm UTC](https://discourse.julialang.org/t/way-to-make-radial-gradient-in-2-d-pgfplot/13271 "2018-08-13T13:46:53Z")

</div>

Let’s say you wanted to make a plot like: How would you make a radial gradient using Plots.jl with the PGFPlots backend?

---

## [System Dynamics packages](https://discourse.julialang.org/t/system-dynamics-packages/8222)

<div class="topic-metadata">

**Author:** [@smldis](https://discourse.julialang.org/u/smldis)\
**Replies:** 12\
**Last updated:** [August 12, 2018, 8:21pm UTC](https://discourse.julialang.org/t/system-dynamics-packages/8222 "2018-08-12T20:21:52Z")

</div>

Continuing the discussion from this post in 1.0 adoption path?: System dynamics ecosystem in julia needs coordination and development efforts. If anybody wants to discuss in this post I am kind of interested too. @jlper…

---

## [JuMP implementation of "Real-Time Order Fulfillment Decisions"](https://discourse.julialang.org/t/jump-implementation-of-real-time-order-fulfillment-decisions/13303)

<div class="topic-metadata">

**Author:** [@clausherther](https://discourse.julialang.org/u/clausherther)\
**Replies:** 0\
**Last updated:** [August 12, 2018, 3:28pm UTC](https://discourse.julialang.org/t/jump-implementation-of-real-time-order-fulfillment-decisions/13303 "2018-08-12T15:28:32Z")

</div>

Hello, I was wondering if anyone knew of a JuMP (or any other, really) implementation of this paper: “Benefits of Reevaluating Real-Time Order Fulfillment Decisions”. https://www.researchgate.net/publication/220391384…

---

## [Clean Numerical Simulation (Taylor Series Method)](https://discourse.julialang.org/t/clean-numerical-simulation-taylor-series-method/12893)

<div class="topic-metadata">

**Author:** [@m4r35n357](https://discourse.julialang.org/u/m4r35n357)\
**Replies:** 9\
**Last updated:** [August 12, 2018, 1:41pm UTC](https://discourse.julialang.org/t/clean-numerical-simulation-taylor-series-method/12893 "2018-08-12T13:41:30Z")

</div>

Hi all, very new to Julia, and not in academia BTW. I have been investigating the opinions and findings in this paper, and I thought I would contribute some Julia code which performs the aforementioned simulation to pot…

---

## [Observables for widgets that get rebuilt?](https://discourse.julialang.org/t/observables-for-widgets-that-get-rebuilt/13285)

<div class="topic-metadata">

**Author:** [@Raf](https://discourse.julialang.org/u/Raf)\
**Replies:** 0\
**Last updated:** [August 12, 2018, 10:12am UTC](https://discourse.julialang.org/t/observables-for-widgets-that-get-rebuilt/13285 "2018-08-12T10:12:01Z")

</div>

So I have a bunch of models that output the same plots, with parameter control. I want to be able to swap the models and the controller sliders update. Then using the new sliders should update the plots. I have it mostl…

---

## [Help with Makie at Hackfest?](https://discourse.julialang.org/t/help-with-makie-at-hackfest/13247)

<div class="topic-metadata">

**Author:** [@johnh](https://discourse.julialang.org/u/johnh)\
**Replies:** 1\
**Last updated:** [August 11, 2018, 7:08pm UTC](https://discourse.julialang.org/t/help-with-makie-at-hackfest/13247 "2018-08-11T19:08:56Z")

</div>

If anyone is good at using Makie could I have some hints and tips at the Hackfest?

---

## [2D Plots... some missing possibilities?](https://discourse.julialang.org/t/2d-plots-some-missing-possibilities/12894)

<div class="topic-metadata">

**Author:** [@BLI](https://discourse.julialang.org/u/BLI)\
**Replies:** 12\
**Last updated:** [August 10, 2018, 3:40pm UTC](https://discourse.julialang.org/t/2d-plots-some-missing-possibilities/12894 "2018-08-10T15:40:07Z")

</div>

I’m continuing to explore Plots (mainly 2D for now). Some useful things are missing, others are perhaps only hard to find? Figure title – has not been implemented yet. There are some workarounds, but these are somewhat…

---

## [Model solving - method error using NLsolve](https://discourse.julialang.org/t/model-solving-method-error-using-nlsolve/13035)

<div class="topic-metadata">

**Author:** [@axelle](https://discourse.julialang.org/u/axelle)\
**Replies:** 15\
**Last updated:** [August 10, 2018, 9:12am UTC](https://discourse.julialang.org/t/model-solving-method-error-using-nlsolve/13035 "2018-08-10T09:12:14Z")

</div>

Hi, I am stuck with a code (solving an economic model) and would need some help. I am sorry if the mistake is obvious, but I am new to this and would really need some help there. I have the following code: function bk\_o…

---

## [Error in JuMP @NLexpression](https://discourse.julialang.org/t/error-in-jump-nlexpression/13147)

<div class="topic-metadata">

**Author:** [@Arrigo\_Benedetti](https://discourse.julialang.org/u/Arrigo_Benedetti)\
**Replies:** 0\
**Last updated:** [August 10, 2018, 6:34am UTC](https://discourse.julialang.org/t/error-in-jump-nlexpression/13147 "2018-08-10T06:34:09Z")

</div>

I am running the following code with Julia 0.6.4: using JuMP using NLopt using FastGaussQuadrature using SpecialFunctions function FisherInf(args...) N = 32 t = \[x\[i\] for i=1:N\] lambda = \[args\[i\] for i…

---

## [Flux: GPU slower than CPU?](https://discourse.julialang.org/t/flux-gpu-slower-than-cpu/13072)

<div class="topic-metadata">

**Author:** [@rkat](https://discourse.julialang.org/u/rkat)\
**Replies:** 7\
**Last updated:** [August 10, 2018, 5:55am UTC](https://discourse.julialang.org/t/flux-gpu-slower-than-cpu/13072 "2018-08-10T05:55:08Z")

</div>

I have a loss function which I’ve defined using Flux. On my boring old CPU on one batch, it performs like this: julia\> @time loss(x, y) 0.071089 seconds (402 allocations: 20.235 MiB, 4.37% gc time) 1.3092503770133925 …

---

## [JLD and FileIO StackOverflowError: error when saving EPL fixtures data](https://discourse.julialang.org/t/jld-and-fileio-stackoverflowerror-error-when-saving-epl-fixtures-data/13139)

<div class="topic-metadata">

**Author:** [@xiaodai](https://discourse.julialang.org/u/xiaodai)\
**Replies:** 0\
**Last updated:** [August 10, 2018, 3:05am UTC](https://discourse.julialang.org/t/jld-and-fileio-stackoverflowerror-error-when-saving-epl-fixtures-data/13139 "2018-08-10T03:05:40Z")

</div>

I run the below code, which should work on your computer if you have all the packages installed. Basically it downloads some CSV from football-data.co.uk and tries to save it as a JLD file. The last line always gives an …

---

## [CLBlast, a tuned OpenCL BLAS library](https://discourse.julialang.org/t/clblast-a-tuned-opencl-blas-library/13082)

<div class="topic-metadata">

**Author:** [@ranocha](https://discourse.julialang.org/u/ranocha)\
**Replies:** 6\
**Last updated:** [August 9, 2018, 2:54pm UTC](https://discourse.julialang.org/t/clblast-a-tuned-opencl-blas-library/13082 "2018-08-09T14:54:22Z")

</div>

I’ve written a wrapper for CLBlast, a " tuned OpenCL BLAS library", which can be found at GitHub - ranocha/CLBlast.jl: Julia wrapper of CLBlast, a "tuned OpenCL BLAS library".. Most parts seem to work and there is a perf…

---

## [SatelliteToolbox.jl and ReferenceFrameRotations.jl has now full support for Julia 0.7 and 1.0](https://discourse.julialang.org/t/satellitetoolbox-jl-and-referenceframerotations-jl-has-now-full-support-for-julia-0-7-and-1-0/13054)

<div class="topic-metadata">

**Author:** [@Ronis\_BR](https://discourse.julialang.org/u/Ronis_BR)\
**Replies:** 0\
**Last updated:** [August 9, 2018, 2:38am UTC](https://discourse.julialang.org/t/satellitetoolbox-jl-and-referenceframerotations-jl-has-now-full-support-for-julia-0-7-and-1-0/13054 "2018-08-09T02:38:07Z")

</div>

Hi guys! I just tagged v0.3.0 of both SatelliteToolbox.jl and ReferenceFrameRotations.jl. With those new version, full support for v0.7 and v1.0 is added (no more warnings! :tada::tada::tada:). Unfortunately, together …

---

## [Deleting @NLconstraint from JuMP model](https://discourse.julialang.org/t/deleting-nlconstraint-from-jump-model/13002)

<div class="topic-metadata">

**Author:** [@chkwon](https://discourse.julialang.org/u/chkwon)\
**Replies:** 3\
**Last updated:** [August 8, 2018, 3:19pm UTC](https://discourse.julialang.org/t/deleting-nlconstraint-from-jump-model/13002 "2018-08-08T15:19:59Z")

</div>

JuMP v0.18: When I want to delete all @NLconstraints added to the JuMP model, is it safe and sufficient to do the following? m.nlpdata.nlconstr = Array{JuMP.GenericRangeConstraint{JuMP.NonlinearExprData},1}(0)

---

## [Memory allocation in multi-thread vs single-thread](https://discourse.julialang.org/t/memory-allocation-in-multi-thread-vs-single-thread/12999)

<div class="topic-metadata">

**Author:** [@aabrkay](https://discourse.julialang.org/u/aabrkay)\
**Replies:** 0\
**Last updated:** [August 7, 2018, 3:23pm UTC](https://discourse.julialang.org/t/memory-allocation-in-multi-thread-vs-single-thread/12999 "2018-08-07T15:23:38Z")

</div>

Dear All, I have a question about the significant difference in memory allocation between the multi-threads and single threads. I was experimenting Walk on Sphere method for Poisson equation, which is an embarrassingly p…

---

## [Gadfly Performance Issue](https://discourse.julialang.org/t/gadfly-performance-issue/12982)

<div class="topic-metadata">

**Author:** [@GantZA](https://discourse.julialang.org/u/GantZA)\
**Replies:** 2\
**Last updated:** [August 7, 2018, 11:39am UTC](https://discourse.julialang.org/t/gadfly-performance-issue/12982 "2018-08-07T11:39:55Z")

</div>

I’m trying to create a QQ plot for about 400 000 inter-arrival times, but the resulting plot takes far too long to open or julia crashes. Is there a way to plot more efficiently? plot\_test = plot( x = test1\[:ita\], …

---

## [ANN: MetaFields.jl and Flatten.jl - packages for writing composable models](https://discourse.julialang.org/t/ann-metafields-jl-and-flatten-jl-packages-for-writing-composable-models/12849)

<div class="topic-metadata">

**Author:** [@Raf](https://discourse.julialang.org/u/Raf)\
**Replies:** 10\
**Last updated:** [August 7, 2018, 10:14am UTC](https://discourse.julialang.org/t/ann-metafields-jl-and-flatten-jl-packages-for-writing-composable-models/12849 "2018-08-07T10:14:35Z")

</div>

I’ve been putting together a few biophysical models over the last year, and have constantly encountered what I feel to be the core software problem in my field: sharing high performance code in general ways never happens…

---

## [Has anyone made scatter plots with hexagonal binning?](https://discourse.julialang.org/t/has-anyone-made-scatter-plots-with-hexagonal-binning/12909)

<div class="topic-metadata">

**Author:** [@djsegal](https://discourse.julialang.org/u/djsegal)\
**Replies:** 7\
**Last updated:** [August 7, 2018, 8:01am UTC](https://discourse.julialang.org/t/has-anyone-made-scatter-plots-with-hexagonal-binning/12909 "2018-08-07T08:01:59Z")

</div>

Something came up in my research where I’m doing a scatter plot with way too many points. Has anyone had experience making hexagonally binned plots like the ones seen here: // preferably with pgfplots (or plots.jl), …

---

## [JuMP joins NumFOCUS](https://discourse.julialang.org/t/jump-joins-numfocus/12954)

<div class="topic-metadata">

**Author:** [@miles.lubin](https://discourse.julialang.org/u/miles.lubin)\
**Replies:** 0\
**Last updated:** [August 6, 2018, 6:53pm UTC](https://discourse.julialang.org/t/jump-joins-numfocus/12954 "2018-08-06T18:53:51Z")

</div>

Enjoy JuliaCon!

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

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