Fixed prestore mapping.

[phys/su2clebsch.git] / su2clebsch.hpp

diff --git a/su2clebsch.hpp b/su2clebsch.hpp
index 689c0a1424fb1e7b94902b4749270081a4d5b563..955868c04f9871279d5c58308bf73466fac81a7b 100644 (file)
--- a/su2clebsch.hpp
+++ b/su2clebsch.hpp
@@ -1,11 +1,27 @@
 #ifndef SU2CLEBSCH_HPP
 #define SU2CLEBSCH_HPP
 
+using namespace std;
+
 namespace su2clebsch {
 
+  struct prestore {
+  private:
+    unsigned int max2j;
+    double *cgcdata;
+  public:
+    prestore(unsigned int _max2j);
+    ~prestore();
+
+    friend double cgc (int factor1_2j, int factor1_2m, int factor2_2j, int factor2_2m,
+                      int irrep_2j, int irrep_2m, const prestore& Prestore);
+  };
+  
   double cgc (int factor1_2j, int factor1_2m, int factor2_2j, int factor2_2m,
              int irrep_2j, int irrep_2m);
 
+  double cgc (int factor1_2j, int factor1_2m, int factor2_2j, int factor2_2m,
+             int irrep_2j, int irrep_2m, const prestore& Prestore);
 };
 
 #endif