access gauss-legendre quadrature weights from fortran

Issue #5 resolved
Former user created an issue

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)

  1. Nathanaël Schaeffer 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.

  2. Log in to comment