# The Higher Education and Research forge

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

# SCM Repository

 author breuils Tue, 24 Jul 2018 00:03:59 +0000 (02:03 +0200) committer breuils Tue, 24 Jul 2018 00:03:59 +0000 (02:03 +0200)

index cefa409..797ece8 100644 (file)
@@ -28,9 +28,9 @@ namespace c3ga2{
/// \return a multivector corresponding to a point of DCGA
template<typename T>
c3ga2::Mvec<T> point(const T &x, const T &y, const T &z){
-       c3ga2::Mvec<T> c3ga2Point;
+          c3ga2::Mvec<T> c3ga2Point;

-       c3ga2::Mvec<T> mvCGA1;
+          c3ga2::Mvec<T> mvCGA1;
c3ga2::Mvec<T> mvCGA2;
mvCGA1[c3ga2::E01] = mvCGA2[c3ga2::E02] = 1.0;
mvCGA1[c3ga2::E1]  = mvCGA2[c3ga2::E4]  = x;
@@ -38,9 +38,9 @@ namespace c3ga2{
mvCGA1[c3ga2::E3]  = mvCGA2[c3ga2::E6]  = z;
mvCGA1[c3ga2::Ei1] = mvCGA2[c3ga2::Ei2] = 0.5*(x*x+y*y+z*z);

-       c3ga2Point = mvCGA1^mvCGA2;
+          c3ga2Point = mvCGA1^mvCGA2;

-       return c3ga2Point;
+          return c3ga2Point;
}

@@ -52,29 +52,29 @@ namespace c3ga2{
/// \return the multivector corresponding to the Txsquared
template<typename T>
c3ga2::Mvec<T> Tx2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E14] = -1.0;
+          c3ga2::Mvec<T> mv;
+          mv[c3ga2::E14] = -1.0;

-       return mv;
+          return mv;
}

/// \brief build Tysquared
/// \return the multivector corresponding to the Tysquared
template<typename T>
c3ga2::Mvec<T> Ty2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E25] = -1.0;
-
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E25] = -1.0;
+
+       return mv;
}

/// \brief build Tzsquared
/// \return the multivector corresponding to the Tzsquared
template<typename T>
c3ga2::Mvec<T> Tz2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E36] = -1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E36] = -1.0;
+       return mv;
}

@@ -82,10 +82,10 @@ namespace c3ga2{
/// \return the multivector corresponding to the operator Txy
template<typename T>
c3ga2::Mvec<T> Txy(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E15] = -0.5;
-       mv[c3ga2::E24] = -0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E15] = -0.5;
+       mv[c3ga2::E24] = -0.5;
+       return mv;
}

@@ -93,20 +93,20 @@ namespace c3ga2{
/// \return the multivector corresponding to the operator Tzx
template<typename T>
c3ga2::Mvec<T> Tzx(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E16] = -0.5;
-       mv[c3ga2::E34] = -0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E16] = -0.5;
+       mv[c3ga2::E34] = -0.5;
+       return mv;
}

/// \brief build the cross term operator Tyz
/// \return the multivector corresponding to the operator Tyz
template<typename T>
c3ga2::Mvec<T> Tyz(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E35] = -0.5;
-       mv[c3ga2::E26] = -0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E35] = -0.5;
+       mv[c3ga2::E26] = -0.5;
+       return mv;
}

@@ -114,10 +114,10 @@ namespace c3ga2{
/// \return the multivector corresponding to the operator Tx
template<typename T>
c3ga2::Mvec<T> Tx(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E1i2] = 0.5;
-       mv[c3ga2::Ei14] = 0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E1i2] = 0.5;
+       mv[c3ga2::Ei14] = 0.5;
+       return mv;
}

@@ -125,29 +125,29 @@ namespace c3ga2{
/// \return the multivector corresponding to the operator Ty
template<typename T>
c3ga2::Mvec<T> Ty(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E2i2] = 0.5;
-       mv[c3ga2::Ei15] = 0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E2i2] = 0.5;
+       mv[c3ga2::Ei15] = 0.5;
+       return mv;
}

/// \brief build the single term operator Tz
/// \return the multivector corresponding to the operator Tz
template<typename T>
c3ga2::Mvec<T> Tz(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E3i2] = 0.5;
-       mv[c3ga2::Ei16] = 0.5;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E3i2] = 0.5;
+       mv[c3ga2::Ei16] = 0.5;
+       return mv;
}

/// \brief build the unit term operator Tz
/// \return the multivector corresponding to the operator T1
template<typename T>
c3ga2::Mvec<T> T1(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::Ei1i2] = -1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::Ei1i2] = -1.0;
+       return mv;
}

@@ -155,49 +155,49 @@ namespace c3ga2{
/// \return the multivector corresponding to the operator Tt4
template<typename T>
c3ga2::Mvec<T> Tt4(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E0102] = -4.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E0102] = -4.0;
+       return mv;
}

/// \brief build the cyclide component extraction operator
/// \return the multivector corresponding to the operator Tt2
template<typename T>
c3ga2::Mvec<T> Tt2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::Ei102] = -1.0;
-       mv[c3ga2::E01i2] = -1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::Ei102] = -1.0;
+       mv[c3ga2::E01i2] = -1.0;
+       return mv;
}

/// \brief build the cyclide component extraction operator
/// \return the multivector corresponding to the operator Txt2
template<typename T>
c3ga2::Mvec<T> Txt2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E102] = 1.0;
-       mv[c3ga2::E014] = 1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E102] = 1.0;
+       mv[c3ga2::E014] = 1.0;
+       return mv;
}

/// \brief build the cyclide component extraction operator
/// \return the multivector corresponding to the operator Tyt2
template<typename T>
c3ga2::Mvec<T> Tyt2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E202] = 1.0;
-       mv[c3ga2::E015] = 1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E202] = 1.0;
+       mv[c3ga2::E015] = 1.0;
+       return mv;
}

/// \brief build the cyclide component extraction operator
/// \return the multivector corresponding to the operator Tzt2
template<typename T>
c3ga2::Mvec<T> Tzt2(){
-       c3ga2::Mvec<T> mv;
-       mv[c3ga2::E302] = 1.0;
-       mv[c3ga2::E016] = 1.0;
-       return mv;
+       c3ga2::Mvec<T> mv;
+       mv[c3ga2::E302] = 1.0;
+       mv[c3ga2::E016] = 1.0;
+       return mv;
}

@@ -205,7 +205,7 @@ namespace c3ga2{
/// with ax^2 + by^2 + cz^2 + dxy + exz + fyz + gxw + hyw + izw + jw^2 = 0
/// \brief construct the quadric whose coefficients are {a,b,c,d,e,f,g,h,i,j}
template<typename T>