# Specific Domains

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

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

**Page:** 81

---

## [How to copy view of CuArray to Array efficiently?](https://discourse.julialang.org/t/how-to-copy-view-of-cuarray-to-array-efficiently/120960)

<div class="topic-metadata">

**Author:** [@0samuraiE](https://discourse.julialang.org/u/0samuraiE)\
**Replies:** 4\
**Last updated:** [October 6, 2024, 11:08am UTC](https://discourse.julialang.org/t/how-to-copy-view-of-cuarray-to-array-efficiently/120960 "2024-10-06T11:08:53Z")

</div>

Hi all. For some reason, I encountered a GPU allocation problem when I copied a CuArray view. I tried several copying methods, but the third result was much slower. Should I avoid copying the CuArray view to the host? …

---

## [How to use GPU acceleration to solve linear equation Ax=b](https://discourse.julialang.org/t/how-to-use-gpu-acceleration-to-solve-linear-equation-ax-b/120841)

<div class="topic-metadata">

**Author:** [@automaticvehiclerook](https://discourse.julialang.org/u/automaticvehiclerook)\
**Replies:** 9\
**Last updated:** [October 5, 2024, 5:35am UTC](https://discourse.julialang.org/t/how-to-use-gpu-acceleration-to-solve-linear-equation-ax-b/120841 "2024-10-05T05:35:07Z")

</div>

In my experiment,I want to solve the equation Ax=b on GPU to accelerate the process.I have tried the iterative solver from Kryvol.jl and cuSOLVER,but the GPU is still slower than CPU.Is that normal? using CUDA using Spa…

---

## [TIMES.jl prototype](https://discourse.julialang.org/t/times-jl-prototype/120932)

<div class="topic-metadata">

**Author:** [@olejandro](https://discourse.julialang.org/u/olejandro)\
**Replies:** 3\
**Last updated:** [October 5, 2024, 1:02am UTC](https://discourse.julialang.org/t/times-jl-prototype/120932 "2024-10-05T01:02:30Z")

</div>

Dear all, Together with other members of the TIMES community, I’ve been working on a small TIMES.jl prototype translating a small number of TIMES equations formulated in GAMS to Julia. Basically, this is a continuation…

---

## [I find R's DSL for defining regression a bit mind blowing. GLM.jl and other have tried to "adopt" them. What would an ideal DSL for defining regression look like in Julia?](https://discourse.julialang.org/t/i-find-rs-dsl-for-defining-regression-a-bit-mind-blowing-glm-jl-and-other-have-tried-to-adopt-them-what-would-an-ideal-dsl-for-defining-regression-look-like-in-julia/120786)

<div class="topic-metadata">

**Author:** [@xiaodai](https://discourse.julialang.org/u/xiaodai)\
**Replies:** 10\
**Last updated:** [October 4, 2024, 5:44pm UTC](https://discourse.julialang.org/t/i-find-rs-dsl-for-defining-regression-a-bit-mind-blowing-glm-jl-and-other-have-tried-to-adopt-them-what-would-an-ideal-dsl-for-defining-regression-look-like-in-julia/120786 "2024-10-04T17:44:45Z")

</div>

In R glm(x~-1+y+z) fits a GLM model without an intercept and using y and z as data, so there will be two coefficients from the model fit. I think GLM.jl has adopted the ~ DSL. But I’ve always found R’s DSL on regression…

---

## [Turing.jl Memory Issues](https://discourse.julialang.org/t/turing-jl-memory-issues/120925)

<div class="topic-metadata">

**Author:** [@willm](https://discourse.julialang.org/u/willm)\
**Replies:** 0\
**Last updated:** [October 4, 2024, 3:45pm UTC](https://discourse.julialang.org/t/turing-jl-memory-issues/120925 "2024-10-04T15:45:44Z")

</div>

Hi All, I designed a relatively simple hierarchical model for my research. To fit the model I initially wanted to use MCMC (NUTS). The other researchers in the lab tend to use R, so I started with rstan. Rstan worked pe…

---

## [Colour next to the label in the legend of a mesh3d in Plots.jl](https://discourse.julialang.org/t/colour-next-to-the-label-in-the-legend-of-a-mesh3d-in-plots-jl/120896)

<div class="topic-metadata">

**Author:** [@Valery](https://discourse.julialang.org/u/Valery)\
**Replies:** 1\
**Last updated:** [October 4, 2024, 9:33am UTC](https://discourse.julialang.org/t/colour-next-to-the-label-in-the-legend-of-a-mesh3d-in-plots-jl/120896 "2024-10-04T09:33:15Z")

</div>

Why the colour of the surfaces does not appear next to the text label in the legend box when I use the code below ? How to have it appear ? x = \[0, 0, 0, 0, 1, 1, 1, 1\] y = \[0, 0, 1, 1, 0, 0, 1, 1\] z = \[0, 1, 0, 1, 0…

---

## [Allow multiple datetime formats with CSV.jl](https://discourse.julialang.org/t/allow-multiple-datetime-formats-with-csv-jl/120891)

<div class="topic-metadata">

**Author:** [@Mastomaki](https://discourse.julialang.org/u/Mastomaki)\
**Replies:** 0\
**Last updated:** [October 4, 2024, 6:55am UTC](https://discourse.julialang.org/t/allow-multiple-datetime-formats-with-csv-jl/120891 "2024-10-04T06:55:30Z")

</div>

I would like to read CSV files with a time column where the format of the time can be one of several alternatives. CSV.file() does not automatically recognize the formats. Now I do like follows but is there are better w…

---

## [How does Dualization.jl solve a SOS but not compute it?](https://discourse.julialang.org/t/how-does-dualization-jl-solve-a-sos-but-not-compute-it/119966)

<div class="topic-metadata">

**Author:** [@ashefa](https://discourse.julialang.org/u/ashefa)\
**Replies:** 2\
**Last updated:** [October 4, 2024, 6:24am UTC](https://discourse.julialang.org/t/how-does-dualization-jl-solve-a-sos-but-not-compute-it/119966 "2024-10-04T06:24:43Z")

</div>

Hi, While solving a conic optimization problem, I needed to compute its dual, so I used Dualization.jl. However, it does not support the constraints modeled by SumOfSquares.jl, resulting in the following error: Constra…

---

## [Error on PEtab related to SciMLSensitivity package](https://discourse.julialang.org/t/error-on-petab-related-to-scimlsensitivity-package/120871)

<div class="topic-metadata">

**Author:** [@Nick1](https://discourse.julialang.org/u/Nick1)\
**Replies:** 0\
**Last updated:** [October 3, 2024, 9:19pm UTC](https://discourse.julialang.org/t/error-on-petab-related-to-scimlsensitivity-package/120871 "2024-10-03T21:19:59Z")

</div>

Hi, I am running a parameter estimation on PEtab and I get the following error when I use InterpolatingAdjoint(autojacvec = ReverseDiffVJP(true)) or InterpolatingAdjoint(autojacvec = EnzymeVJP())) type FakeIntegrator …

---

## [Makie - can I exclude some plots from automatic axis limit calculation?](https://discourse.julialang.org/t/makie-can-i-exclude-some-plots-from-automatic-axis-limit-calculation/120870)

<div class="topic-metadata">

**Author:** [@nchisholm](https://discourse.julialang.org/u/nchisholm)\
**Replies:** 1\
**Last updated:** [October 3, 2024, 8:56pm UTC](https://discourse.julialang.org/t/makie-can-i-exclude-some-plots-from-automatic-axis-limit-calculation/120870 "2024-10-03T20:56:55Z")

</div>

When plotting multiple plots into an Axis in Makie, is it possible to exclude some of them from being accounted for when automatic axis limit calculation is performed? Something like fig = Figure() ax = Axis(fig\[1,1\]) …

---

## [Implementing an equivalent of python's xarray (labeled nd-arrays) in julia](https://discourse.julialang.org/t/implementing-an-equivalent-of-pythons-xarray-labeled-nd-arrays-in-julia/1741)

<div class="topic-metadata">

**Author:** [@nbren12](https://discourse.julialang.org/u/nbren12)\
**Replies:** 15\
**Last updated:** [October 3, 2024, 3:29pm UTC](https://discourse.julialang.org/t/implementing-an-equivalent-of-pythons-xarray-labeled-nd-arrays-in-julia/1741 "2024-10-03T15:29:58Z")

</div>

One of the python packages I use the most is an excellent library called xarray, which is very convenient for handing multi-dimensional data. This package provides a DataArray object which associates a n-d array with nam…

---

## [Issue with gridlines when cliping is used in GMT's grdimage](https://discourse.julialang.org/t/issue-with-gridlines-when-cliping-is-used-in-gmts-grdimage/120844)

<div class="topic-metadata">

**Author:** [@k945](https://discourse.julialang.org/u/k945)\
**Replies:** 0\
**Last updated:** [October 3, 2024, 8:50am UTC](https://discourse.julialang.org/t/issue-with-gridlines-when-cliping-is-used-in-gmts-grdimage/120844 "2024-10-03T08:50:20Z")

</div>

How do I mask out the ocean when using GMT? Doing this using grdimage does not appear to be working. Here is the image without clipping applied However when I set cliping in grdimage to :land I get the following …

---

## [CuArray + GLMakie](https://discourse.julialang.org/t/cuarray-glmakie/52461)

<div class="topic-metadata">

**Author:** [@LaurentPlagne](https://discourse.julialang.org/u/LaurentPlagne)\
**Replies:** 14\
**Last updated:** [October 3, 2024, 8:55am UTC](https://discourse.julialang.org/t/cuarray-glmakie/52461 "2024-10-03T08:55:57Z")

</div>

Hi, I wonder if it would make sense (from the performance point of view) to be able to use GLMakie.jl with CuArrays from CUDA.jl ?

---

## [How to profile Julia MPI code?](https://discourse.julialang.org/t/how-to-profile-julia-mpi-code/57136)

<div class="topic-metadata">

**Author:** [@PolarizedPoutine](https://discourse.julialang.org/u/PolarizedPoutine)\
**Replies:** 17\
**Last updated:** [October 3, 2024, 12:16am UTC](https://discourse.julialang.org/t/how-to-profile-julia-mpi-code/57136 "2024-10-03T00:16:05Z")

</div>

How do people usually profile Julia MPI code to find bottlenecks and reduce communication overhead? As an MPI noob I thought there would be some profiler tools out (kinda like cuprof and NVIDIA Nsight) there but they se…

---

## [Does Julia recognize block, banded, and block-banded matrices?](https://discourse.julialang.org/t/does-julia-recognize-block-banded-and-block-banded-matrices/120821)

<div class="topic-metadata">

**Author:** [@rsenne](https://discourse.julialang.org/u/rsenne)\
**Replies:** 6\
**Last updated:** [October 2, 2024, 6:49pm UTC](https://discourse.julialang.org/t/does-julia-recognize-block-banded-and-block-banded-matrices/120821 "2024-10-02T18:49:24Z")

</div>

For a simple example of why I need this, let’s say I want to compute the entropy of a multivariate Gaussian. One can do this as follows: \\text{Entropy} \\propto -\\frac{1}{2}\*\\log|-H| where -H is the negative Hessian at …

---

## [Flux: combine two neural networks to run simultaneously?](https://discourse.julialang.org/t/flux-combine-two-neural-networks-to-run-simultaneously/120776)

<div class="topic-metadata">

**Author:** [@greatpet](https://discourse.julialang.org/u/greatpet)\
**Replies:** 4\
**Last updated:** [October 2, 2024, 12:21pm UTC](https://discourse.julialang.org/t/flux-combine-two-neural-networks-to-run-simultaneously/120776 "2024-10-02T12:21:34Z")

</div>

For example, I have an input which is a length-10 vector. I want the first neural network to output a length-3 vector based on the first 5 elements of the input, and the second neural network to output a length-4 vector …

---

## [Benders Subproblem Infeasibility](https://discourse.julialang.org/t/benders-subproblem-infeasibility/102622)

<div class="topic-metadata">

**Author:** [@bcoulier](https://discourse.julialang.org/u/bcoulier)\
**Replies:** 6\
**Last updated:** [October 1, 2024, 9:46pm UTC](https://discourse.julialang.org/t/benders-subproblem-infeasibility/102622 "2024-10-01T21:46:49Z")

</div>

Hello everyone, I’m currently working on an optimization problem where I’m using Benders decomposition to determine the placement of EV charging stations. The master problem determines the placement of charging stations…

---

## [Why does GLM.jl fail to fit my logistic regression model but scikit learn and R has no issues? I think I found the reason](https://discourse.julialang.org/t/why-does-glm-jl-fail-to-fit-my-logistic-regression-model-but-scikit-learn-and-r-has-no-issues-i-think-i-found-the-reason/118553)

<div class="topic-metadata">

**Author:** [@xiaodai](https://discourse.julialang.org/u/xiaodai)\
**Replies:** 5\
**Last updated:** [October 1, 2024, 2:43pm UTC](https://discourse.julialang.org/t/why-does-glm-jl-fail-to-fit-my-logistic-regression-model-but-scikit-learn-and-r-has-no-issues-i-think-i-found-the-reason/118553 "2024-10-01T14:43:54Z")

</div>

I have been trying to use GLM.jl to fit some logistic regression models and the fit failed. I don’t remember the exact error message but it was something about singularity. That annoyed me for a couple of years but I ne…

---

## [Problem with weight retrieval in Directed Graphs](https://discourse.julialang.org/t/problem-with-weight-retrieval-in-directed-graphs/120784)

<div class="topic-metadata">

**Author:** [@paul20](https://discourse.julialang.org/u/paul20)\
**Replies:** 3\
**Last updated:** [October 1, 2024, 2:41pm UTC](https://discourse.julialang.org/t/problem-with-weight-retrieval-in-directed-graphs/120784 "2024-10-01T14:41:44Z")

</div>

I’ve been using the SimpleWeightedDiGraph from the SimpleWeightedGraphs package for creating directed weighted graphs. Recently, I encountered a problem related to weight retrieval of edges. When initializing a SimpleWe…

---

## [My system is overdetermined. Are there tools to find the culprit equation?](https://discourse.julialang.org/t/my-system-is-overdetermined-are-there-tools-to-find-the-culprit-equation/119947)

<div class="topic-metadata">

**Author:** [@Datseris](https://discourse.julialang.org/u/Datseris)\
**Replies:** 4\
**Last updated:** [October 1, 2024, 1:29pm UTC](https://discourse.julialang.org/t/my-system-is-overdetermined-are-there-tools-to-find-the-culprit-equation/119947 "2024-10-01T13:29:57Z")

</div>

Recently I’ve updated a very big project to the latest version of MTK. I now get this warning: ┌ Warning: Initialization system is overdetermined. 1 equations for 0 unknowns. Initialization will default to using least …

---

## [Adding uncertainty on timestep of Bayesian ODE solve through \`saveat\`](https://discourse.julialang.org/t/adding-uncertainty-on-timestep-of-bayesian-ode-solve-through-saveat/120424)

<div class="topic-metadata">

**Author:** [@japhir](https://discourse.julialang.org/u/japhir)\
**Replies:** 9\
**Last updated:** [October 1, 2024, 5:11am UTC](https://discourse.julialang.org/t/adding-uncertainty-on-timestep-of-bayesian-ode-solve-through-saveat/120424 "2024-10-01T05:11:17Z")

</div>

I’d like to do a Bayesian inversion of parameters to an ODE model. I would like to add some uncertainty on the input time series. Is there any way to do this? I’ve tried to add a parameter tau ~ Uniform(-60, 60) and used…

---

## [CheckGradient In Manopt.jl](https://discourse.julialang.org/t/checkgradient-in-manopt-jl/120743)

<div class="topic-metadata">

**Author:** [@fimiller](https://discourse.julialang.org/u/fimiller)\
**Replies:** 0\
**Last updated:** [September 30, 2024, 11:25pm UTC](https://discourse.julialang.org/t/checkgradient-in-manopt-jl/120743 "2024-09-30T23:25:08Z")

</div>

I am performing an optimization problem and wish to make sure that the gradient I have is correct. When using the check\_gradient function, I find that the gradient is only correct some of the time, (about 90% of the tim…

---

## [Benders decomposition](https://discourse.julialang.org/t/benders-decomposition/120734)

<div class="topic-metadata">

**Author:** [@wahadhasan](https://discourse.julialang.org/u/wahadhasan)\
**Replies:** 1\
**Last updated:** [September 30, 2024, 8:21pm UTC](https://discourse.julialang.org/t/benders-decomposition/120734 "2024-09-30T20:21:28Z")

</div>

Hi there, can anybody tell me how to take the sub-gradient here in case of an (==) constraint or (\>=) constraint instead of standard format (\<=)? Is it the same procedure or need to convert it into standard format first? …

---

## [Copy metadata between DataFrames](https://discourse.julialang.org/t/copy-metadata-between-dataframes/107293)

<div class="topic-metadata">

**Author:** [@alex-s-gardner](https://discourse.julialang.org/u/alex-s-gardner)\
**Replies:** 2\
**Last updated:** [September 30, 2024, 7:32pm UTC](https://discourse.julialang.org/t/copy-metadata-between-dataframes/107293 "2024-09-30T19:32:45Z")

</div>

Is there a method to copy all of the metadata from one DataFrame to another?

---

## [(Dis)Advantages of quintic (fifth order) vs. cubic (third order) Hermite interpolation?](https://discourse.julialang.org/t/dis-advantages-of-quintic-fifth-order-vs-cubic-third-order-hermite-interpolation/120383)

<div class="topic-metadata">

**Author:** [@Marco\_Masi](https://discourse.julialang.org/u/Marco_Masi)\
**Replies:** 18\
**Last updated:** [September 30, 2024, 5:23pm UTC](https://discourse.julialang.org/t/dis-advantages-of-quintic-fifth-order-vs-cubic-third-order-hermite-interpolation/120383 "2024-09-30T17:23:51Z")

</div>

After solving an ordinary differential equation (ODE) of a dynamical system and obtaining the solution in time with its first and second derivatives, I fit the data points with a QuinticHermiteSpline(ddu, du, u, t). It a…

---

## [How to map rectangular lat/lon region on non-rectangular projection with GeoMakie?](https://discourse.julialang.org/t/how-to-map-rectangular-lat-lon-region-on-non-rectangular-projection-with-geomakie/103140)

<div class="topic-metadata">

**Author:** [@haffner](https://discourse.julialang.org/u/haffner)\
**Replies:** 5\
**Last updated:** [September 30, 2024, 4:06pm UTC](https://discourse.julialang.org/t/how-to-map-rectangular-lat-lon-region-on-non-rectangular-projection-with-geomakie/103140 "2024-09-30T16:06:04Z")

</div>

I’m trying to map a rectangular lat/lon polygon on a non-rectangular map projection using GeoMakie, When I use four corner points to define the region, GeoMakie draws straight lines between them. The example below show…

---

## [Quiver plots with PlotlyJS](https://discourse.julialang.org/t/quiver-plots-with-plotlyjs/120695)

<div class="topic-metadata">

**Author:** [@VivMendes](https://discourse.julialang.org/u/VivMendes)\
**Replies:** 6\
**Last updated:** [September 30, 2024, 3:50pm UTC](https://discourse.julialang.org/t/quiver-plots-with-plotlyjs/120695 "2024-09-30T15:50:57Z")

</div>

Hi. I am using a code written by @empet here to produce a quiver plot using PlotlyJS. It works very well if the ranges of the two variables are the same (in a system of differential equations). For example, the code belo…

---

## [How to make ODESystem with DataInterpolations latex output more readable?](https://discourse.julialang.org/t/how-to-make-odesystem-with-datainterpolations-latex-output-more-readable/120700)

<div class="topic-metadata">

**Author:** [@langestefan](https://discourse.julialang.org/u/langestefan)\
**Replies:** 7\
**Last updated:** [September 30, 2024, 9:57am UTC](https://discourse.julialang.org/t/how-to-make-odesystem-with-datainterpolations-latex-output-more-readable/120700 "2024-09-30T09:57:03Z")

</div>

I have defined my model in ModelingToolkit.jl as such: # define model parameters @parameters Ci Ria Aw # define sub-model components @named C\_i = Capacitor(C=Ci) # thermal capacity of interior @named R\_ia = Resistor…

---

## [How do I extract more than a pvalue HypothesisTests.jl?](https://discourse.julialang.org/t/how-do-i-extract-more-than-a-pvalue-hypothesistests-jl/120703)

<div class="topic-metadata">

**Author:** [@Rosejoycrocker](https://discourse.julialang.org/u/Rosejoycrocker)\
**Replies:** 1\
**Last updated:** [September 30, 2024, 7:50am UTC](https://discourse.julialang.org/t/how-do-i-extract-more-than-a-pvalue-hypothesistests-jl/120703 "2024-09-30T07:50:29Z")

</div>

Hi, I’m using HypothesisTests.jl to do a series of 2 Sample Kolmogorov Smirnov tests. This means I’m looping through a large series of data sets and applying the test KS\_object = ApproximateTwoSampleKSTest(data\_set1, da…

---

## [Why re-use factorization(A) rather than inv(A)](https://discourse.julialang.org/t/why-re-use-factorization-a-rather-than-inv-a/113234)

<div class="topic-metadata">

**Author:** [@photor](https://discourse.julialang.org/u/photor)\
**Replies:** 7\
**Last updated:** [September 30, 2024, 5:45am UTC](https://discourse.julialang.org/t/why-re-use-factorization-a-rather-than-inv-a/113234 "2024-09-30T05:45:00Z")

</div>

What is the advantage of first decomposing and then left dividing every time over inverting once and multiplying the inverse every time?

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

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