MEDMEM Industrialization 2008

+    void reverseApply(double* destPt, const double* srcPt) const;

diff --git a/src/INTERP_KERNEL/TetraAffineTransform.cxx b/src/INTERP_KERNEL/TetraAffineTransform.cxx
index e018176a32709bb760d438ed5c4fd707f55e5687..6a57bf49a9918d6a15daa90ee93f64f5c3fffab3 100644 (file)
--- a/src/INTERP_KERNEL/TetraAffineTransform.cxx
+++ b/src/INTERP_KERNEL/TetraAffineTransform.cxx
@@ -35,6 +35,9 @@ namespace INTERP_KERNEL
             _linearTransform[3*j + i] = (pts[i+1])[j] - (pts[0])[j];
           }
       }
+
+    // remember _linearTransform for the reverse transformation
+    memcpy( _backLinearTransform, _linearTransform, 9*sizeof(double));
     
     calculateDeterminant();
     
@@ -127,6 +130,48 @@ namespace INTERP_KERNEL
       }
   }
 
+  /**
+   * Calculates the reverse transform of point srcPt and stores the result in destPt.
+   * If destPt == srcPt, then srcPt is overwritten safely.
+   *
+   * @param destPt  double[3] in which to store the transformed point
+   * @param srcPt   double[3] containing coordinates of points to be transformed
+   */
+  void TetraAffineTransform::reverseApply(double* destPt, const double* srcPt) const
+  {
+    double* dest = destPt;
+    
+    // are we self-allocating ?
+    const bool selfAllocation = (destPt == srcPt);
+    
+    if(selfAllocation)
+      {
+        // alloc temporary memory
+        dest = new double[3];
+       
+        LOG(6, "Info : Self-affectation in TetraAffineTransform::reverseApply");
+      }
+    
+    for(int i = 0 ; i < 3 ; ++i)
+      {
+        // matrix - vector multiplication
+        dest[i] = _backLinearTransform[3*i] * srcPt[0] + _backLinearTransform[3*i + 1] * srcPt[1] + _backLinearTransform[3*i + 2] * srcPt[2];
+
+        // translation
+        dest[i] -= _translation[i];
+      }
+
+    if(selfAllocation)
+      {
+        // copy result back to destPt
+        for(int i = 0 ; i < 3 ; ++i)
+          {
+            destPt[i] = dest[i];
+          }
+        delete[] dest;
+      }
+  }
+
   /**
    * Returns the determinant of the linear part A of the transform.
    *

diff --git a/src/INTERP_KERNEL/TetraAffineTransform.hxx b/src/INTERP_KERNEL/TetraAffineTransform.hxx
index 835bdbb2880a8b5cb274fbaa3a2a83716ac2e3d7..138e53875b7139d7a40a043a86062afadf71c5aa 100644 (file)
--- a/src/INTERP_KERNEL/TetraAffineTransform.hxx
+++ b/src/INTERP_KERNEL/TetraAffineTransform.hxx
@@ -21,6 +21,8 @@ namespace INTERP_KERNEL
 
     void apply(double* destPt, const double* srcPt) const;
 
+    void reverseApply(double* destPt, const double* srcPt) const;
+
     double determinant() const;
 
     void dump() const;
@@ -49,6 +51,8 @@ namespace INTERP_KERNEL
     /// The determinant of the matrix A is calculated at the construction of the object and cached.
     double _determinant;
 
+    /// 3x3 matrix AT is transposed A matrix used for y -> AT(y - b) transformation
+    double _backLinearTransform[9];
   };