-
assigned issue to
- marked as enhancement
access gauss-legendre quadrature weights from fortran
Issue #5
resolved
Perhaps something like this:
{{{ /// fills the given array with the gaussian quadrature weights.
void shtns_gauss_wts_(double *wts) { for (int i=0; i<sht_data->nlat_2; i++) wts[i] = sht_data->wg[i]; } }}}
could be added to the fortran API in SHT.c?
Comments (3)
-
repo owner -
repo owner Last revision (f3160213e25b) provides the function shtns_gauss_wts() to both C and Fortran. Note that the gauss weights are stored after some renormalization, but these function return the "raw" gauss weights. Please confirm that this is ok for you, so I can close the issue.
-
repo owner - changed status to resolved
- Log in to comment