shamrock.external.coala.limiter.gammafunction

shamrock.external.coala.limiter.gammafunction(eps, nbins, kpol, massgrid, massbins, gij)

Function to compute the limiter coefficient to ensure positivity of the numerical solution (Zhang and Shu 2010)

only for DG scheme k>0, piecewise polynomial approximation

Parameters

epsscalar, type -> float

minimum value for mass distribution approximation gij

nbinsscalar, type -> integer

number of dust bins

kpolscalar, type -> integer

degree of polynomials for approximation

massgrid1D array (dim = nbins+1), type -> float

grid of masses given borders value of mass bins

massbins1D array (dim = nbins), type -> float

arithmetic mean value of massgrid for each mass bins

gij2D array (dim = (nbins,kpol+1)), type -> float

components of g on the polynomial basis

Returns

tab_gamma1D array (dim = nbins), type -> float

limiter coefficient in each bin