solv_numpy#

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

numpy implementation of forward-euler integration.

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)\)

Returns:

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

Return type:

numpy.array()