solv_MKL_sparse#

MCEq.solvers.solv_MKL_sparse(nsteps, dX, rho_inv, int_m, dec_m, phi, grid_idcs)[source]#

Intel MKL sparse BLAS implementation of forward-euler integration.

Function requires that the path to the MKL runtime library libmkl_rt.[so/dylib] defined in the config file.

Parameters:

nsteps (int) – number of integration steps
dX (numpy.array() [nsteps]) – vector of step-sizes \(\Delta X_i\) in g/cm**2
rho_inv (numpy.array() [nsteps]) – vector of density values \(\frac{1}{\rho(X_i)}\)
int_m (numpy.array()) – interaction matrix (1) in dense or sparse representation
dec_m (numpy.array()) – decay matrix (2) in dense or sparse representation
phi (numpy.array()) – initial state vector \(\Phi(X_0)\)
grid_idcs (list) – indices at which longitudinal solutions have to be saved.

Returns:

state vector \(\Phi(X_{nsteps})\) after integration

Return type:

numpy.array()