SHGLQ Interface

public interface SHGLQ

Calls

interface~~shglq~~CallsGraph interface~shglq SHGLQ proc~shglq_dp SHGLQ_DP interface~shglq->proc~shglq_dp dcos dcos proc~shglq_dp->dcos

Called by

interface~~shglq~~CalledByGraph interface~shglq SHGLQ proc~decompose DECOMPOSE proc~decompose->interface~shglq proc~decompose~2 DECOMPOSE proc~decompose~2->interface~shglq proc~run_gwi_sweep RUN_GWI_SWEEP proc~run_gwi_sweep->interface~shglq proc~run_part_a RUN_PART_A proc~run_part_a->proc~decompose proc~run_part_a~2 RUN_PART_A proc~run_part_a~2->proc~decompose~2 proc~run_part_b RUN_PART_B proc~run_part_b->proc~decompose proc~run_part_b~2 RUN_PART_B proc~run_part_b~2->proc~decompose~2 program~gwi_gwj_sweep GWI_GWJ_SWEEP program~gwi_gwj_sweep->proc~run_gwi_sweep program~vsh_decomposition VSH_DECOMPOSITION program~vsh_decomposition->proc~run_part_a program~vsh_decomposition->proc~run_part_b program~vsh_decomposition_tor VSH_DECOMPOSITION_TOR program~vsh_decomposition_tor->proc~run_part_a~2 program~vsh_decomposition_tor->proc~run_part_b~2

Module Procedures

private subroutine SHGLQ_DP(ZERO, W, LMAX)

Gauss-Legendre quadrature nodes and weights, exact for integrating any polynomial of degree using only points -- the natural quadrature for spherical harmonic transforms, since is a polynomial in times , and the transform integral separates exactly into this -quadrature times a uniform sum over (itself exact for finite Fourier content). Nodes are the zeros of , found by Newton iteration (converges to machine precision in under 10 steps); weights are .

Read more…

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(out) :: ZERO(LMAX+1)
real(kind=dp), intent(out) :: W(LMAX+1)
integer(kind=i4), intent(in) :: LMAX

Requested exactness is degree ;