Questions tagged [julia]
Julia is a high-level, high-performance dynamic programming language for technical computing. It addresses the two-language problem by combining the ease of use of languages such as R and Python with the performance of C and Fortran.
32 questions
1
vote
2
answers
248
views
The Minimum Area Enclosing Circle
Given a set of points in 2D, how can I derive the minimum area circle which encloses all points?
Related questions (By @RodrigodeAzevedo):
How Do I Find the Minimum Area Ellipse that Encloses a Set ...
- 767
0
votes
1
answer
165
views
How to define the Jacobian of PDEs for use with ODE solvers -- Case: The incompressible Navier-Stokes equations
Question
To solve system of equations resulting from the semi-discretization of PDEs (e.g., using finite element methods), one can use for example SciML.ODEFunction ...
- 269
1
vote
1
answer
199
views
Using abs2 function for automatic differentiation
In a Julia code to approximate the characteristic function of $]-\infty,0]$ (being $1$ on this interval and zero elsewhere) I saw this implementation :
...
- 498
3
votes
0
answers
101
views
Computing Lyapunov spectrum numerically in 2 dimensions: code yields exponents with positive sum?
I implemented a QR-decomposition based algorithm for estimating the Lyapunov spectrum for the following two-dimensional discrete map, described in Sprott (1993):
$x_{n+1}=a_0+a_1x_n+a_2x^2_n+a_3x_ny_n+...
- 51
1
vote
0
answers
83
views
Dogleg Trust Region Method Fail
I'm implementing in Julia the Dogleg Trust Region Method.
I've got this:
...
0
votes
0
answers
77
views
Are there any efficient numerical methods to solve the recurrence relation of a function？
First, the recurrence relation is:
$$\pi_{k+1}(\omega)=(\omega-\alpha_k)\pi_{k}(\omega)-\beta_k\pi_{k-1}(\omega),$$
where
$$\alpha_k=\frac{\int_0^\Lambda\omega\pi_{k}^2(\omega)h^2(\omega)\text{d}\...
- 31
0
votes
1
answer
268
views
Auto differentiation with JAX in python and ForwardDiff.jl in Julia give matrices with different values but same structure. Are the results right?
Using ForwardDiff in Julia gives me the following values for 2 matrices
I have tried to use JAX to auto differentiate the same function in Python but I am getting
the following results for A and B ...
- 11
2
votes
1
answer
83
views
Solving $(I-Q)x={\bf 1}$ for sub-stochastic sparse $Q$ of dimension 5M $\times$ 5M
I have a (right) sub-stochastic CSC sparse matrix $Q$ of dimension 5 million, with 200 million nonzero entries, which is a nonzero percentage of 0.0008%, so it is indeed extremely sparse. It is not ...
- 523
5
votes
3
answers
502
views
Packages suitable for numerical optimization of functions with discontinuous gradient at the point of minimum
Are there packages for numerical optimization in julia or python, or in any other system for scientific computing, capable of taking into account the discontinuity of gradient at the minimum point? ...
- 153
0
votes
0
answers
96
views
non zero dirichlet boundary condition entered in weak form
i am trying to write a julia code for linear elasticity in my case i dont have body force and traction but there is a nonzero drichlet bc(ubc) i want to engage the bc in weak form in linear part. is ...
- 1
0
votes
1
answer
81
views
Optimization of the log-absolute: reformulating to DCP-compliant on Julia
I am trying to reformulate this optimization problem in order to get a DCP-complaint expression on Julia (I am using the ...
- 205
1
vote
0
answers
99
views
Implementation of integration schemes for ordinary differential equations in Python and peformance comparison
I look for a book/manual where I can find implementations of different integration schemes for ordinary differential equations (like 4-th order Runge-Kutta) in Python with Numba.
To be more specific, ...
- 181
2
votes
1
answer
461
views
solving a Algebraic Differential Equation in Julia using modelingToolKit.JL
I'm trying to solve a differential algebraic equation in Julia's modelingTookKit.JL, where the vector field has the form f(X) = 0.
I found an example of a DAE in the below link
modelingToolkit.JL DAE](...
- 53
3
votes
1
answer
545
views
Solving DAE in Julia using GPUs
I'm trying to solve a Differential Algebraic Equation (DAE) in Julia which is very computationally expensive using GPUs. I'm brand new to Julia and don't have much experience coding with GPUs. The ...
- 53
6
votes
1
answer
616
views
PETSc-like library for Julia
I want to build an application for Material Point Method (and probably other meshfree methods too) in Julia and I am looking for library for direct and iterative solvers that can help me with it. One ...
- 163