The Higher Education and Research forge

Home My Page Projects Code Snippets Project Openings Garamon
Summary Activity SCM

SCM Repository

authorbreuils <stephane.breuils@u-pem.fr>
Mon, 23 Jul 2018 22:30:50 +0000 (00:30 +0200)
committerbreuils <stephane.breuils@u-pem.fr>
Mon, 23 Jul 2018 22:30:50 +0000 (00:30 +0200)
plugin/st3gaTools.hpp

index fd572e2..44f84bc 100644 (file)
@@ -27,7 +27,7 @@ namespace st3ga{
     template<typename T> 
     st3ga::Mvec<T> sigma_i(const unsigned int i){
         st3ga::Mvec<T> mv;
-        mv[1+1<<i]  = 1.0; // compute gamma 0 \wedge gamma i
+        mv[1+1<<i]  = -1.0; // compute gamma 0 \wedge gamma i
         return mv;
     }
 
@@ -36,7 +36,7 @@ namespace st3ga{
     template<typename T>
     st3ga::Mvec<T> sigma_1(){
         st3ga::Mvec<T> mv;
-        mv[st3ga::E01] = 1.0;
+        mv[st3ga::E01] = -1.0;
         return mv;
     }
 
@@ -45,15 +45,15 @@ namespace st3ga{
     template<typename T>
     st3ga::Mvec<T> sigma_2(){
         st3ga::Mvec<T> mv;
-        mv[st3ga::E02] = 1.0;
+        mv[st3ga::E02] = -1.0;
         return mv;
     }
 
     /// \brief build the bivectors basis used, third timelike bivector 
     template<typename T>
-    st3ga::Mvec<T> sigma3(){
+    st3ga::Mvec<T> sigma_3(){
         st3ga::Mvec<T> mv;
-        mv[st3ga::E03] = 1.0;
+        mv[st3ga::E03] = -1.0;
         return mv;
     }
 
@@ -67,26 +67,6 @@ namespace st3ga{
     }
 
 
-    // Check commutativity of the new basis, see Equation (9.5) of Section 9.1 of the reference
-    template<typename T>
-    bool isSigma_i_jCommutative(const T &epsilon = 1.0e-7){
-
-        //  compare the sign of sigma_i_1 wedge sigma_k_2 with sigma_k_2 wedge sigma_i_1
-        for(unsigned int i=1;i<3;++i){ // first particle index
-            for(unsigned int k=1;k<3;++k){ // second particle index
-                if(i!=k){
-                    st3ga::Mvec<T> sigma1 = sigma_j_i<T>(1,i); // first particle blade
-                    st3ga::Mvec<T> sigma2 = sigma_j_i<T>(2,k); // second particle blade
-
-                    st3ga::Mvec<T> resultDirect = sigma1 ^ sigma2;
-                    st3ga::Mvec<T> resultPermuted = sigma2 ^ sigma1;
-                    if(std::abs(coefficient(resultDirect, i, k) - coefficient(resultPermuted, i, k)) > epsilon) 
-                        return false;
-                }
-            }
-        }   
-        return true;
-    }