Roll back an integration by mistake

untabify

Merge from V5_1_3_BR branch (07/12/09)

*** empty log message ***

Integration of Trio non regression test.

Integration of Trio modifications.

Protection of acos against values out of [-1.,1.]

Bug on 2 tests division by 0.

Suppression of bug on empty mesh.

Suppression of bug on empty mesh.

Activate floating point exceptions.

*** empty log message ***

*** empty log message ***

Polys management non conservative interpolators.

Polys management in MEDMEM.

Suppression of useless attributes.

*** empty log message ***

Dealing polyhedra for 3D P0P0 non conservative interpolation.

Dealing polyhedra for 3D P0P0 non conservative interpolation.

*** empty log message ***

non conservative Interp 3D

removing bug. MEDMEM test coming soon.

3D non conservative INTERP

Factorization IntersectorPlanar Intersector3D coming soon.

Adding 3D non conservative interp.

-pedantic -Wall compilation

-pedantic -Wall compilation

-pedantic -Wall compilation

-pedantic -Wall compilation

*** empty log message ***

*** empty log message ***

Adding non conservative interpolators for 2D and 3DSurf.
Suppression of precision degradation for P1P0Bary.

Adding non conservative interpolators for 2D and 3DSurf.

correct previous integration (Porting to Python 2.6)

Porting to Python 2.6 - add coding page specification for Python scripts

in test3DInterpP1P0Bary_1(), increase compare tolerance up to 1e-3

use tolerance at double comparison

0020525: [CEA] test MEDMEM fails

-      if ( srcMesh.getTypeOfElement( i ) != NORM_TETRA4 )
+      if ( srcMesh.getTypeOfElement( OTT<ConnType,numPol>::indFC( i )) != NORM_TETRA4 )

0020525: [CEA] test MEDMEM fails

-      if ( meshS.getTypeOfElement( i ) != NORM_TRI3 )
+      if ( meshS.getTypeOfElement( OTT<ConnType,numPol>::indFC( i )) != NORM_TRI3 )

windows port

windows port

windows port

windows port

windows port

Adding some missing headers in install.

-    std::pair< int, double[12] > subTetraNodes[24]; // a node of sub tetra and its coordinates
+    std::pair< int, std::vector<double> > subTetraNodes[24]; // a node of sub tetra and its coordinates

bug mantis 20457

fix memory leak in testInterpolationCC()

 const double * MEDCouplingNormalizedCartesianMesh<SPACEDIM,MESHDIM>::getCoordsAlongAxis(int axis) const
 {
-  return _mesh->getCoordsAt(axis)->getPointer();
+  return _mesh->getCoordsAt(axis)->getConstPointer();
 }

0020441: [CEA 349] P1P1 barycentric interpolators

0020441: [CEA 349] P1P1 barycentric interpolators
   {
+    // SPEC:
+    // "Limitation. For the P1P0 barycentric improvement only triangle source cells in 2D and
+    // tetrahedrons in 3D will be supported by interpolators. If a non
+    // triangle/tetrahedron source cell is detected an INTERP_KERNEL::Exception should be thrown."
+
+    // Check types of source elements here rather than in intersectCells() since a wrong type can be
+    // found late after a long time of calculation.
+
+    const unsigned long numSrcElems = meshS.getNumberOfElements();
+    for(unsigned long i = 0 ; i < numSrcElems ; ++i)
+      if ( meshS.getTypeOfElement( i ) != NORM_... )
+        throw INTERP_KERNEL::Exception("P1P0 barycentric algorithm works only with triangular source meshes");
   }

0020441: [CEA 349] P1P1 barycentric interpolators
  do not call intersectCells() if no source cells found
-        intersector->intersectCells(targetIdx,intersectElems,result);
+        if ( !intersectElems.empty() )
+          intersector->intersectCells(targetIdx,intersectElems,result);

0020441: [CEA 349] P1P1 barycentric interpolators

+    CPPUNIT_TEST( test3DInterpP1P1_1 );

remove extra lines

0020441: [CEA 349] P1P1 barycentric interpolators

+    else if(methC=="P1P1")
+      intersector=new PolyhedronIntersectorP1P1<MyMeshType,MatrixType>(targetMesh, srcMesh, getSplittingPolicy());

0020441: [CEA 349] P1P1 barycentric interpolators

+Intersector3DP1P1.hxx \
+Intersector3DP1P1.txx \
+PolyhedronIntersectorP1P1.hxx \
+PolyhedronIntersectorP1P1.txx \

new gcc and windows port

0020495: EDF : MEDMEM Documentation

0020440: [CEA 349] P1P0 barycentric interpolators

+    void test2DInterpP1P0Bary_1();
+    void test3DSurfInterpP1P0Bary_1();
+    void test3DInterpP1P0Bary_1();
   private:
+    MEDCouplingUMesh *build3DSourceMesh_2();
+    MEDCouplingUMesh *build3DTargetMesh_2();

0020440: [CEA 349] P1P0 barycentric interpolators

+    CPPUNIT_TEST( test_UnitTetraIntersectionBary_12 );

0020440: [CEA 349] P1P0 barycentric interpolators

+    double intersectWith(const QuadraticPolygon& other, double* barycenter) const;

0020440: [CEA 349] P1P0 barycentric interpolators

Fixing bugs

