# Specific Domains

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

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

**Page:** 508

---

## [Some trick in NLconstraint to avoid NaN or Inf for zero initial condition ?](https://discourse.julialang.org/t/some-trick-in-nlconstraint-to-avoid-nan-or-inf-for-zero-initial-condition/9831)

<div class="topic-metadata">

**Author:** [@cfrodriguez](https://discourse.julialang.org/u/cfrodriguez)\
**Replies:** 3\
**Last updated:** [February 26, 2019, 2:33pm UTC](https://discourse.julialang.org/t/some-trick-in-nlconstraint-to-avoid-nan-or-inf-for-zero-initial-condition/9831 "2019-02-26T14:33:57Z")

</div>

Hello all, I am new using JuMP. In my code, I added a thermodynamic model (UNIQUAC). When xNT1 = 0.001 it work well. But when the value is almost zero (1e-11) it doesn’t work and JuMP give me “Warning: Cutting back alph…

---

## [Using CUSOLVER in CuArrays.jl](https://discourse.julialang.org/t/using-cusolver-in-cuarrays-jl/21020)

<div class="topic-metadata">

**Author:** [@eliassno](https://discourse.julialang.org/u/eliassno)\
**Replies:** 12\
**Last updated:** [February 26, 2019, 2:21pm UTC](https://discourse.julialang.org/t/using-cusolver-in-cuarrays-jl/21020 "2019-02-26T14:21:35Z")

</div>

I’m writing a GPU-version of gesv!, using CuArrays.jl. The following works: using CuArrays, LinearAlgebra, Test function gpugesv!(A,b) A, ipiv = CuArrays.CUSOLVER.getrf!(A) CuArrays.CUSOLVER.getrs!('N'…

---

## [Smallest Eigenvalue of a Moler Matrix \[Again\]](https://discourse.julialang.org/t/smallest-eigenvalue-of-a-moler-matrix-again/21112)

<div class="topic-metadata">

**Author:** [@Hans\_W\_Borchers](https://discourse.julialang.org/u/Hans_W_Borchers)\
**Replies:** 6\
**Last updated:** [February 26, 2019, 12:43pm UTC](https://discourse.julialang.org/t/smallest-eigenvalue-of-a-moler-matrix-again/21112 "2019-02-26T12:43:12Z")

</div>

I am referring to help I received from Andreas Noack 10 months ago. https://discourse.julialang.org/t/smallest-eigenvalue-of-a-moler-matrix/10510 I wanted to compute the smallest eigenvalue of Moler matrices with high …

---

## [Flux not able to resolve strong observation, based on standard example that works](https://discourse.julialang.org/t/flux-not-able-to-resolve-strong-observation-based-on-standard-example-that-works/21177)

<div class="topic-metadata">

**Author:** [@ianshmean](https://discourse.julialang.org/u/ianshmean)\
**Replies:** 2\
**Last updated:** [February 25, 2019, 6:09pm UTC](https://discourse.julialang.org/t/flux-not-able-to-resolve-strong-observation-based-on-standard-example-that-works/21177 "2019-02-25T18:09:14Z")

</div>

I don’t understand why the housing example found here, which functions as expected, is not able to train on alternative data where I’m intentionally including the target in the observations to make it really easy. using…

---

## ["A Comparison of Three Programming Languages for a Full-Fledged Next-Generation Sequencing Tool" selects...Go?](https://discourse.julialang.org/t/a-comparison-of-three-programming-languages-for-a-full-fledged-next-generation-sequencing-tool-selects-go/21126)

<div class="topic-metadata">

**Author:** [@jballanc](https://discourse.julialang.org/u/jballanc)\
**Replies:** 5\
**Last updated:** [February 25, 2019, 5:40pm UTC](https://discourse.julialang.org/t/a-comparison-of-three-programming-languages-for-a-full-fledged-next-generation-sequencing-tool-selects-go/21126 "2019-02-25T17:40:09Z")

</div>

Was wondering if anyone has had a chance to look at this paper posted to bioRxiv yesterday? I’ve only just briefly skimmed the paper and associated code, but I find the results…interesting. The authors appear to have tak…

---

## [Getting the optimal solution](https://discourse.julialang.org/t/getting-the-optimal-solution/21167)

<div class="topic-metadata">

**Author:** [@Rajat\_Sanyal](https://discourse.julialang.org/u/Rajat_Sanyal)\
**Replies:** 2\
**Last updated:** [February 25, 2019, 2:22pm UTC](https://discourse.julialang.org/t/getting-the-optimal-solution/21167 "2019-02-25T14:22:13Z")

</div>

Let us consider the following optimization problem: using JuMP, Ipopt function Opti\_model(A,b) m = Model(with\_optimizer(Ipopt.Optimizer)) @variable(m,x\[1:5\],lower\_bound=0,upper\_bound=1) @objective(m,Max,sum(…

---

## [SVD of Matrix w/ updated columns](https://discourse.julialang.org/t/svd-of-matrix-w-updated-columns/21048)

<div class="topic-metadata">

**Author:** [@DrPapa](https://discourse.julialang.org/u/DrPapa)\
**Replies:** 7\
**Last updated:** [February 22, 2019, 5:26pm UTC](https://discourse.julialang.org/t/svd-of-matrix-w-updated-columns/21048 "2019-02-22T17:26:20Z")

</div>

I have a M\\times N matrix A where N\>\>M. I need to compute the SVD of sub-matrices using a sliding window over the columns, i.e. svd(A\[:,1:k\]), svd(A\[:,2:k+1\]), svd(A\[:,3:k+2\])… Does anybody know of an efficient way to u…

---

## [Flux custom loss function](https://discourse.julialang.org/t/flux-custom-loss-function/21127)

<div class="topic-metadata">

**Author:** [@jstrube](https://discourse.julialang.org/u/jstrube)\
**Replies:** 3\
**Last updated:** [February 24, 2019, 9:45pm UTC](https://discourse.julialang.org/t/flux-custom-loss-function/21127 "2019-02-24T21:45:11Z")

</div>

I’m trying to run a straightforward regression problem through Flux. I’ve modified the loss function from the example, from a simple mse to something like a chi-square. The problem has to do with the calibration of a cal…

---

## [Rank-revealing QR/SVD Decompositions](https://discourse.julialang.org/t/rank-revealing-qr-svd-decompositions/21155)

<div class="topic-metadata">

**Author:** [@Hans\_W\_Borchers](https://discourse.julialang.org/u/Hans_W_Borchers)\
**Replies:** 0\
**Last updated:** [February 24, 2019, 8:27pm UTC](https://discourse.julialang.org/t/rank-revealing-qr-svd-decompositions/21155 "2019-02-24T20:27:23Z")

</div>

Are there implementations of rank-revealing QR/SVD decompositions, resp. QR and SVD updating / downdating mechanisms in any Julia package? See for instance: \[1\] Li and Zheng, A Rank-revealing Method with Updating, Down…

---

## [VegaLite rescale data plots](https://discourse.julialang.org/t/vegalite-rescale-data-plots/21143)

<div class="topic-metadata">

**Author:** [@nik.ltin](https://discourse.julialang.org/u/nik.ltin)\
**Replies:** 0\
**Last updated:** [February 24, 2019, 1:32pm UTC](https://discourse.julialang.org/t/vegalite-rescale-data-plots/21143 "2019-02-24T13:32:54Z")

</div>

Hello all, I’m using VegaLite for bar plots with the @vlplot(:bar...) command. I’m looking at a dataset that I’d like to rescale to go from 1 to 100. I tried adding in the command julia\> y={:df1, scale={domain=\[0,100\]…

---

## [VegaLite rescale plots domain](https://discourse.julialang.org/t/vegalite-rescale-plots-domain/21120)

<div class="topic-metadata">

**Author:** [@nik.ltin](https://discourse.julialang.org/u/nik.ltin)\
**Replies:** 0\
**Last updated:** [February 23, 2019, 7:44pm UTC](https://discourse.julialang.org/t/vegalite-rescale-plots-domain/21120 "2019-02-23T19:44:23Z")

</div>

Hello all, I’m using VegaLite for bar plots with the @vlplot(:bar...) command. I’m looking at a dataset that I’d like to rescale to go from 1 to 100. I tried adding in the command y={:df1, scale={domain=\[0,100\], nice=fal…

---

## [Storing mouse click event coordinates in PyPlot window](https://discourse.julialang.org/t/storing-mouse-click-event-coordinates-in-pyplot-window/21117)

<div class="topic-metadata">

**Author:** [@Carsten](https://discourse.julialang.org/u/Carsten)\
**Replies:** 1\
**Last updated:** [February 23, 2019, 9:01pm UTC](https://discourse.julialang.org/t/storing-mouse-click-event-coordinates-in-pyplot-window/21117 "2019-02-23T21:01:58Z")

</div>

I want to get access to the cursor positions within a PyPlot window and store all clicked values into arrays. I was able to print out the x,y values but I failed to store them into arrays. The code first I used is: fu…

---

## [Name styles of variables, parameters, sets, etc...](https://discourse.julialang.org/t/name-styles-of-variables-parameters-sets-etc/21089)

<div class="topic-metadata">

**Author:** [@haokuny](https://discourse.julialang.org/u/haokuny)\
**Replies:** 1\
**Last updated:** [February 24, 2019, 12:30am UTC](https://discourse.julialang.org/t/name-styles-of-variables-parameters-sets-etc/21089 "2019-02-24T00:30:55Z")

</div>

When building a JuMP model, one thing I feel very important is to distinguish between variables and parameters in each equation. It’s important because, for example, the term “var. \* par.” is linear but the term “var. \* …

---

## [Help with my MINLP](https://discourse.julialang.org/t/help-with-my-minlp/19927)

<div class="topic-metadata">

**Author:** [@muriloasouza](https://discourse.julialang.org/u/muriloasouza)\
**Replies:** 8\
**Last updated:** [February 23, 2019, 6:59pm UTC](https://discourse.julialang.org/t/help-with-my-minlp/19927 "2019-02-23T18:59:52Z")

</div>

Hello, sorry for this long post and problem… but if any help is trully appreciated. I need help with my MINLP, since the results i’m getting is not what i was expecting… This is an electrical engineering problem, where …

---

## [Sum dimensions of matrix](https://discourse.julialang.org/t/sum-dimensions-of-matrix/21098)

<div class="topic-metadata">

**Author:** [@nik.ltin](https://discourse.julialang.org/u/nik.ltin)\
**Replies:** 11\
**Last updated:** [February 23, 2019, 2:25pm UTC](https://discourse.julialang.org/t/sum-dimensions-of-matrix/21098 "2019-02-23T14:25:20Z")

</div>

Hi there, A bit new to Julia so this might be a very basic question: I have a file with a 3 dimensional array (x,y,z), and I’d like to sum certain values from the z dimension to later on pull from and graph. For example…

---

## [Make changes of the new model effecting the old one](https://discourse.julialang.org/t/make-changes-of-the-new-model-effecting-the-old-one/21092)

<div class="topic-metadata">

**Author:** [@haokuny](https://discourse.julialang.org/u/haokuny)\
**Replies:** 2\
**Last updated:** [February 23, 2019, 6:17am UTC](https://discourse.julialang.org/t/make-changes-of-the-new-model-effecting-the-old-one/21092 "2019-02-23T06:17:09Z")

</div>

Let’s say I have built a model called mymodel and I want to make some extension, meaning to add some new variables and constraints, but I want to do it to my new model. So I mymodelnew = mymodel. Then the problem is, wh…

---

## [Distributed nested in Pmap](https://discourse.julialang.org/t/distributed-nested-in-pmap/21081)

<div class="topic-metadata">

**Author:** [@Alessio\_Bellisomi](https://discourse.julialang.org/u/Alessio_Bellisomi)\
**Replies:** 2\
**Last updated:** [February 22, 2019, 5:54pm UTC](https://discourse.julialang.org/t/distributed-nested-in-pmap/21081 "2019-02-22T17:54:59Z")

</div>

Hello, I am fairly new to Julia, I am trying to understand a little better how the workload is distributed, so I prepared this little script to help me out. An external pmap splits the work between different processes,…

---

## [UpwindOperator from DiffEqOperators.jl](https://discourse.julialang.org/t/upwindoperator-from-diffeqoperators-jl/21053)

<div class="topic-metadata">

**Author:** [@boutor2](https://discourse.julialang.org/u/boutor2)\
**Replies:** 2\
**Last updated:** [February 22, 2019, 8:38am UTC](https://discourse.julialang.org/t/upwindoperator-from-diffeqoperators-jl/21053 "2019-02-22T08:38:36Z")

</div>

I have been playing around with DifferentialEquations.jl, mainly for solving reaction diffusion PDE’s. I now wish to add an advection term to my model. For this purpose I discovered DiffEqOperators.jl. UpwindOperator sou…

---

## [Initializing Sparse Matrices with CuArrays.jl](https://discourse.julialang.org/t/initializing-sparse-matrices-with-cuarrays-jl/21062)

<div class="topic-metadata">

**Author:** [@carlomontec](https://discourse.julialang.org/u/carlomontec)\
**Replies:** 4\
**Last updated:** [February 22, 2019, 7:05am UTC](https://discourse.julialang.org/t/initializing-sparse-matrices-with-cuarrays-jl/21062 "2019-02-22T07:05:41Z")

</div>

Hello, total newbie with GPU computing here. I have seen that is possible to use sparse matrices with CuArrays.jl, but I have not found a way to initialize them. I tried to with the example described at https://github.…

---

## [Use Symplectic Integrator for ODE for time reversability constraint](https://discourse.julialang.org/t/use-symplectic-integrator-for-ode-for-time-reversability-constraint/20955)

<div class="topic-metadata">

**Author:** [@mleprovost](https://discourse.julialang.org/u/mleprovost)\
**Replies:** 6\
**Last updated:** [February 21, 2019, 8:39pm UTC](https://discourse.julialang.org/t/use-symplectic-integrator-for-ode-for-time-reversability-constraint/20955 "2019-02-21T20:39:59Z")

</div>

Hi everyone, I would like to use Symplectic Integrators such as VelocityVerlet or VerletLeapfrog from DifferentialEquations.jl to ensure the time-reversibility of the equation of motion. Nonetheless, my problem can no…

---

## [About the optimization problem choice in JuMP v0.19](https://discourse.julialang.org/t/about-the-optimization-problem-choice-in-jump-v0-19/21018)

<div class="topic-metadata">

**Author:** [@ErwanPiR](https://discourse.julialang.org/u/ErwanPiR)\
**Replies:** 8\
**Last updated:** [February 21, 2019, 2:20pm UTC](https://discourse.julialang.org/t/about-the-optimization-problem-choice-in-jump-v0-19/21018 "2019-02-21T14:20:29Z")

</div>

I hope this isn’t a dumb question… With the precedent version of JuMP, to define the kind of optimization problem (LP, MIP,…), the syntax was Model(solver=GLPKSolverMIP(args…)) or Model(solver=GLPKSolverLP()) I …

---

## [Package for genome graphics?](https://discourse.julialang.org/t/package-for-genome-graphics/20550)

<div class="topic-metadata">

**Author:** [@aalexandersson](https://discourse.julialang.org/u/aalexandersson)\
**Replies:** 11\
**Last updated:** [February 21, 2019, 1:27am UTC](https://discourse.julialang.org/t/package-for-genome-graphics/20550 "2019-02-21T01:27:14Z")

</div>

Are there any graphing capabilities in Julia 1.x specifically for genomics (“genome graphs”), for example a genome viewer? Ben Ward illustrated the Julia package GenomeBrowsers in his presentation at JuliaCon 2018. It h…

---

## [Confidence Interval for certain Contrast and model residual](https://discourse.julialang.org/t/confidence-interval-for-certain-contrast-and-model-residual/20339)

<div class="topic-metadata">

**Author:** [@PharmCat](https://discourse.julialang.org/u/PharmCat)\
**Replies:** 13\
**Last updated:** [February 20, 2019, 6:13pm UTC](https://discourse.julialang.org/t/confidence-interval-for-certain-contrast-and-model-residual/20339 "2019-02-20T18:13:44Z")

</div>

I have a model like this: ols = lm(@formula(LnVar ~ Seq+Per+Trt+Subj), df) I try to get confidence interval for Trt contrast, LSM and model residual. How to obtain it? In R you just do confint(lmobj, c("Trt"), level=…

---

## [Recommended way to do work stealing in Julia 1.0?](https://discourse.julialang.org/t/recommended-way-to-do-work-stealing-in-julia-1-0/15573)

<div class="topic-metadata">

**Author:** [@orenbenkiki](https://discourse.julialang.org/u/orenbenkiki)\
**Replies:** 6\
**Last updated:** [February 20, 2019, 1:48pm UTC](https://discourse.julialang.org/t/recommended-way-to-do-work-stealing-in-julia-1-0/15573 "2019-02-20T13:48:02Z")

</div>

I have a Python system which is a mixture of C++ code, pandas/numpy code running on top of BLAS/MKL, and plain Python code binding the whole thing together. “It ain’t a pretty sight”. It can be viewed as a nested tree (…

---

## [Keyword arguments in PlotRecipes](https://discourse.julialang.org/t/keyword-arguments-in-plotrecipes/20780)

<div class="topic-metadata">

**Author:** [@moesphere](https://discourse.julialang.org/u/moesphere)\
**Replies:** 4\
**Last updated:** [February 20, 2019, 7:45am UTC](https://discourse.julialang.org/t/keyword-arguments-in-plotrecipes/20780 "2019-02-20T07:45:52Z")

</div>

Following minimal example: struct PlotType end @recipe function f(::PlotType; title="Title", ylabel="ylabel") x := \[1,2,3\] y := rand(3,5) end The title is plotted fine, I also can change the title when calling …

---

## [Error with Ipopt ](https://discourse.julialang.org/t/error-with-ipopt/20989)

<div class="topic-metadata">

**Author:** [@Y\_Li](https://discourse.julialang.org/u/Y_Li)\
**Replies:** 2\
**Last updated:** [February 20, 2019, 1:32am UTC](https://discourse.julialang.org/t/error-with-ipopt/20989 "2019-02-20T01:32:35Z")

</div>

Hi, I am using JuMP and IPOPT to solve an non-linear optimization, below is my code using JuMP using Ipopt model =Model(solver=IpoptSolver()) @variable(model, x ) @variable(model, y ) @variable(model, z ) #define our ob…

---

## [ADAM optimizer error in Flux](https://discourse.julialang.org/t/adam-optimizer-error-in-flux/20947)

<div class="topic-metadata">

**Author:** [@microlifecc](https://discourse.julialang.org/u/microlifecc)\
**Replies:** 1\
**Last updated:** [February 19, 2019, 6:01pm UTC](https://discourse.julialang.org/t/adam-optimizer-error-in-flux/20947 "2019-02-19T18:01:23Z")

</div>

Hi! I am defining a neural network as: n\_hidden = 20 model = Chain(Dense(n\_inputs, n\_hidden, relu), Dense(n\_hidden, n\_outputs, identity), softmax) L(x,y) = Flux.crossentropy(model(x), y) ps = Flux.params(…

---

## [GLM.jl with unknown column names](https://discourse.julialang.org/t/glm-jl-with-unknown-column-names/20692)

<div class="topic-metadata">

**Author:** [@george](https://discourse.julialang.org/u/george)\
**Replies:** 4\
**Last updated:** [February 19, 2019, 5:47pm UTC](https://discourse.julialang.org/t/glm-jl-with-unknown-column-names/20692 "2019-02-19T17:47:29Z")

</div>

Hello, I’m using Julia for some basic statistical analysis. I’d like to run a simple linear regression on a small dataset using the GLM.jl package. My data looks something like this: │ Row │ stat.id12312.somemetadata …

---

## [JuMP: Catch error if solver stops early due to network error or running out of memory ?](https://discourse.julialang.org/t/jump-catch-error-if-solver-stops-early-due-to-network-error-or-running-out-of-memory/20962)

<div class="topic-metadata">

**Author:** [@Stuart\_Rogers](https://discourse.julialang.org/u/Stuart_Rogers)\
**Replies:** 1\
**Last updated:** [February 19, 2019, 1:32am UTC](https://discourse.julialang.org/t/jump-catch-error-if-solver-stops-early-due-to-network-error-or-running-out-of-memory/20962 "2019-02-19T01:32:55Z")

</div>

I am solving a large MILP by calling Gurobi through JuMP v0.18. Gurobi issues an error, which I believe is due to it running out of RAM. Is it possible for JuMP to catch the error and collect the best known solution?

---

## [JuMP v0.19 has been released](https://discourse.julialang.org/t/jump-v0-19-has-been-released/20878)

<div class="topic-metadata">

**Author:** [@blegat](https://discourse.julialang.org/u/blegat)\
**Replies:** 6\
**Last updated:** [February 18, 2019, 8:22pm UTC](https://discourse.julialang.org/t/jump-v0-19-has-been-released/20878 "2019-02-18T20:22:59Z")

</div>

We are thrilled to announce the release of JuMP v0.19. Work on this release started nearly two years ago at the first JuMP-dev meetup. As introduced in an earlier post and at the second JuMP-dev meetup, JuMP has undergo…

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

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