vb_spherical_grad

 spherical gradient dot
 [usage]
   gradF_Q = vb_spherical_grad( X, Q, N, Rmax, f_g )
 [input]
      X : coordinates of current (NP x 3)
      Q : normal vector of current (NP x 3)
      N : maximum order of spherical function
   Rmax : radial
    f_g : target current flag 'f'(F=internal) or 'g'(G=external)
 [output]
   grad : grad_x * Qx + grad_y * Qy + grad_z * Qz

 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function grad = vb_spherical_grad(X,Q,N,Rmax,f_g)
0002 % spherical gradient dot
0003 % [usage]
0004 %   gradF_Q = vb_spherical_grad( X, Q, N, Rmax, f_g )
0005 % [input]
0006 %      X : coordinates of current (NP x 3)
0007 %      Q : normal vector of current (NP x 3)
0008 %      N : maximum order of spherical function
0009 %   Rmax : radial
0010 %    f_g : target current flag 'f'(F=internal) or 'g'(G=external)
0011 % [output]
0012 %   grad : grad_x * Qx + grad_y * Qy + grad_z * Qz
0013 %
0014 % Copyright (C) 2011, ATR All Rights Reserved.
0015 % License : New BSD License(see VBMEG_LICENSE.txt)
0016 
0017 switch f_g
0018   case  'f'
0019     grad = vb_spherical_grad_F(X, Q, N, Rmax);
0020   case  'g'
0021     grad = vb_spherical_grad_G(X, Q, N, Rmax);
0022 end
0023 
0024 return;