shrewd.m

Computes SHREWD weights for a given two- or three-angle spherical grid. See the paper by Eden and Levitt for details on now the algorithm works: http://dx.doi.org/10.1006/jmre.1998.1427

Syntax

weights=shrewd(alphas,betas,gammas,max_rank,max_error)

Arguments

alphas - alpha Euler angles (ZYZ active) of the

              grid, in radians, set to all-zeros for
              two-angle grids

      betas - beta Euler angles (ZYZ active) of the
              grid, in radians

     gammas - gamma Euler angles (ZYZ active) of the
              grid, in radians

   max_rank - maximum spherical rank to take into consi-
              deration when minimizing residuals

  max_error - maximum residual absolute error per spheri-
              cal function

Outputs

    weights - a vector of grid weights for each
              [alpha beta gamma] point supplied.

ilya.kuprov@weizmann.ac.il

Outputs

weights - a vector of grid weights for each

              [alpha beta gamma] point supplied.

See also

Integration grids, Appendix I: powder grids

Version 2.6, authors: Ilya Kuprov