slag_basis.m Checked in

Author: Boris Leistedt

src/main/c/flag_spherbessel.c

@@ -84,11 +84,11 @@ void flag_sbesselslag(double *sbesselslag, int ell, double *kvalues, int Nk, int
         if( j == 0 ){
           weight =  pow(2, 2) * gsl_sf_fact(ell + 2) * (n + 1.0) * (n + 2.0) / 2.0 ; // 
         }else{
-          weight = - 2 * (ell + 2 + j) * (n - j + 1) * weight / (j * (j + 2)); //  
+          weight = - 2 * (n - j + 1) * weight / ((j + 2)); //  
         }
         //sbesselslag[n * Nk + k] += pow((n+1)*(n+2), -0.5) * weight * mulk_coefs[j+2] ;
         printf(" %3.0e ",weight);
-        vals[j] = pow((n+1)*(n+2), -0.5) * weight * mulk_coefs[j+2];
+        vals[j] = pow((n+1)*(n+2), -0.5) * ( weight * mulk_coefs[j+2] );
       }
       result = 0.0, hold = 0.0, temp = 0.0, y = 0.0;
       for(j = 0; j <= n/2; j++){

src/main/matlab/slag_basis.m

@@ -0,0 +1,21 @@
+function f = slag_basis(N, nodes, tau)
+
+% slag_basis_mex - Compute the spherical Laguerre basis functions
+% up to order N on a grid of radii nodes (with rescaling factor tau)
+%
+% Default usage :
+%
+%   f = slag_basis(N, nodes, tau)
+%   plot(f)
+%
+% where N is the maximum order/mode of the spherical Laguerre functions,
+% nodes is the radial grig on which to compute the basis function,
+% tau is the rescaling factor.
+%
+% FLAG package to perform 3D Fourier-Laguerre Analysis
+% Copyright (C) 2012  Boris Leistedt & Jason McEwen
+% See LICENSE.txt for license details
+
+f = slag_basis_mex(N, nodes, tau);
+  
+end