shamrock.external.coala.L2_proj.L2proj#

shamrock.external.coala.L2_proj.L2proj(eps, nbins, kpol, 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 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
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

gij2D array (dim = (nbins,kpol+1)), type -> float: initial components of g on the polynomial basis