# Specific Domains

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

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

**Page:** 434

---

## [Handling Correlated shocks in HMM](https://discourse.julialang.org/t/handling-correlated-shocks-in-hmm/43355)

<div class="topic-metadata">

**Author:** [@balaji1975](https://discourse.julialang.org/u/balaji1975)\
**Replies:** 4\
**Last updated:** [July 21, 2020, 11:43am UTC](https://discourse.julialang.org/t/handling-correlated-shocks-in-hmm/43355 "2020-07-21T11:43:13Z")

</div>

Is there a way to handle the below in Turing? y(t) = f(y(t-1), s(t)) + g(y(t-1), s(t)) \* e1(t) s(t) = p(s(t-1)) + q(s(t-1)) \* e2(t) y(t) is the data (Tx1) and s(t) is the (hidden) state variable. e1 & e2 are correlate…

---

## [Convert a VariableRef in JuMP objective function or constraint](https://discourse.julialang.org/t/convert-a-variableref-in-jump-objective-function-or-constraint/43408)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 1\
**Last updated:** [July 21, 2020, 1:08am UTC](https://discourse.julialang.org/t/convert-a-variableref-in-jump-objective-function-or-constraint/43408 "2020-07-21T01:08:42Z")

</div>

I have a model with decision variables like this model = Model(Cbc.Optimizer) @variable(model, x\[1:num\_groups, 1:size(data, 1)\], Bin) and I need to be able to use the values of x to return a submatrix from another mat…

---

## [Extending the \`@formula\` syntax to generate functions](https://discourse.julialang.org/t/extending-the-formula-syntax-to-generate-functions/42882)

<div class="topic-metadata">

**Author:** [@BrendonChau](https://discourse.julialang.org/u/BrendonChau)\
**Replies:** 3\
**Last updated:** [July 20, 2020, 6:23pm UTC](https://discourse.julialang.org/t/extending-the-formula-syntax-to-generate-functions/42882 "2020-07-20T18:23:42Z")

</div>

I’m working on a project trying to use @formula in a fairly novel way, for example, if I wanted to use GLM.jl to fit a linear model with two predictors and their interaction I would use something like this: f = @formula…

---

## [Interested in using TaylorSeries.jl as an implementation of high-order automatic differentiation / AD , as presented in the book by W. Tucker](https://discourse.julialang.org/t/interested-in-using-taylorseries-jl-as-an-implementation-of-high-order-automatic-differentiation-ad-as-presented-in-the-book-by-w-tucker/43356)

<div class="topic-metadata">

**Author:** [@Marc.Cox](https://discourse.julialang.org/u/Marc.Cox)\
**Replies:** 5\
**Last updated:** [July 20, 2020, 4:08pm UTC](https://discourse.julialang.org/t/interested-in-using-taylorseries-jl-as-an-implementation-of-high-order-automatic-differentiation-ad-as-presented-in-the-book-by-w-tucker/43356 "2020-07-20T16:08:06Z")

</div>

Below is my (initial) solution to an issue I was having with TaylorSeries.jl , which leads to another question \>\> And more generally I’m interested in using TaylorSeries.jl as an implementation of high-order automatic di…

---

## [Help setting up a problem with JuMP and/or Clustering.jl](https://discourse.julialang.org/t/help-setting-up-a-problem-with-jump-and-or-clustering-jl/43252)

<div class="topic-metadata">

**Author:** [@mthelm85](https://discourse.julialang.org/u/mthelm85)\
**Replies:** 5\
**Last updated:** [July 19, 2020, 8:22pm UTC](https://discourse.julialang.org/t/help-setting-up-a-problem-with-jump-and-or-clustering-jl/43252 "2020-07-19T20:22:14Z")

</div>

Let’s say I want to divide all of the counties within a state into 4 groups and I want all of the counties in a group to be geographically contiguous. Furthermore, when I’m done, I want the total population of each of th…

---

## [JuMP Release Cadence](https://discourse.julialang.org/t/jump-release-cadence/43338)

<div class="topic-metadata">

**Author:** [@vtjeng](https://discourse.julialang.org/u/vtjeng)\
**Replies:** 2\
**Last updated:** [July 19, 2020, 5:54pm UTC](https://discourse.julialang.org/t/jump-release-cadence/43338 "2020-07-19T17:54:51Z")

</div>

Is there somewhere where I can track when the next release of JuMP is planned for? In this case, I’d like to know when https://github.com/jump-dev/JuMP.jl/pull/2275 will make it into a release (I’ve been holding off of …

---

## [A question about BilevelJuMP.jl](https://discourse.julialang.org/t/a-question-about-bileveljump-jl/43336)

<div class="topic-metadata">

**Author:** [@abcde](https://discourse.julialang.org/u/abcde)\
**Replies:** 1\
**Last updated:** [July 19, 2020, 5:49pm UTC](https://discourse.julialang.org/t/a-question-about-bileveljump-jl/43336 "2020-07-19T17:49:11Z")

</div>

using JuMP, BilevelJuMP, Xpress model = BilevelModel() @variable(Lower(model), x) @variable(Upper(model), y) @objective(Upper(model), Min, 3x + y) @constraints(Upper(model), begin x \<= 5 y \<= 8 y \>= 0 end)…

---

## [How to use dilate in Julia](https://discourse.julialang.org/t/how-to-use-dilate-in-julia/43241)

<div class="topic-metadata">

**Author:** [@sudo](https://discourse.julialang.org/u/sudo)\
**Replies:** 9\
**Last updated:** [July 19, 2020, 1:58pm UTC](https://discourse.julialang.org/t/how-to-use-dilate-in-julia/43241 "2020-07-19T13:58:05Z")

</div>

Hi! I am looking for something similar to python’s dilation function in skimage.morphology. I want to achieve precisely the following functionality: \>\>\> import numpy as np \>\>\> from skimage.morphology import dilation \>\>\> …

---

## [ForwardDiffSensitivity() and TrackerAdjoint() give different answers](https://discourse.julialang.org/t/forwarddiffsensitivity-and-trackeradjoint-give-different-answers/43254)

<div class="topic-metadata">

**Author:** [@fernando-duarte](https://discourse.julialang.org/u/fernando-duarte)\
**Replies:** 3\
**Last updated:** [July 18, 2020, 10:57pm UTC](https://discourse.julialang.org/t/forwarddiffsensitivity-and-trackeradjoint-give-different-answers/43254 "2020-07-18T22:57:16Z")

</div>

using DiffEqSensitivity, OrdinaryDiffEq, Zygote, StochasticDiffEq # similar to documentation, but out-of-place function fiip(u,p,t) du1 = dx = p\[1\]\*u\[1\] - p\[2\]\*u\[1\]\*u\[2\] du2 = dy = -p\[3\]\*u\[2\] + p\[4\]\*u\[1\]\*u\[2\] ret…

---

## [Error: JuMP constraint vs Array constraint (same size) --\> ArgumentError: Broadcasting](https://discourse.julialang.org/t/error-jump-constraint-vs-array-constraint-same-size-argumenterror-broadcasting/42812)

<div class="topic-metadata">

**Author:** [@PabloGalaz](https://discourse.julialang.org/u/PabloGalaz)\
**Replies:** 2\
**Last updated:** [July 18, 2020, 7:27pm UTC](https://discourse.julialang.org/t/error-jump-constraint-vs-array-constraint-same-size-argumenterror-broadcasting/42812 "2020-07-18T19:27:55Z")

</div>

Hello everyone, I’m using Julia v1.4.2 and JuMP v0.21.2 to solve a optimization problem. I want to compare a decision variable with an array of zeros and a calculated array (with other array and another decision variabl…

---

## [\[ANN\] CSV.jl 0.7 Release](https://discourse.julialang.org/t/ann-csv-jl-0-7-release/42162)

<div class="topic-metadata">

**Author:** [@quinnj](https://discourse.julialang.org/u/quinnj)\
**Replies:** 38\
**Last updated:** [July 18, 2020, 7:18pm UTC](https://discourse.julialang.org/t/ann-csv-jl-0-7-release/42162 "2020-07-18T19:18:17Z")

</div>

I’m happy to announce a new 0.7 release of the CSV.jl package. This is the last anticipated release before an official 1.0 release, and can be considered a 1.0 “trial run”. Major changes include: Deprecations CSV.read(…

---

## [Asking again about tick labels in Makie](https://discourse.julialang.org/t/asking-again-about-tick-labels-in-makie/43304)

<div class="topic-metadata">

**Author:** [@Noel\_Araujo](https://discourse.julialang.org/u/Noel_Araujo)\
**Replies:** 0\
**Last updated:** [July 18, 2020, 7:08pm UTC](https://discourse.julialang.org/t/asking-again-about-tick-labels-in-makie/43304 "2020-07-18T19:08:20Z")

</div>

Hellos :vulcan\_salute: I come upon this post of how to change tick labels with Makie, but is it not working because NT is not defined. Looking at Makie documentation, the solution should be the function: Function xtic…

---

## [Unfortunate interaction between ImageCore and Gadfly](https://discourse.julialang.org/t/unfortunate-interaction-between-imagecore-and-gadfly/43287)

<div class="topic-metadata">

**Author:** [@NicholasWMRitchie](https://discourse.julialang.org/u/NicholasWMRitchie)\
**Replies:** 1\
**Last updated:** [July 18, 2020, 3:39pm UTC](https://discourse.julialang.org/t/unfortunate-interaction-between-imagecore-and-gadfly/43287 "2020-07-18T15:39:36Z")

</div>

I’ve discovered an unfortunate interaction in which just loading ImageCore causes Gadfly to fail. This took me days to debug :frowning: as initially I assumed the error was in my code. It was only when I guessed that t…

---

## [JuAFEM mesh issue](https://discourse.julialang.org/t/juafem-mesh-issue/41467)

<div class="topic-metadata">

**Author:** [@bocc](https://discourse.julialang.org/u/bocc)\
**Replies:** 19\
**Last updated:** [July 18, 2020, 1:07pm UTC](https://discourse.julialang.org/t/juafem-mesh-issue/41467 "2020-07-18T13:07:36Z")

</div>

Hello, sorry in advance for posting this here, as this is most likely not a JuAFEM issue, but a problem with an external mesh. Also apologies for not being able to create a copy-pastable MWE. Any insight is appreciated …

---

## [FiniteDiff.finite\_difference\_gradient changes input when it errors?](https://discourse.julialang.org/t/finitediff-finite-difference-gradient-changes-input-when-it-errors/43279)

<div class="topic-metadata">

**Author:** [@amrods](https://discourse.julialang.org/u/amrods)\
**Replies:** 3\
**Last updated:** [July 18, 2020, 12:45pm UTC](https://discourse.julialang.org/t/finitediff-finite-difference-gradient-changes-input-when-it-errors/43279 "2020-07-18T12:45:01Z")

</div>

I ran into a surprising thing today. Consider the function sumlog: function sumlog(x) sum(log(i) for i in x) end I’m interested in taking the gradient of sumlog at x: using FiniteDiff x = \[1e-6, 1e-6\] FiniteDiff…

---

## [\[ANN\] Announcing Torch.jl](https://discourse.julialang.org/t/ann-announcing-torch-jl/42390)

<div class="topic-metadata">

**Author:** [@dhairyagandhi96](https://discourse.julialang.org/u/dhairyagandhi96)\
**Replies:** 14\
**Last updated:** [July 17, 2020, 6:59pm UTC](https://discourse.julialang.org/t/ann-announcing-torch-jl/42390 "2020-07-17T18:59:14Z")

</div>

I am excited to share some new development in the Julia ML stack. We are announcing Torch.jl which wraps the existing PyTorch and its CUDA kernels and makes them available through Julia. Torch.jl also exposes the kernels…

---

## [Projecting Histograms in 3D](https://discourse.julialang.org/t/projecting-histograms-in-3d/43149)

<div class="topic-metadata">

**Author:** [@Ian\_Slagle](https://discourse.julialang.org/u/Ian_Slagle)\
**Replies:** 6\
**Last updated:** [July 17, 2020, 5:42pm UTC](https://discourse.julialang.org/t/projecting-histograms-in-3d/43149 "2020-07-17T17:42:56Z")

</div>

Is it possible to make a plot like this in Julia? I’ve looked into plot3d, but I can’t get histograms to work with it.

---

## [Julia Vscode new version issue with JuMP](https://discourse.julialang.org/t/julia-vscode-new-version-issue-with-jump/43134)

<div class="topic-metadata">

**Author:** [@afazevedo](https://discourse.julialang.org/u/afazevedo)\
**Replies:** 4\
**Last updated:** [July 17, 2020, 10:35am UTC](https://discourse.julialang.org/t/julia-vscode-new-version-issue-with-jump/43134 "2020-07-17T10:35:05Z")

</div>

I’m having a problem with SparsesAxisArray JuMP on Julia for Vscode. When I use @variable with this condition below, it reproduces an error using Gurobi, JuMP model = Model(Gurobi.Optimizer) n = 5 V = 1:n V0= 1:n+…

---

## [A package/example for a simple neural network?](https://discourse.julialang.org/t/a-package-example-for-a-simple-neural-network/43007)

<div class="topic-metadata">

**Author:** [@zxjroger](https://discourse.julialang.org/u/zxjroger)\
**Replies:** 4\
**Last updated:** [July 17, 2020, 8:41am UTC](https://discourse.julialang.org/t/a-package-example-for-a-simple-neural-network/43007 "2020-07-17T08:41:06Z")

</div>

I try to find an easy-to-use package or example for training a forward-feed single-layer neural network. In R, the neuralnet and nnet packages are easy to use for a layman like me. But, I couldn’t find an equivalent pack…

---

## [ArchGDAL - memory issue with vector data](https://discourse.julialang.org/t/archgdal-memory-issue-with-vector-data/42827)

<div class="topic-metadata">

**Author:** [@julien\_goo](https://discourse.julialang.org/u/julien_goo)\
**Replies:** 3\
**Last updated:** [July 17, 2020, 6:20am UTC](https://discourse.julialang.org/t/archgdal-memory-issue-with-vector-data/42827 "2020-07-17T06:20:29Z")

</div>

Hello all, I’m starting to use geo data in Julia with ArcGDAL, after using the R sf package for a while. My use case is simple : I have WGS 84 positions and I want to guess the timezone they belong to. I have the world…

---

## [Help JuMP: we want to hear from people using JuMP to drive impactful changes to their business](https://discourse.julialang.org/t/help-jump-we-want-to-hear-from-people-using-jump-to-drive-impactful-changes-to-their-business/43191)

<div class="topic-metadata">

**Author:** [@odow](https://discourse.julialang.org/u/odow)\
**Replies:** 0\
**Last updated:** [July 16, 2020, 7:15pm UTC](https://discourse.julialang.org/t/help-jump-we-want-to-hear-from-people-using-jump-to-drive-impactful-changes-to-their-business/43191 "2020-07-16T19:15:28Z")

</div>

Dear all, JuMP is in the running for an award and we need your help. Have you used JuMP (either prototyping, or in deployment) to drive impactful changes to your business or organization? If so, would you be interest…

---

## [Help understanding MOI bridges, Even Proof-of-concept doesn't work](https://discourse.julialang.org/t/help-understanding-moi-bridges-even-proof-of-concept-doesnt-work/43152)

<div class="topic-metadata">

**Author:** [@votroto](https://discourse.julialang.org/u/votroto)\
**Replies:** 3\
**Last updated:** [July 16, 2020, 1:52pm UTC](https://discourse.julialang.org/t/help-understanding-moi-bridges-even-proof-of-concept-doesnt-work/43152 "2020-07-16T13:52:26Z")

</div>

Hi, I’ve been playing around with MOI bridges, but I can’t get even the simplest thing to work. Let’s imagine I want to create a bridge that simply ignores any PSD constraints: module MWE using MathOptInterface const M…

---

## [Locally solved](https://discourse.julialang.org/t/locally-solved/43179)

<div class="topic-metadata">

**Author:** [@mayar](https://discourse.julialang.org/u/mayar)\
**Replies:** 2\
**Last updated:** [July 16, 2020, 1:49pm UTC](https://discourse.julialang.org/t/locally-solved/43179 "2020-07-16T13:49:10Z")

</div>

hello there, I have used Ipopt with JuMP to solve an NLP problem and I got locally solved as a termination status…I tried changing the starting point but the results never changed What can I do to get an Optimal Soluti…

---

## [Plotting weighted graphs](https://discourse.julialang.org/t/plotting-weighted-graphs/43127)

<div class="topic-metadata">

**Author:** [@ohmsweetohm1](https://discourse.julialang.org/u/ohmsweetohm1)\
**Replies:** 3\
**Last updated:** [July 16, 2020, 12:14pm UTC](https://discourse.julialang.org/t/plotting-weighted-graphs/43127 "2020-07-16T12:14:45Z")

</div>

Is it possible to plot a SimpleWeightedGraphs.SimpleWeightedDiGraph? Using GraphPlot I get an error.

---

## [Need help with Implementing a simple planner in Julia in the context of the blocks world](https://discourse.julialang.org/t/need-help-with-implementing-a-simple-planner-in-julia-in-the-context-of-the-blocks-world/42928)

<div class="topic-metadata">

**Author:** [@Moreno](https://discourse.julialang.org/u/Moreno)\
**Replies:** 17\
**Last updated:** [July 16, 2020, 10:25am UTC](https://discourse.julialang.org/t/need-help-with-implementing-a-simple-planner-in-julia-in-the-context-of-the-blocks-world/42928 "2020-07-16T10:25:02Z")

</div>

Consider a collection of blocks and several tables. A block A can be in the box or on one of the tables at a given time. While on a table, a block A can be on the surface of the table or on top of another block B. The …

---

## [Drop incomplete groups from a DataFrame](https://discourse.julialang.org/t/drop-incomplete-groups-from-a-dataframe/43131)

<div class="topic-metadata">

**Author:** [@Tamas\_Papp](https://discourse.julialang.org/u/Tamas_Papp)\
**Replies:** 5\
**Last updated:** [July 16, 2020, 8:48am UTC](https://discourse.julialang.org/t/drop-incomplete-groups-from-a-dataframe/43131 "2020-07-16T08:48:16Z")

</div>

After removing some rows from a dataframe, I would like to keep only groups (grouping on a variable) which have “complete” observations (households always have 2 members in this data in the dataset I start with, so for …

---

## [DifferentialEquations.jl and systems with heterogenous units](https://discourse.julialang.org/t/differentialequations-jl-and-systems-with-heterogenous-units/2952)

<div class="topic-metadata">

**Author:** [@helgee](https://discourse.julialang.org/u/helgee)\
**Replies:** 17\
**Last updated:** [July 16, 2020, 1:31am UTC](https://discourse.julialang.org/t/differentialequations-jl-and-systems-with-heterogenous-units/2952 "2020-07-16T01:31:08Z")

</div>

Hi everybody, I am trying to solve the following system of ODE with DifferentialEquations.jl function newton!(t, y, dy, μ) r = norm(y\[1:3\]) dy\[1:3\] = y\[4:6\] dy\[4:6\] = -μ \* y\[1:3\] / r^3 end Without units th…

---

## [JuMP.all\_constraints not listing GenericQuadExpr constraints](https://discourse.julialang.org/t/jump-all-constraints-not-listing-genericquadexpr-constraints/43103)

<div class="topic-metadata">

**Author:** [@adow031](https://discourse.julialang.org/u/adow031)\
**Replies:** 8\
**Last updated:** [July 15, 2020, 11:16pm UTC](https://discourse.julialang.org/t/jump-all-constraints-not-listing-genericquadexpr-constraints/43103 "2020-07-15T23:16:35Z")

</div>

I have a quadratic program, and was trying to list all the constraints. However, all\_constraints, wouldn’t list any of the constraints of type GenericQuadExpr. Below is a small example exhibiting the problem. m = Model(…

---

## [Summary statistics for arbitrary variables within Turing model](https://discourse.julialang.org/t/summary-statistics-for-arbitrary-variables-within-turing-model/43120)

<div class="topic-metadata">

**Author:** [@balaji1975](https://discourse.julialang.org/u/balaji1975)\
**Replies:** 2\
**Last updated:** [July 15, 2020, 10:30pm UTC](https://discourse.julialang.org/t/summary-statistics-for-arbitrary-variables-within-turing-model/43120 "2020-07-15T22:30:38Z")

</div>

Hi - Is there a way to get summary stats on theta defined in the model below? @model schools(y, s) = begin N = length(y) eta = Array{Float64}(undef, N) theta = Array{Float64}(undef, N) mu ~ Uniform(-100, 100) tau…

---

## [Differentiable @unpack-ing of function arguments](https://discourse.julialang.org/t/differentiable-unpack-ing-of-function-arguments/43146)

<div class="topic-metadata">

**Author:** [@mefob](https://discourse.julialang.org/u/mefob)\
**Replies:** 9\
**Last updated:** [July 15, 2020, 10:27pm UTC](https://discourse.julialang.org/t/differentiable-unpack-ing-of-function-arguments/43146 "2020-07-15T22:27:47Z")

</div>

What is the recommended approach to @unpack parameters in a differential equation when using gradient wrt these parameters: I could not get named tuples or dicts to work, this is what I tried: using DiffEqSensitivi…

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

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