Start debugging 3D interpolation error on OCTA12 in target mesh

diff --git a/src/INTERP_KERNEL/SplitterTetra.hxx b/src/INTERP_KERNEL/SplitterTetra.hxx
index dc082002ceb5df22e54093cb177e8eef29add831..2f3c5f0b8ac0d4ce475b7d91869638655cae1905 100644 (file)
--- a/src/INTERP_KERNEL/SplitterTetra.hxx
+++ b/src/INTERP_KERNEL/SplitterTetra.hxx
@@ -363,7 +363,7 @@ namespace INTERP_KERNEL
     
     /// Smallest volume of the intersecting elements in the transformed space that will be returned as non-zero. 
     /// Since the scale is always the same in the transformed space (the target tetrahedron is unitary), this number is independent of the scale of the meshes.
-    static const double SPARSE_TRUNCATION_LIMIT=1.0e-14;
+    static const double SPARSE_TRUNCATION_LIMIT;
   };
 
   /**

diff --git a/src/INTERP_KERNEL/SplitterTetra.txx b/src/INTERP_KERNEL/SplitterTetra.txx
index 8364bb1b4a5782d9a547cf9e7cf820c0c47d7e31..6364c79d9e96ae2f5d4c7bd54a76514bd7348a52 100644 (file)
--- a/src/INTERP_KERNEL/SplitterTetra.txx
+++ b/src/INTERP_KERNEL/SplitterTetra.txx
@@ -38,6 +38,9 @@
 
 namespace INTERP_KERNEL
 {
+  template<class MyMeshType> 
+  const double SplitterTetra<MyMeshType>::SPARSE_TRUNCATION_LIMIT=1.0e-14;
+
   /*!
    * output is expected to be allocated with 24*sizeof(void*) in order to store the 24 tetras.
    * These tetras have to be deallocated.