# Specific Domains

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

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

**Page:** 288

---

## [Easy plotting of vector valued function in Plots.jl](https://discourse.julialang.org/t/easy-plotting-of-vector-valued-function-in-plots-jl/79788)

<div class="topic-metadata">

**Author:** [@Joseph\_Bradley](https://discourse.julialang.org/u/Joseph_Bradley)\
**Replies:** 4\
**Last updated:** [April 21, 2022, 4:59pm UTC](https://discourse.julialang.org/t/easy-plotting-of-vector-valued-function-in-plots-jl/79788 "2022-04-21T16:59:28Z")

</div>

Hi, I’d like to plot a function from R -\> R^N in Plots.jl. My code is: using Plots N = 9 v₁ = ones(N) v₂ = \[1, -8, ones(N-2)...\] z₀ = 1.0 function z(t; v₁ = v₁, v₂ = v₂, z₀ = z₀) (z₀ / 9.0) .\* (v₁ - exp(-9.0\*t)\*v₂) …

---

## [Adding labels to surfaces in the same figure](https://discourse.julialang.org/t/adding-labels-to-surfaces-in-the-same-figure/79708)

<div class="topic-metadata">

**Author:** [@gnicolosi](https://discourse.julialang.org/u/gnicolosi)\
**Replies:** 4\
**Last updated:** [April 21, 2022, 3:53pm UTC](https://discourse.julialang.org/t/adding-labels-to-surfaces-in-the-same-figure/79708 "2022-04-21T15:53:43Z")

</div>

Hi all, how do I add labels to the surfaces below? Using label= as when plotting curves in 2d does not work. There’s very few information on attributes to surface() and I can’t find a solution anywhere else for this simp…

---

## [Group Selection by Condition in DataFramesMeta.jl](https://discourse.julialang.org/t/group-selection-by-condition-in-dataframesmeta-jl/79784)

<div class="topic-metadata">

**Author:** [@Rahul](https://discourse.julialang.org/u/Rahul)\
**Replies:** 1\
**Last updated:** [April 21, 2022, 3:01pm UTC](https://discourse.julialang.org/t/group-selection-by-condition-in-dataframesmeta-jl/79784 "2022-04-21T15:01:46Z")

</div>

I am looking to select groups that follow a condition using DataFramesMeta.jl. Reference example in R that I am looking to reproduce: require(data.table) df = data.table(id = c(rep(1, 10), rep(2, 10), rep(3, 10)), …

---

## [Unable to transform a quadratic constraint into a second-order cone constraint because the quadratic constraint is not strongly convex](https://discourse.julialang.org/t/unable-to-transform-a-quadratic-constraint-into-a-second-order-cone-constraint-because-the-quadratic-constraint-is-not-strongly-convex/79767)

<div class="topic-metadata">

**Author:** [@Eric.O](https://discourse.julialang.org/u/Eric.O)\
**Replies:** 4\
**Last updated:** [April 21, 2022, 1:13pm UTC](https://discourse.julialang.org/t/unable-to-transform-a-quadratic-constraint-into-a-second-order-cone-constraint-because-the-quadratic-constraint-is-not-strongly-convex/79767 "2022-04-21T13:13:38Z")

</div>

I am using Mosek solver to solve the model below. I get the model showing but the solver is not able to display the results. Am guessing is an issue with the second order cone formulation. I will appreciate your guidance…

---

## [ODE Integrator Function](https://discourse.julialang.org/t/ode-integrator-function/79756)

<div class="topic-metadata">

**Author:** [@BinaryBot](https://discourse.julialang.org/u/BinaryBot)\
**Replies:** 6\
**Last updated:** [April 21, 2022, 10:53am UTC](https://discourse.julialang.org/t/ode-integrator-function/79756 "2022-04-21T10:53:05Z")

</div>

My ODE runs from tspan(1,120). It needs an integrator function that can do the following: Store the solutions at timestep 2,3,4…120 (i.e integral numbers) in a list. Although step size does not have to be 1, it just h…

---

## [Why does setting JULIA\_EXLCUSIVE=1 make MKL run single-threaded?](https://discourse.julialang.org/t/why-does-setting-julia-exlcusive-1-make-mkl-run-single-threaded/79766)

<div class="topic-metadata">

**Author:** [@jishnub](https://discourse.julialang.org/u/jishnub)\
**Replies:** 0\
**Last updated:** [April 21, 2022, 7:22am UTC](https://discourse.julialang.org/t/why-does-setting-julia-exlcusive-1-make-mkl-run-single-threaded/79766 "2022-04-21T07:22:49Z")

</div>

I am working on a cluster with 28 cores allocated to me. I start Julia using MKL\_DYNAMIC="FALSE" JULIA\_EXCLUSIVE=1 julia although I’m uncertain if this is the right approach. Eventually, I aim to run multiple julia th…

---

## [Writing a QP to MPS?](https://discourse.julialang.org/t/writing-a-qp-to-mps/79764)

<div class="topic-metadata">

**Author:** [@gleyland](https://discourse.julialang.org/u/gleyland)\
**Replies:** 1\
**Last updated:** [April 21, 2022, 5:40am UTC](https://discourse.julialang.org/t/writing-a-qp-to-mps/79764 "2022-04-21T05:40:15Z")

</div>

Hi, As far as I can tell, you can’t write a QP to an MPS file with write\_to\_file, but I’ve told that MPS does support QPs. Are either of those correct? (ie, should I be able to write QP to MPS, but I’m doing it wrong, …

---

## [Problem simulating Matlab ODE15s in Julia](https://discourse.julialang.org/t/problem-simulating-matlab-ode15s-in-julia/79687)

<div class="topic-metadata">

**Author:** [@BinaryBot](https://discourse.julialang.org/u/BinaryBot)\
**Replies:** 8\
**Last updated:** [April 20, 2022, 5:28pm UTC](https://discourse.julialang.org/t/problem-simulating-matlab-ode15s-in-julia/79687 "2022-04-20T17:28:33Z")

</div>

Here is the code I have, using Pkg using PyCall using CSV; using DataFrames; using Dierckx;using DifferentialEquations;using LSODA;using MATLABDiffEq; using ParameterizedFunctions;using IfElse; p = \[ 1.35e-2 #k1 …

---

## [Reading ENVI style multispectral images with archGDAL](https://discourse.julialang.org/t/reading-envi-style-multispectral-images-with-archgdal/79702)

<div class="topic-metadata">

**Author:** [@Arrigo\_Benedetti](https://discourse.julialang.org/u/Arrigo_Benedetti)\
**Replies:** 4\
**Last updated:** [April 20, 2022, 2:09pm UTC](https://discourse.julialang.org/t/reading-envi-style-multispectral-images-with-archgdal/79702 "2022-04-20T14:09:13Z")

</div>

I have several multispectral images with separate ENVI style header text files and I would like to load them with archGDAL.jl. I looked at the documentation of archGDAL.read() but I could not find how to do it. Maybe arc…

---

## [Plots.jl: Change formating of logscale ticks to normal formatting](https://discourse.julialang.org/t/plots-jl-change-formating-of-logscale-ticks-to-normal-formatting/79733)

<div class="topic-metadata">

**Author:** [@loonatick-src](https://discourse.julialang.org/u/loonatick-src)\
**Replies:** 2\
**Last updated:** [April 20, 2022, 12:08pm UTC](https://discourse.julialang.org/t/plots-jl-change-formating-of-logscale-ticks-to-normal-formatting/79733 "2022-04-20T12:08:33Z")

</div>

Consider the following plot using Plots pyplot() radix\_range\_xticks = -4:4 radix\_range = LinRange(-4, 4, 100) xs = 2.0 .^ radix\_range\_xticks pf = 50.0 bw = 25.0 roofline(x, bw, pf) = min(bw\*x, pf) ys = roofline.(xs, bw…

---

## [Warmstart in \`COSMO.jl\` with \`JuMP\`](https://discourse.julialang.org/t/warmstart-in-cosmo-jl-with-jump/79727)

<div class="topic-metadata">

**Author:** [@wujinq](https://discourse.julialang.org/u/wujinq)\
**Replies:** 1\
**Last updated:** [April 20, 2022, 9:58am UTC](https://discourse.julialang.org/t/warmstart-in-cosmo-jl-with-jump/79727 "2022-04-20T09:58:37Z")

</div>

I’m currently using COSMO.jl to do a large scale SDP with JuMP as the interface. Is it possible to store the status of model for example every 100000 iterations, and if max\_iter is reached, to start from the starting sta…

---

## [Calculate area of a raster](https://discourse.julialang.org/t/calculate-area-of-a-raster/79223)

<div class="topic-metadata">

**Author:** [@ahmoreira](https://discourse.julialang.org/u/ahmoreira)\
**Replies:** 4\
**Last updated:** [April 20, 2022, 8:30am UTC](https://discourse.julialang.org/t/calculate-area-of-a-raster/79223 "2022-04-20T08:30:43Z")

</div>

Hi all, I am trying to calculate the area of urban land cover only for cells in which 100% of the cell is urban cover. I am calculating the area as the product of the resolution times the number of cells, but it is a bi…

---

## [Manual RNN on GPU](https://discourse.julialang.org/t/manual-rnn-on-gpu/79580)

<div class="topic-metadata">

**Author:** [@kailukowiak](https://discourse.julialang.org/u/kailukowiak)\
**Replies:** 3\
**Last updated:** [April 20, 2022, 12:29am UTC](https://discourse.julialang.org/t/manual-rnn-on-gpu/79580 "2022-04-20T00:29:59Z")

</div>

I’m trying to redo this tutorial in Flux. The neural net runs fine on a CPU but when I try transferring it to the GPU I get an error which is caused by indexing a CUDA array. Here’s an MRE: ## using Flux using Flux: on…

---

## [Solving SDP problem is 100x slower than Matlab CVX](https://discourse.julialang.org/t/solving-sdp-problem-is-100x-slower-than-matlab-cvx/79704)

<div class="topic-metadata">

**Author:** [@biona001](https://discourse.julialang.org/u/biona001)\
**Replies:** 5\
**Last updated:** [April 19, 2022, 9:28pm UTC](https://discourse.julialang.org/t/solving-sdp-problem-is-100x-slower-than-matlab-cvx/79704 "2022-04-19T21:28:35Z")

</div>

I’m trying to solve the following SDP problem \\begin{align\*} maximize \\ & \\sum\_j s\_j\\\\ \\text{subject to } \\ & 2\\Sigma - diag(s) \\succeq 0,\\\\ & 0 \\le s\_j \\le 1 \\end{align\*} Julia’s SCS solver with Convex.jl is 100x sl…

---

## [Bayesian regression with parametrized basis functions in Turing.jl](https://discourse.julialang.org/t/bayesian-regression-with-parametrized-basis-functions-in-turing-jl/71725)

<div class="topic-metadata">

**Author:** [@mvsoom](https://discourse.julialang.org/u/mvsoom)\
**Replies:** 7\
**Last updated:** [April 19, 2022, 5:55pm UTC](https://discourse.julialang.org/t/bayesian-regression-with-parametrized-basis-functions-in-turing-jl/71725 "2022-04-19T17:55:14Z")

</div>

I’m trying to migrate a simple Bayesian model handwritten in Python to Turing.jl for a more modern approach. The model is a classical parametrized basis function expansion: a 1D (speech) signal f(t) is expanded into par…

---

## [Best practices for weak and uninformative priors in Turing](https://discourse.julialang.org/t/best-practices-for-weak-and-uninformative-priors-in-turing/79666)

<div class="topic-metadata">

**Author:** [@charperflow](https://discourse.julialang.org/u/charperflow)\
**Replies:** 5\
**Last updated:** [April 19, 2022, 5:29pm UTC](https://discourse.julialang.org/t/best-practices-for-weak-and-uninformative-priors-in-turing/79666 "2022-04-19T17:29:45Z")

</div>

I am working through a question for a class taught in R and I am trying to recreate the class examples in Julia. Now that we have started doing bayesian regression I have hit a bit of a learning curve, because I am doing…

---

## [Calculating crossentropy with varying sized vectors on GPU](https://discourse.julialang.org/t/calculating-crossentropy-with-varying-sized-vectors-on-gpu/79700)

<div class="topic-metadata">

**Author:** [@Casper](https://discourse.julialang.org/u/Casper)\
**Replies:** 0\
**Last updated:** [April 19, 2022, 4:15pm UTC](https://discourse.julialang.org/t/calculating-crossentropy-with-varying-sized-vectors-on-gpu/79700 "2022-04-19T16:15:08Z")

</div>

Hello, I’m attempting to calculate the cross entropy of a bunch of samples of different length. using CUDA entropy(π̂, π) =-sum(π .\* log.(π̂ )) m = 3 n = 3 w = ones(m) |\> gpu a1 = rand(Float32, m, n) a1 = softmax(a…

---

## [AWSS3.jl other cloud provider](https://discourse.julialang.org/t/awss3-jl-other-cloud-provider/79505)

<div class="topic-metadata">

**Author:** [@bienpierre](https://discourse.julialang.org/u/bienpierre)\
**Replies:** 2\
**Last updated:** [April 19, 2022, 3:22pm UTC](https://discourse.julialang.org/t/awss3-jl-other-cloud-provider/79505 "2022-04-19T15:22:47Z")

</div>

Hello, Is it possible to manage S3 buckets with AWSS3.jl from other cloud providers. Several cloud providers allow using aws cli ? Thanks for answering.

---

## [Application of DiffEqFlux to real data, with useful practical techniques](https://discourse.julialang.org/t/application-of-diffeqflux-to-real-data-with-useful-practical-techniques/79691)

<div class="topic-metadata">

**Author:** [@evolbio](https://discourse.julialang.org/u/evolbio)\
**Replies:** 0\
**Last updated:** [April 19, 2022, 1:33pm UTC](https://discourse.julialang.org/t/application-of-diffeqflux-to-real-data-with-useful-practical-techniques/79691 "2022-04-19T13:33:12Z")

</div>

Self-promotion … The following article may be a useful guide for fitting real time series data to ODE and NODE models with DiffEqFlux.jl. \[2204.07833\] Optimizing differential equations to fit data and predict outcomes T…

---

## [Replicating Matlab ode15s results in Julia](https://discourse.julialang.org/t/replicating-matlab-ode15s-results-in-julia/79280)

<div class="topic-metadata">

**Author:** [@BinaryBot](https://discourse.julialang.org/u/BinaryBot)\
**Replies:** 15\
**Last updated:** [April 19, 2022, 12:31pm UTC](https://discourse.julialang.org/t/replicating-matlab-ode15s-results-in-julia/79280 "2022-04-19T12:31:32Z")

</div>

I’m new to Julia and I have been working on ODE solvers for my first project .I recently got my hands onto some data which was generated using Matlab. I have been trying to replicate the results using the Julia Different…

---

## [Gaussian Processes and Infinite Neural Networks with Julia](https://discourse.julialang.org/t/gaussian-processes-and-infinite-neural-networks-with-julia/79677)

<div class="topic-metadata">

**Author:** [@bonevbs](https://discourse.julialang.org/u/bonevbs)\
**Replies:** 0\
**Last updated:** [April 19, 2022, 10:31am UTC](https://discourse.julialang.org/t/gaussian-processes-and-infinite-neural-networks-with-julia/79677 "2022-04-19T10:31:52Z")

</div>

I recently wanted to explore the connection between GPs and Infinite Neural Networks and took this as an excuse to become familiar with the Distributions package and Flux.jl. I thought others here might find this interes…

---

## [Use of intermediate variables to condition likelihood](https://discourse.julialang.org/t/use-of-intermediate-variables-to-condition-likelihood/78313)

<div class="topic-metadata">

**Author:** [@3mtmuhPnQv](https://discourse.julialang.org/u/3mtmuhPnQv)\
**Replies:** 1\
**Last updated:** [April 19, 2022, 8:58am UTC](https://discourse.julialang.org/t/use-of-intermediate-variables-to-condition-likelihood/78313 "2022-04-19T08:58:18Z")

</div>

I’d like to model prices as following a kind of stochastic difference equation, where jumps from one time step to the next follow p\_{t} = p\_{t - 1} + c \\exp(- b t\_0) \\Delta t + X(t) X ~ \\text{student\_t}(3, 0, \\sqrt{\\De…

---

## [How can I improve this parallel routine on mutually independent tasks?](https://discourse.julialang.org/t/how-can-i-improve-this-parallel-routine-on-mutually-independent-tasks/79669)

<div class="topic-metadata">

**Author:** [@zxjroger](https://discourse.julialang.org/u/zxjroger)\
**Replies:** 1\
**Last updated:** [April 19, 2022, 6:04am UTC](https://discourse.julialang.org/t/how-can-i-improve-this-parallel-routine-on-mutually-independent-tasks/79669 "2022-04-19T06:04:23Z")

</div>

I intend to solve N mutually independent tasks in parallel and each task is solved on a single CPU. I find an example code online and revised it as in the following MWE. I was wondering: Why there is a for loop over a…

---

## [Querying optimizer to set for another JuMP model](https://discourse.julialang.org/t/querying-optimizer-to-set-for-another-jump-model/79672)

<div class="topic-metadata">

**Author:** [@math\_opt](https://discourse.julialang.org/u/math_opt)\
**Replies:** 15\
**Last updated:** [April 19, 2022, 5:07am UTC](https://discourse.julialang.org/t/querying-optimizer-to-set-for-another-jump-model/79672 "2022-04-19T05:07:09Z")

</div>

I have a question (probably very naive though) on how to query an optimizer used for a particular model, so that I can use it in another JuMP model. For example, suppose user decides to use GLPK for the first model as sh…

---

## [Turing.jl identical likelihood but different \`MAP()\` result?](https://discourse.julialang.org/t/turing-jl-identical-likelihood-but-different-map-result/79671)

<div class="topic-metadata">

**Author:** [@jling](https://discourse.julialang.org/u/jling)\
**Replies:** 0\
**Last updated:** [April 19, 2022, 3:09am UTC](https://discourse.julialang.org/t/turing-jl-identical-likelihood-but-different-map-result/79671 "2022-04-19T03:09:58Z")

</div>

sample script: https://gist.github.com/Moelf/ab2f10ce45873bb91eaec60ace8388c0 you will need \]dev https://github.com/JuliaHEP/LiteHF.jl \]add Turing, Optim to compare, flip this line: https://github.com/JuliaHEP/LiteH…

---

## [Dummy variable change optimization result](https://discourse.julialang.org/t/dummy-variable-change-optimization-result/79646)

<div class="topic-metadata">

**Author:** [@dscheirman](https://discourse.julialang.org/u/dscheirman)\
**Replies:** 2\
**Last updated:** [April 18, 2022, 10:10pm UTC](https://discourse.julialang.org/t/dummy-variable-change-optimization-result/79646 "2022-04-18T22:10:20Z")

</div>

Hi, I have established a JuMP model (gurobi optimizer) and one constraint seems to be not respected: @constraint(obj,OP\_t == sum(sum(sum(((opex\_trans\[i,1,k,1\]\*Nsh\[i,1,k,t\]) + (opex\_trans\[i,1,k,2\]\*Np\[i,1,k,t\])) for i in…

---

## [Does GLM return an estimate for the error term (residual standard deviation)?](https://discourse.julialang.org/t/does-glm-return-an-estimate-for-the-error-term-residual-standard-deviation/79614)

<div class="topic-metadata">

**Author:** [@charperflow](https://discourse.julialang.org/u/charperflow)\
**Replies:** 2\
**Last updated:** [April 18, 2022, 1:48pm UTC](https://discourse.julialang.org/t/does-glm-return-an-estimate-for-the-error-term-residual-standard-deviation/79614 "2022-04-18T13:48:57Z")

</div>

I was under the impression that if I called lm(@formula(y~x),data) from the GLM.jl package was fitting a simple linear regression of the form: y=a+bx\_i+\\epsilon\_i The results from the regression being: \\hat{y}=\\hat{a}…

---

## [Using Makie for contour plot with only one edge](https://discourse.julialang.org/t/using-makie-for-contour-plot-with-only-one-edge/79603)

<div class="topic-metadata">

**Author:** [@Sciemon](https://discourse.julialang.org/u/Sciemon)\
**Replies:** 3\
**Last updated:** [April 17, 2022, 4:30pm UTC](https://discourse.julialang.org/t/using-makie-for-contour-plot-with-only-one-edge/79603 "2022-04-17T16:30:39Z")

</div>

I want a 2D contour plot which shows the boundary where “matrix” becomes larger than zero. I do: CairoMakie.contour(matrix,levels = \[0.0\]) The last step to make me happy would be to set the color of the boundary line. …

---

## [Help speeding up mixed integer SDP solve](https://discourse.julialang.org/t/help-speeding-up-mixed-integer-sdp-solve/23889)

<div class="topic-metadata">

**Author:** [@ericphanson](https://discourse.julialang.org/u/ericphanson)\
**Replies:** 3\
**Last updated:** [April 16, 2022, 5:26am UTC](https://discourse.julialang.org/t/help-speeding-up-mixed-integer-sdp-solve/23889 "2022-04-16T05:26:30Z")

</div>

I have a big mixed-integer SDP that I’d like to solve. Essentially my difficulty is it takes too long to solve-- using Pajarito, it doesn’t even finish the first MILP solve in 12 hours with access to 216 cores and 6 TB o…

---

## [Early Stopping Using Float Cutoff on Dual](https://discourse.julialang.org/t/early-stopping-using-float-cutoff-on-dual/79516)

<div class="topic-metadata">

**Author:** [@chelseas](https://discourse.julialang.org/u/chelseas)\
**Replies:** 2\
**Last updated:** [April 15, 2022, 11:27pm UTC](https://discourse.julialang.org/t/early-stopping-using-float-cutoff-on-dual/79516 "2022-04-15T23:27:31Z")

</div>

I cannot figure out how to do the following. I would like for JuMP to stop once the dual bound (lower bound for minimization) is above some value. OR similarly, once the dual bound (upper bound for maximization) is belo…

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

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