No commit message

No commit message

Geo/STensor3.h

@@ -245,6 +245,15 @@ class STensor3 {
     ithis.setMat(m);
     return ithis;
   }
+  STensor3 transpose () const
+  {
+    STensor3 ithis;
+    for (int i = 0; i < 3; i++)
+      for (int j = 0; j < 3; j++)
+        ithis(i,j) = (*this)(j,i);
+    return ithis;
+  }
+
   STensor3 operator + (const STensor3 &other) const
   {
     STensor3 res(*this);

Numeric/fullMatrix.h

@@ -84,12 +84,12 @@ class fullVector
       _data = new scalar[_r];
       _own_data = true;
       if(resetValue)
-        setAll(0.);
+        setAll(scalar());
       return true;
     }
     _r = r;
     if(resetValue)
-      setAll(0.);
+      setAll(scalar());
     return false;
   }
   void setAsProxy(const fullVector<scalar> &original, int r_start, int r)