solv_spacc_sparse

MCEq.solvers.solv_spacc_sparse(nsteps, dX, rho_inv, spacc_int_m, spacc_dec_m, phi, grid_idcs)

Apple Accelerate (vecLib) implementation.

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