New version of library for vector version of Fortran implementation.

CMakeLists.txt

@@ -2,7 +2,7 @@
 project(wignerSymbols)
 set (wignerSymbols_VERSION_MAJOR 0)
 set (wignerSymbols_VERSION_MINOR 1)
-set (wignerSymbols_VERSION_RELEASE 2)
+set (wignerSymbols_VERSION_RELEASE 3)
 
 # CMake config
 cmake_minimum_required(VERSION 2.8)

README.md

@@ -4,7 +4,7 @@ Wigner Symbols
 
 A C++ ensemble of functions to compute the Wigner 3j- and 6j- symbols. It reimplements the algorithm designed
 by Schulten and Gordon in C++, but also contains the original Fortran implementation. 
-It can act either as a complelte C++ replacement to the original Fortran implementation, 
+It can act either as a complete C++ replacement to the original Fortran implementation, 
 or a C++ interface to them. See the API docs for details. 
 It can either compute an array of Wigner 3j or 6j symbols, or a single
 coefficient. It also computes the Clebsch-Gordan coefficients.

src/wignerSymbols-cpp.cpp

@@ -13,7 +13,7 @@ std::vector<double> wigner3j(double l2, double l3,
  double m1, double m2, double m3)
 {
  // We compute the numeric limits of double precision. 
- double huge = sqrt(std::numeric_limits<double>::max()/40.0);
+ double huge = sqrt(std::numeric_limits<double>::max()/20.0);
  double srhuge = sqrt(huge);
  double tiny = std::numeric_limits<double>::min();
  double srtiny = sqrt(tiny);

tests/testWigner.cpp

@@ -31,6 +31,8 @@ double timingWignerSymbolsF(double l2, double l3, double m1, double m2, double m
 
 int main(int argc, char* argv[])
 {
+ // We evalaute ThreeJ((1,0),(2,0),(1,0))
+ std::cout << WignerSymbols::wigner3j(1,2,1,0,0,0) << std::endl;
  // We time the computation for different sizes.
  int N = 500;
  double l2(500.0),m1(-10.0),m2(60.0),m3(-50.0);