# Specific Domains

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

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

**Page:** 468

---

## [Barycentric interpolation formula with interval arithmetic](https://discourse.julialang.org/t/barycentric-interpolation-formula-with-interval-arithmetic/33380)

<div class="topic-metadata">

**Author:** [@ArjunNarayanan](https://discourse.julialang.org/u/ArjunNarayanan)\
**Replies:** 0\
**Last updated:** [January 15, 2020, 6:07am UTC](https://discourse.julialang.org/t/barycentric-interpolation-formula-with-interval-arithmetic/33380 "2020-01-15T06:07:16Z")

</div>

Barycentric interpolation is a numerically stable algorithm to evaluate Lagrange interpolating polynomials. An implementation in Julia is available at BarycentricInterpolation.jl. The procedure takes as input a set of n…

---

## [Make numbers from literature keep their reference](https://discourse.julialang.org/t/make-numbers-from-literature-keep-their-reference/33337)

<div class="topic-metadata">

**Author:** [@feanor12](https://discourse.julialang.org/u/feanor12)\
**Replies:** 0\
**Last updated:** [January 14, 2020, 8:06am UTC](https://discourse.julialang.org/t/make-numbers-from-literature-keep-their-reference/33337 "2020-01-14T08:06:38Z")

</div>

Julia is capable of defining new numbers, which can be used in conjunction with regular ones. To keep track of the source of the number(some material parameter) it should be possible to keep the DOI associated with it du…

---

## [Vega-Lite width and height in heatmap in Juno sometimes doesn't work!](https://discourse.julialang.org/t/vega-lite-width-and-height-in-heatmap-in-juno-sometimes-doesnt-work/33262)

<div class="topic-metadata">

**Author:** [@babaq](https://discourse.julialang.org/u/babaq)\
**Replies:** 8\
**Last updated:** [January 13, 2020, 10:57pm UTC](https://discourse.julialang.org/t/vega-lite-width-and-height-in-heatmap-in-juno-sometimes-doesnt-work/33262 "2020-01-13T22:57:15Z")

</div>

Hi, I tested the heatmap samples from the Vega-Lite.jl, such as, d=60 x = \[j for i in -d:d, j in -d:d\] y = \[i for i in -d:d, j in -d:d\] z = x.^2 .+ y.^2 data = DataFrame(x=vec(x'),y=vec(y'),z=vec(z')) data |\> @vlplot(…

---

## [Specifying integer edge weights in SimpleWeightedGraph](https://discourse.julialang.org/t/specifying-integer-edge-weights-in-simpleweightedgraph/33315)

<div class="topic-metadata">

**Author:** [@rkurchin](https://discourse.julialang.org/u/rkurchin)\
**Replies:** 1\
**Last updated:** [January 13, 2020, 9:01pm UTC](https://discourse.julialang.org/t/specifying-integer-edge-weights-in-simpleweightedgraph/33315 "2020-01-13T21:01:43Z")

</div>

Hi, I hope this is the right place to post this. It’s sort of a followup to this question but on a somewhat different aspect. In particular, I want to construct a SimpleWeightedGraph (which will eventually be passed to …

---

## [KernelError: recursion is currently not supported](https://discourse.julialang.org/t/kernelerror-recursion-is-currently-not-supported/33233)

<div class="topic-metadata">

**Author:** [@FujiwaraTakumiEH](https://discourse.julialang.org/u/FujiwaraTakumiEH)\
**Replies:** 4\
**Last updated:** [January 13, 2020, 10:40am UTC](https://discourse.julialang.org/t/kernelerror-recursion-is-currently-not-supported/33233 "2020-01-13T10:40:07Z")

</div>

Background: I changed my task to a simple example for description: I will take every element in array x, make it into a third-order matrix and calculate its determinant, and then assign the value of the determinant to…

---

## [Problem in @NLobjective](https://discourse.julialang.org/t/problem-in-nlobjective/33190)

<div class="topic-metadata">

**Author:** [@GiuliaM](https://discourse.julialang.org/u/GiuliaM)\
**Replies:** 2\
**Last updated:** [January 13, 2020, 10:22am UTC](https://discourse.julialang.org/t/problem-in-nlobjective/33190 "2020-01-13T10:22:00Z")

</div>

Hi everyone! I’m a new Julia user. I am using Julia 1.2 and JuMP v0.20.0 for some optimization problems. I have a constrained nonlinear optimization and I would like to define a nonlinear objective function composed b…

---

## [How to access a layer's parameters in Flux 0.10?](https://discourse.julialang.org/t/how-to-access-a-layers-parameters-in-flux-0-10/33201)

<div class="topic-metadata">

**Author:** [@homocomputeris](https://discourse.julialang.org/u/homocomputeris)\
**Replies:** 2\
**Last updated:** [January 13, 2020, 9:53am UTC](https://discourse.julialang.org/t/how-to-access-a-layers-parameters-in-flux-0-10/33201 "2020-01-13T09:53:20Z")

</div>

In Flux 0.10 par = params(model) returns a Zygote.Params object, which is an iterable but not indexed. How to get the parameters of the 2nd layer, for example?

---

## [Pmap but with control over which workers do the tasks?](https://discourse.julialang.org/t/pmap-but-with-control-over-which-workers-do-the-tasks/33270)

<div class="topic-metadata">

**Author:** [@Alan\_Crawford](https://discourse.julialang.org/u/Alan_Crawford)\
**Replies:** 0\
**Last updated:** [January 12, 2020, 5:51pm UTC](https://discourse.julialang.org/t/pmap-but-with-control-over-which-workers-do-the-tasks/33270 "2020-01-12T17:51:08Z")

</div>

I am wondering how to parallelize an expensive computation but where I specify which workers do the work? By way of background, in my use case each worker does a large computation on worker-specific objects (i.e. store…

---

## [Geometric DL library for Julia](https://discourse.julialang.org/t/geometric-dl-library-for-julia/22901)

<div class="topic-metadata">

**Author:** [@Azamat](https://discourse.julialang.org/u/Azamat)\
**Replies:** 5\
**Last updated:** [January 13, 2020, 6:25am UTC](https://discourse.julialang.org/t/geometric-dl-library-for-julia/22901 "2020-01-13T06:25:34Z")

</div>

I’m currently implementing geometric deep learning library. In particular, I’m extending Flux.jl to graphs and manifolds (taking advantage of the excellent LightGraphs.jl library). It’s working name is FluxGeometric.jl. …

---

## [Nonparametric estimation of spatial segregation in a multivariate point process](https://discourse.julialang.org/t/nonparametric-estimation-of-spatial-segregation-in-a-multivariate-point-process/33248)

<div class="topic-metadata">

**Author:** [@danriggins](https://discourse.julialang.org/u/danriggins)\
**Replies:** 7\
**Last updated:** [January 12, 2020, 4:32pm UTC](https://discourse.julialang.org/t/nonparametric-estimation-of-spatial-segregation-in-a-multivariate-point-process/33248 "2020-01-12T16:32:29Z")

</div>

Hi all, I’ve been trying to use the methods from this paper to estimate spatial segregation between crashes involving pedestrians versus those not in Indianapolis. The authors made this R package to facilitate the necess…

---

## [Flux 0.10.0 does not update function definitions in Jupyter](https://discourse.julialang.org/t/flux-0-10-0-does-not-update-function-definitions-in-jupyter/33261)

<div class="topic-metadata">

**Author:** [@RM314](https://discourse.julialang.org/u/RM314)\
**Replies:** 1\
**Last updated:** [January 12, 2020, 10:39am UTC](https://discourse.julialang.org/t/flux-0-10-0-does-not-update-function-definitions-in-jupyter/33261 "2020-01-12T10:39:42Z")

</div>

I stumbled onto this behaviour recently and I’m unsure if it’s intentional, or what to do about it. If I update a function definition, invocations of gradient() seem to use the old definition even as the rest of the note…

---

## [Please help with the very simplest 4-line Autodiff example: Flux0.9](https://discourse.julialang.org/t/please-help-with-the-very-simplest-4-line-autodiff-example-flux0-9/33257)

<div class="topic-metadata">

**Author:** [@jaynick](https://discourse.julialang.org/u/jaynick)\
**Replies:** 1\
**Last updated:** [January 12, 2020, 4:33am UTC](https://discourse.julialang.org/t/please-help-with-the-very-simplest-4-line-autodiff-example-flux0-9/33257 "2020-01-12T04:33:57Z")

</div>

This is with Flux 0.9, Julia1.3 Meta comment: Julia looks amazing. I keep coming back to it, but everything I try does not seem to work. Bad beginner luck. 1a) This example I want to work. It should find \[0,0\] as the…

---

## [Morphing algorithm in Julia?](https://discourse.julialang.org/t/morphing-algorithm-in-julia/33235)

<div class="topic-metadata">

**Author:** [@LaurentPlagne](https://discourse.julialang.org/u/LaurentPlagne)\
**Replies:** 1\
**Last updated:** [January 11, 2020, 12:03pm UTC](https://discourse.julialang.org/t/morphing-algorithm-in-julia/33235 "2020-01-11T12:03:37Z")

</div>

Hi, I wonder if it already exists a Julia tool to generate pretty transitions between different images. If not, I could try to translate existing implementations from other languages. Any hint or starting points ? La…

---

## [Need help debugging an NLopt optimization](https://discourse.julialang.org/t/need-help-debugging-an-nlopt-optimization/33229)

<div class="topic-metadata">

**Author:** [@mingjia1025](https://discourse.julialang.org/u/mingjia1025)\
**Replies:** 1\
**Last updated:** [January 11, 2020, 8:01am UTC](https://discourse.julialang.org/t/need-help-debugging-an-nlopt-optimization/33229 "2020-01-11T08:01:48Z")

</div>

I’m trying to use the LD\_MMA algoritmn in NLopt to fit my model. I get Xtol reached after 2 iterations and the values of the variables don’t seem to change from the starting values. Actually when I leave out the xtol re…

---

## [Best way to pass auxiliary arrays to Knet loss function](https://discourse.julialang.org/t/best-way-to-pass-auxiliary-arrays-to-knet-loss-function/33227)

<div class="topic-metadata">

**Author:** [@dfenn](https://discourse.julialang.org/u/dfenn)\
**Replies:** 0\
**Last updated:** [January 11, 2020, 1:58am UTC](https://discourse.julialang.org/t/best-way-to-pass-auxiliary-arrays-to-knet-loss-function/33227 "2020-01-11T01:58:10Z")

</div>

Julia beginner here, so I apologize if this is an obvious question. I’m implementing a model using Knet, and I’m writing a loss function that includes a penalty term that relies on information that is provided from outsi…

---

## [Custom adjoint for getindex in Zygote](https://discourse.julialang.org/t/custom-adjoint-for-getindex-in-zygote/33225)

<div class="topic-metadata">

**Author:** [@wgm](https://discourse.julialang.org/u/wgm)\
**Replies:** 0\
**Last updated:** [January 11, 2020, 1:37am UTC](https://discourse.julialang.org/t/custom-adjoint-for-getindex-in-zygote/33225 "2020-01-11T01:37:25Z")

</div>

I am new to Zygote and have a question about defining my own adjoint involving getindex with custom indexing. I’ve defined a type import Base: getindex struct ProbTensor{N} elements::Array{Float64,N} end function g…

---

## [How to scatter plot in Makie with custom marker shape](https://discourse.julialang.org/t/how-to-scatter-plot-in-makie-with-custom-marker-shape/33224)

<div class="topic-metadata">

**Author:** [@babaq](https://discourse.julialang.org/u/babaq)\
**Replies:** 0\
**Last updated:** [January 11, 2020, 1:10am UTC](https://discourse.julialang.org/t/how-to-scatter-plot-in-makie-with-custom-marker-shape/33224 "2020-01-11T01:10:22Z")

</div>

I see Plots support custom shape(defined by vertices) as marker, but how to do this in Makie.

---

## [Save objective values in an optimization loop](https://discourse.julialang.org/t/save-objective-values-in-an-optimization-loop/33192)

<div class="topic-metadata">

**Author:** [@yordiak](https://discourse.julialang.org/u/yordiak)\
**Replies:** 2\
**Last updated:** [January 10, 2020, 2:43pm UTC](https://discourse.julialang.org/t/save-objective-values-in-an-optimization-loop/33192 "2020-01-10T14:43:42Z")

</div>

Hello, here is my code using JuMP, Plots, Gurobi m=Model(with\_optimizer(Gurobi.Optimizer)) @variable(m, x\>=0) @variable(m, y\>=0) k=0 obj\_val=true while k\<=1 φ(x,y)=3\*x+k\*y @objective(m, Max, φ(x,y)) …

---

## [Maximum Likelihood Multivariate Model](https://discourse.julialang.org/t/maximum-likelihood-multivariate-model/33156)

<div class="topic-metadata">

**Author:** [@Maximilian](https://discourse.julialang.org/u/Maximilian)\
**Replies:** 15\
**Last updated:** [January 10, 2020, 12:09am UTC](https://discourse.julialang.org/t/maximum-likelihood-multivariate-model/33156 "2020-01-10T00:09:06Z")

</div>

I am trying to optimize a likelihood function for the multivariate case. This is the code I am using using Optim, LinearAlgebra, ForwardDiff #observed covariance cov\_obs = \[1.3629 0.408729 0.581832 0.…

---

## [Losing some output from redirect\_stdout() on remote workers](https://discourse.julialang.org/t/losing-some-output-from-redirect-stdout-on-remote-workers/33153)

<div class="topic-metadata">

**Author:** [@healyp](https://discourse.julialang.org/u/healyp)\
**Replies:** 2\
**Last updated:** [January 9, 2020, 4:52pm UTC](https://discourse.julialang.org/t/losing-some-output-from-redirect-stdout-on-remote-workers/33153 "2020-01-09T16:52:01Z")

</div>

I’m running some CPLEX (a Math. Programming solver) jobs as remote jobs and trapping each job’s output with the following (more or less): local res open(traceFile, "w") do out # https://discourse.julialang.org/t…

---

## [Stop Julia at certain optimiality gap and return decision variables](https://discourse.julialang.org/t/stop-julia-at-certain-optimiality-gap-and-return-decision-variables/2957)

<div class="topic-metadata">

**Author:** [@Duaa](https://discourse.julialang.org/u/Duaa)\
**Replies:** 12\
**Last updated:** [January 9, 2020, 4:43pm UTC](https://discourse.julialang.org/t/stop-julia-at-certain-optimiality-gap-and-return-decision-variables/2957 "2020-01-09T16:43:04Z")

</div>

I have a MILP problem and I am trying to solve it using CPLEX package through JuMP in Julia Language. However, the problem takes a long time to be solved. Is there a way to stop the program and return decision variables …

---

## [Connection between solver and Julia's JuMP Package](https://discourse.julialang.org/t/connection-between-solver-and-julias-jump-package/33076)

<div class="topic-metadata">

**Author:** [@Josa9321](https://discourse.julialang.org/u/Josa9321)\
**Replies:** 7\
**Last updated:** [January 8, 2020, 11:40pm UTC](https://discourse.julialang.org/t/connection-between-solver-and-julias-jump-package/33076 "2020-01-08T23:40:09Z")

</div>

I’m looking for ways or tips to connect a simple solver with the JuMP Package. I started reading the manual of Math Opt Interface (I don’t know the right terms, but I’ll call it the Backend of the JuMP), also I’m trying …

---

## [Animations - formats other than GIF](https://discourse.julialang.org/t/animations-formats-other-than-gif/32931)

<div class="topic-metadata">

**Author:** [@lwhitefox](https://discourse.julialang.org/u/lwhitefox)\
**Replies:** 7\
**Last updated:** [January 8, 2020, 9:57pm UTC](https://discourse.julialang.org/t/animations-formats-other-than-gif/32931 "2020-01-08T21:57:40Z")

</div>

I’m currently using Plots to generate graphics, including some animations. The docs make it easy to create animated GIFs but for my purposes I’d prefer an alternate output format such as .mov Any tips or suggestions on…

---

## [VegaLite.jl: size/resolution of pdf output](https://discourse.julialang.org/t/vegalite-jl-size-resolution-of-pdf-output/33021)

<div class="topic-metadata">

**Author:** [@BenjaminBorn](https://discourse.julialang.org/u/BenjaminBorn)\
**Replies:** 7\
**Last updated:** [January 8, 2020, 9:08pm UTC](https://discourse.julialang.org/t/vegalite-jl-size-resolution-of-pdf-output/33021 "2020-01-08T21:08:47Z")

</div>

Is there a way to specify with VegaLite (or maybe other packages) the resolution (or file size) of the pdf output? My goal is to use the pdf figures in a LaTeX document and because I have many figures and each is quite l…

---

## [WebIO not detected. (using Interact)](https://discourse.julialang.org/t/webio-not-detected-using-interact/33117)

<div class="topic-metadata">

**Author:** [@Simplex](https://discourse.julialang.org/u/Simplex)\
**Replies:** 0\
**Last updated:** [January 8, 2020, 8:04pm UTC](https://discourse.julialang.org/t/webio-not-detected-using-interact/33117 "2020-01-08T20:04:00Z")

</div>

I’ve an error when I use Interact. I use Anaconda. I’ve try a lot of things but never works. Thanks for your help ! using Interact @manipulate for i=1:10 @show i end

---

## [Compose.jl text error](https://discourse.julialang.org/t/compose-jl-text-error/33105)

<div class="topic-metadata">

**Author:** [@Simplex](https://discourse.julialang.org/u/Simplex)\
**Replies:** 4\
**Last updated:** [January 8, 2020, 7:25pm UTC](https://discourse.julialang.org/t/compose-jl-text-error/33105 "2020-01-08T19:25:58Z")

</div>

I try to use compose.jl to represent a network. But I also want to write text next to it but I have an error using text. I also tried to fix the network to make it not move at each iteration but it didn’t work. Thanks …

---

## [Automatic differentiation of f(x) constructed by decomposition of a matrix M(x)](https://discourse.julialang.org/t/automatic-differentiation-of-f-x-constructed-by-decomposition-of-a-matrix-m-x/33061)

<div class="topic-metadata">

**Author:** [@juandarias](https://discourse.julialang.org/u/juandarias)\
**Replies:** 7\
**Last updated:** [January 8, 2020, 5:36pm UTC](https://discourse.julialang.org/t/automatic-differentiation-of-f-x-constructed-by-decomposition-of-a-matrix-m-x/33061 "2020-01-08T17:36:43Z")

</div>

I have some objective function f(x), e.g. some type of error function for which I want to calculate its gradient and hessian. To construct f(x), I need the eigenvectors of a certain matrix \\textbf{M}(x). I have tried bo…

---

## [Gradient of \`sum(f, x)\` using Zygote, when \`f\` is real valued and \`x isa CuArray{\<:Complex}\`](https://discourse.julialang.org/t/gradient-of-sum-f-x-using-zygote-when-f-is-real-valued-and-x-isa-cuarray-complex/33081)

<div class="topic-metadata">

**Author:** [@Sleort](https://discourse.julialang.org/u/Sleort)\
**Replies:** 1\
**Last updated:** [January 8, 2020, 6:58am UTC](https://discourse.julialang.org/t/gradient-of-sum-f-x-using-zygote-when-f-is-real-valued-and-x-isa-cuarray-complex/33081 "2020-01-08T06:58:01Z")

</div>

While I’m able to differentiate sum(abs2, a) (using Zygote) when a isa CuArray{\<:Complex}, this fails for sum(f, a) with f = abs2 \[see below for details\]. I haven’t been able to figure out exactly why – and how to fix th…

---

## [Faster read only memory](https://discourse.julialang.org/t/faster-read-only-memory/32457)

<div class="topic-metadata">

**Author:** [@Michael\_Benton](https://discourse.julialang.org/u/Michael_Benton)\
**Replies:** 5\
**Last updated:** [January 8, 2020, 6:50am UTC](https://discourse.julialang.org/t/faster-read-only-memory/32457 "2020-01-08T06:50:57Z")

</div>

I’ve got something that works, but it needs to work faster. I’ve tested different methods outside of Julia that can do this much faster, and by comparison with peakflops.jl and a rudimentary calculation of how many opera…

---

## [Question on boundary conditions in DiffEqOperators.jl for 2D-case](https://discourse.julialang.org/t/question-on-boundary-conditions-in-diffeqoperators-jl-for-2d-case/33079)

<div class="topic-metadata">

**Author:** [@Eben60](https://discourse.julialang.org/u/Eben60)\
**Replies:** 0\
**Last updated:** [January 7, 2020, 8:46pm UTC](https://discourse.julialang.org/t/question-on-boundary-conditions-in-diffeqoperators-jl-for-2d-case/33079 "2020-01-07T20:46:47Z")

</div>

I just started to play around with DiffEqOperators. I’d like to extend the example at https://github.com/JuliaDiffEq/DiffEqOperators.jl/blob/master/examples/poisson.jl to a 2D case. Here is the adapted code: using Dif…

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

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