shamrock.external.coala.L2_proj.L2proj_k0

shamrock.external.coala.L2_proj.L2proj_k0(eps, nbins, N0, m0, massgrid, massbins, Q, vecnodes, vecweights)

Function to compute the initial gij coefficient from the initial function g(m) = m*N0/m0*exp(-m/m0)

DG scheme k=0, piecewise constant approximation

Parameters

epsscalar, type -> float

minimum value for mass distribution approximation gij

nbinsscalar, type -> integer

number of dust bins

N0scalar, type -> float

initial total number density of grains

m0scalar, type -> float

initial mean mass of grains

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

Qscalar, type -> integer

number of points for Gauss-Legendre quadrature

vecnodes1D array (dim = Q), type -> float

nodes of the Legendre polynomials

vecweights1D array (dim = Q), type -> float

weights coefficients for the Gauss-Legendre polynomials

Returns

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

initial components of g on the polynomial basis