10 votes

How are scientific computing workflows faring on Apple's M1 hardware

The way I have been measuring whether the eco-system is ready is how things are going with the transition for the homebrew package manager. They have been carefully documenting the progress of ...
Kyle Mandli's user avatar
9 votes

Recommended language/environment for large scale semi-continuous biological models

You should consider giving Julia a try. Let me explain what's going on in the design space right now that would be of interest to you. Full disclosure I am the lead developer of JuliaDiffEq. ...
Chris Rackauckas's user avatar
8 votes

How are scientific computing workflows faring on Apple's M1 hardware

From this, it looks like there is no functional native Fortran compiler yet. If that is really the case, things look bleak. Almost anything that uses linear algebra includes some Fortran code (Lapack),...
Federico Poloni's user avatar
6 votes
Accepted

The difference between mkl_intel_lp64 vs mkl_gf_lp64 in a numerical reproducibility issue with Intel MKL

The various Fortran standards allow a lot of compiler dependent behaviour in terms of function binary interfaces when being called with "complicated" data types such as Fortran90 style arrays and ...
origimbo's user avatar
  • 2,249
5 votes

How are scientific computing workflows faring on Apple's M1 hardware

Iain Sandoe has been working on porting both gcc and gFortran to this architecture. Based on this, François-Xavier Coudert has created an experimental gFortran release for the M1. I provide a few ...
user1677899's user avatar
3 votes

Numerical solving Lotka-Volterra ODE in R

I did not yet succeed to use the sde.sim() function of the SDE package, however, I succeeded to solve the system (with and without noise) using the suggestions of ...
user213544's user avatar
3 votes
Accepted

Solving an SDE with time-dependent parameter in R

The parameter can be any type, so here I pass in a time-dependent function for p and use it in the differential equation: ...
Chris Rackauckas's user avatar
2 votes
Accepted

How to invert a lagrangian polynomial

Without knowing anything about inverse interpolation, this problem should be a straight-forward numerical task given the information about $g(x)$. Firstly, $g(x)$ is said to be monotonously ...
phil's user avatar
  • 54
2 votes

efficient way to compute lagrange polynomials

If you want to know more about efficient numerical procedures to compute Lagrange polynomials, I recommend the following reference. J.-P. Berrut and L. N. Trefethen, “Barycentric Lagrange ...
Nick C.'s user avatar
  • 188
2 votes

(Lack of) Availability of Finite-Difference library for simple 2D PDEs

I know this is a little bit late, but there's a third party library in Python which solve PDEs using FD, it is called findiff. Check it out.
Heitor Nunes Rosa's user avatar
2 votes

(Lack of) Availability of Finite-Difference library for simple 2D PDEs

It seems hard to write a general enough FD library thas has wide applicability, since FD methods are not as easy to write for general domains, unlike FEM which uses unstructured grids, for which there ...
cfdlab's user avatar
  • 3,028
2 votes
Accepted

Recommended language/environment for large scale semi-continuous biological models

We decided to go with C++, because it is free is fast has a freely available and widely supported numerical integration library (boost/odeint) is easy to integrate with our larger simulation driver/...
Simon Woodward's user avatar
2 votes

Solve discontinuous ODE with lsode

If you speak of discontinuous ODE, one would think of a discontinuous right side function. If the discontinuity is only a jump, and only depends on time, then using the ...
Lutz Lehmann's user avatar
  • 6,064
1 vote
Accepted

How to initialize Eigen C++ parameters within for-loop?

You used Eigen::MatrixXd::Zero for initialising BetaH. Why not follow the same for initialising the variables inside the loop?
Chenna K's user avatar
  • 934
1 vote

Problem of half-planes intersection

This is a textbook example of a continuous optimization problem: \begin{align} \min_{x} \quad &f(x)\\ \text{subject to }& g_i(x) \leq 0 \qquad i=1\,...\,m \end{align} With the objective ...
Yann's user avatar
  • 159
1 vote

Draw contour line to represent multiple contours

Without knowing about the expected behavior of the data I would compute a interpolation over a common parameter space and then average the interpolated values
nicoguaro's user avatar
  • 8,500
1 vote

R function or package for carrying out maximum likelihood techniques in random effect models

I suggest trying Bound Optimization By Quadratic Approximation which is sometimes used for random-effect modeling.
Anton Menshov's user avatar
  • 8,672
1 vote
Accepted

reading and processing large GRIB files

I would suggest cdo. CDO can read GRIB files (though I think it only supports GRIB2 files) and has a host of options for computing means, subsetting, regridding etc. If you add the ...
Sam Hatfield's user avatar
1 vote

Sparse matrix inverse with reduced bandwidth

Three things immediately come to mind: R might not take advantage of sparsity when using the solve command to compute the inverse of a matrix. Usually, the inverse ...
Daniel Shapero's user avatar

Only top scored, non community-wiki answers of a minimum length are eligible