Difference between revisions of "Shrewd.m"

From Spinach Documentation Wiki
Jump to: navigation, search
(Created page with "Computes SHREWD weights for a given two- or three-angle spherical grid. ==Syntax== weights=shrewd(alphas,betas,gammas,max_rank,max_error) ==Description== See the paper...")
 
Line 1: Line 1:
Computes SHREWD weights for a given two- or three-angle spherical grid.
+
{{DISPLAYTITLE:function.m}} __NOTOC__
 +
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==
 
==Syntax==
  
 
     weights=shrewd(alphas,betas,gammas,max_rank,max_error)
 
     weights=shrewd(alphas,betas,gammas,max_rank,max_error)
 
==Description==
 
See the paper by Eden and Levitt for details on now the algorithm works: http://dx.doi.org/10.1006/jmre.1998.1427 - it is not entirely clear, in the magnetic resonance context, whether this procedure is even necessary.
 
  
 
==Arguments==
 
==Arguments==
  
       alphas - alpha Euler angles of the grid, in radians
+
       alphas - alpha Euler angles (ZYZ active) of the  
 +
              grid, in radians
 
   
 
   
       betas - beta Euler angles of the grid, in radians
+
       betas - beta Euler angles (ZYZ active) of the  
 +
              grid, in radians
 
   
 
   
       gammas - gamma Euler angles of the grid,in radians,
+
       gammas - gamma Euler angles (ZYZ active) of the
              set to all-zeros for two-angle grids
+
              grid,in radians, set to all-zeros for
 +
              two-angle grids
 
   
 
   
 
     max_rank - maximum spherical rank to take into consi-
 
     max_rank - maximum spherical rank to take into consi-
Line 23: Line 24:
 
               cal function
 
               cal function
  
==Returns==
+
==Outputs==
The output is a vector of grid weights for each [alpha beta gamma] point supplied.
 
 
==Examples==
 
See kernel/grids directory for a long list of one-, two- and three angle grids.
 
  
==Notes==
+
    weights - a vector of grid weights for each
for a given arrangement of angles, this is the most consistent weight selection procedure in the literature. Lebedev grids satisfy SHREWD condition by definition, all other grids may be rebalanced towards better description of low spherical ranks by the use of this procedure.
+
              [alpha beta gamma] point supplied.
  
 
==See also==
 
==See also==
[[get_hull.m]], [[grid_kron.m]], [[grid_test.m]], [[repulsion.m]]
+
[[Kernel_utilities#Integration_grids|Integration grids]], [[Appendix I: powder grids]]
  
  
''Version 1.10, authors: [[Ilya Kuprov]]''
+
''Version 2.6, authors: [[Ilya Kuprov]]''

Revision as of 20:15, 29 June 2021

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

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

     gammas - gamma Euler angles (ZYZ active) of the
              grid,in radians, set to all-zeros for
              two-angle grids

   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.

See also

Integration grids, Appendix I: powder grids


Version 2.6, authors: Ilya Kuprov