-    UnitTetraIntersectionBary();
+    UnitTetraIntersectionBary(bool isTetraInversed=false);

-    void init();
+    void init(bool isTetraInversed=false);

+    std::vector< std::vector< double > >  _polyNormals;

0020440: [CEA 349] P1P0 barycentric interpolators
0020441: [CEA 349] P1P1 barycentric interpolators

-    double intersectSourceCell(typename MyMeshType::MyConnType srcCell);
+    double intersectSourceCell(typename MyMeshType::MyConnType srcCell, double* baryCentre=0);
+
+    double intersectTetra(const double** tetraCorners);

+    void clearVolumesCache();

remove an extra line

0020440: [CEA 349] P1P0 barycentric interpolators

-    std::string methC(method);
+    std::string methC = InterpolationOptions::filterInterpolationMethod(method);
     if(methC=="P0P0")
       intersector=new PolyhedronIntersector<MyMeshType,MatrixType>(targetMesh, srcMesh, getSplittingPolicy());
     else if(methC=="P0P1")
       intersector=new PolyhedronIntersectorP0P1<MyMeshType,MatrixType>(targetMesh, srcMesh, getSplittingPolicy());
     else if(methC=="P1P0")
       intersector=new PolyhedronIntersectorP1P0<MyMeshType,MatrixType>(targetMesh, srcMesh, getSplittingPolicy());
+    else if(methC=="P1P0Bary")
+      intersector=new PolyhedronIntersectorP1P0Bary<MyMeshType,MatrixType>(targetMesh, srcMesh, getSplittingPolicy());

0020440: [CEA 349] P1P0 barycentric interpolators

+Intersector3DP1P0Bary.hxx \
+Intersector3DP1P0Bary.txx \
+PlanarIntersectorP1P0Bary.hxx \
+PlanarIntersectorP1P0Bary.txx \
+PolyhedronIntersectorP1P0Bary.hxx \
+PolyhedronIntersectorP1P0Bary.txx \

0020440: [CEA 349] P1P0 barycentric interpolators

-    std::string meth(method);
+    std::string meth = InterpolationOptions::filterInterpolationMethod(method);
     if(meth=="P0P0")
       {
       }
+    else if(meth=="P1P0Bary")
+    {
+      switch (InterpolationOptions::getIntersectionType())
+      {
+      case Triangulation:
+        intersector=new TriangulationIntersector<MyMeshType,MatrixType,PlanarIntersectorP1P0Bary>(myMeshT,myMeshS,_dim_caracteristic,
+      case Convex:
+        intersector=new ConvexIntersector<MyMeshType,MatrixType,PlanarIntersectorP1P0Bary>(myMeshT,myMeshS,_dim_caracteristic,
+      case Geometric2D:
+        intersector=new Geometric2DIntersector<MyMeshType,MatrixType,PlanarIntersectorP1P0Bary>(myMeshT, myMeshS, _dim_caracteristic,

0020440: [CEA 349] P1P0 barycentric interpolators
"Class InterpolationOptions should have an additional method : setP1P0BaryMethod(bool) that will set the _P1P0_bary_method attribute...An additional protected non virtual method treating string representation of interpolation method will be added:"

+    void setP1P0BaryMethod(bool isP1P0) { _P1P0_bary_method=isP1P0; }
+    bool getP1P0BaryMethod() const { return _P1P0_bary_method; }
+    std::string filterInterpolationMethod(const std::string& meth) const;

0020440: [CEA 349] P1P0 barycentric interpolators
"Concerning engine intersectors TriangleIntersector,Geometric2DIntersector and ConvexIntersector have to be extended with and additional method:"

+    double intersectGeoBary(const std::vector<double>& targetCell, bool targetCellQuadratic, const double *sourceCell, std::vector<double>& res);

Merge from BR_renumbering

EHPOC: L1.1.1: Structrured/structrured interpolation

Implementation of 1D-3D cartesian/cartesian interpolation

+MEDCouplingNormalizedCartesianMesh.hxx MEDCouplingNormalizedCartesianMesh.txx

EHPOC: L1.1.1: Structrured/structrured interpolation

Implementation of 1D-3D cartesian/cartesian interpolation

+    CPPUNIT_TEST( testInterpolationCC );

EHPOC: L1.1.1: Structrured/structrured interpolation

Implementation of 1D-3D cartesian/cartesian interpolation

+InterpolationCC.hxx                     \
+InterpolationCC.txx                     \

Fix problem of compilation: avoid linkage errors

correcting mistake in polyhedra connectivity documentation

*** empty log message ***

*** empty log message ***

*** empty log message ***

NO TABS PLEASE !!!

*** empty log message ***

Add one test.

delete tabs.

delete tabs.

P0->P1 parllel.

Moved the Remapper tests to MEDMEMCppTests

moved remapper tests to INTERP_KERNELTests

removed remapper tests

*** empty log message ***

Some new functionalities able to deal with dynamic types as classic one.

Some new functionalities : descending and reverse descending connectivities.

Added management of vector fields (multiple component)
Corrected the indentation

Issue 0020194: EDF 977 ALL: Get rid of warnings PACKAGE_VERSION already defined

AG : Debug for MEDMEM DataStructure.

Const forgotten.

*** empty log message ***

Debug

No tabs.

*** empty log message ***

*** empty log message ***

Suppression of some trace.