@@ -25,18 +25,23 @@ Look at the `libsharp-code/libsharp/sharp_cxx.h` file and add if necessary
sharp_geom_info* get_geom_info() { return ginfo;}
sharp_alm_info* get_alm_info() { return ainfo; }
* OpenBLAS (not Darwin):
* OpenBLAS (see note Darwin below):
If the variable `CBLAS` is defined then the code of laguerreTransform.cc use the BLAS matrix multiplication highly optimized. See installation in the Readme.md using
**WARNING** on Darwin the Makefile use the `Accelerate Framework` features which are a bit more effective than OpenBLAS.
**WARNING**: on Darwin the Makefile use the `Accelerate Framework` features which are a bit more effective than OpenBLAS.
# Compilation
> edit Makefile and adapat to local platform :
1) adapt to the type of Machine MacOSX (Darwin) vs Linux
2) adapt the location of the libsharp library (SHARPDIR)
3) on Linux adapt the location of BLAS/OpenBLAS
> make
4) then the make depends on the plateform and if you have BLAS
```
> make MACH=Linux BLAS=1 : on Linux with BLAS ON
> make : (default) on Darwin and use Accelerate Framework
```
> make check
runs a computation of the Nodes & Weights for N=1024
and compare to the xxx.txt.TEST files to produce a
...
...
@@ -68,14 +73,15 @@ If the variable `CBLAS` is defined then the code of laguerreTransform.cc use the
* ./Objs/lagsht_testsuite -n 1024
gives 2 files `lagNodes-1024-Func.txt` and `lagWeights-1024-Func.txt` with the Nodes & Weights computed
* ./Objs/lagsht_testsuite -t 3 -l 1024 -n 128
perform a Laguerre SH transform with the default geometry (Gauss) using Lmax = 1024 and Nmax = 128. The program generates a set of flmn coefficients, then performs a Synthesis followed by an Analysis operations. Finally it computes the maximal absolute and relative errors comparing the original flmn complex values to the results of the Analysis. The timing of the different processes is shown in a sorted tree.
perform a Laguerre SH transform with the default geometry (Gauss) using Lmax = 1024 and Nmax = 128. The program generates a set of flmn coefficients, then performs a Synthesis followed by an Analysis operations. Finally it computes the maximal absolute and relative errors comparing the original flmn complex values to the results of the Analysis. The timing of the different processes is shown in a sorted tree (here on Darwin Mac Os X with Accelerate Framework).
````
./Objs/lagsht_testsuite -t 3 -l 1024 -n 128
Max Memory size: 8589 MBytes
___________ MultiSphericalLaguerreTransform TEST _____________
LaguerreTransform start....
Sum roots = 16640 -> diff with theory: -4.9738e-14
Sum weights = 0.657058 -> diff with theory: 1.58752e-13
Sum roots = 16640 -> diff with theory: -5.68434e-14
Sum weights = 0.657058 -> diff with theory: 1.79737e-13
