1 // Copyright (C) 2007-2008 CEA/DEN, EDF R&D
3 // This library is free software; you can redistribute it and/or
4 // modify it under the terms of the GNU Lesser General Public
5 // License as published by the Free Software Foundation; either
6 // version 2.1 of the License.
8 // This library is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 // Lesser General Public License for more details.
13 // You should have received a copy of the GNU Lesser General Public
14 // License along with this library; if not, write to the Free Software
15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 #include "Interpolation3D.hxx"
20 #include "Interpolation3D.txx"
21 #include "MeshTestToolkit.txx"
23 #include "VectorUtils.hxx"
25 #include "MEDMEM_Mesh.hxx"
26 #include "MEDNormalizedUnstructuredMesh.hxx"
31 using namespace MEDMEM;
32 using namespace MED_EN;
36 * Test program which takes two meshes and calculates their intersection matrix.
38 * USAGE : PerfTest mesh1 mesh2
39 * where mesh1 and mesh2 are the names of two meshes located in
40 * the files mesh1.med, mesh2.med in {$MED_ROOT_DIR}/share/salome/resources/med/
47 * \brief Specialization of MeshTestToolkit for the purposes of performance testing.
50 class PerfTestToolkit : public MeshTestToolkit<3,3>
56 * Calculates the intersection matrix for two meshes.
57 * Outputs the names of the meshes intersected, the number of elements in each mesh,
58 * the number of matrix elements and the number of non-zero matrix elements, etc.
59 * These values help to determine how well the filtering algorithm is working.
61 * @param mesh1path the path to the file containing the source mesh, relative to {$MED_ROOT_DIR}/share/salome/resources/med/
62 * @param mesh1 the name of the source mesh
63 * @param mesh2path the path to the file containing the target mesh, relative to {$MED_ROOT_DIR}/share/salome/resources/med/
64 * @param mesh2 the name of the target mesh
65 * @param m intersection matrix in which to store the result of the intersection
67 void calcIntersectionMatrix(const char* mesh1path, const char* mesh1, const char* mesh2path, const char* mesh2, IntersectionMatrix& m)
69 const string dataBaseDir = getenv("MED_ROOT_DIR");
70 const string dataDir = dataBaseDir + "/share/salome/resources/med/";
72 LOG(1, std::endl << "=== -> intersecting src = " << mesh1 << ", target = " << mesh2 );
74 LOG(5, "Loading " << mesh1 << " from " << mesh1path);
75 const MESH sMesh(MED_DRIVER, dataDir+mesh1path, mesh1);
78 LOG(5, "Loading " << mesh2 << " from " << mesh2path);
79 const MESH tMesh(MED_DRIVER, dataDir+mesh2path, mesh2);
81 MEDNormalizedUnstructuredMesh<3,3> sMesh_wrapper(&sMesh);
82 MEDNormalizedUnstructuredMesh<3,3> tMesh_wrapper(&tMesh);
84 Interpolation3D interpolator;
85 interpolator.interpolateMeshes(sMesh_wrapper, tMesh_wrapper,m,"P0P0");
87 std::pair<int, int> eff = countNumberOfMatrixEntries(m);
88 LOG(1, eff.first << " of " << numTargetElems * numSrcElems << " intersections calculated : ratio = "
89 << double(eff.first) / double(numTargetElems * numSrcElems));
90 LOG(1, eff.second << " non-zero elements of " << eff.first << " total : filter efficiency = "
91 << double(eff.second) / double(eff.first));
93 LOG(1, "Intersection calculation done. " << std::endl );
98 * Counts the number of elements in an intersection matrix, and the number of these which are non-zero.
100 * @param m the intersection matrix
101 * @return pair<int, int> containing as its first element the number of elements in m and as its second element the
102 * number these which are non-zero
104 std::pair<int,int> countNumberOfMatrixEntries(const IntersectionMatrix& m)
109 for(IntersectionMatrix::const_iterator iter = m.begin() ; iter != m.end() ; ++iter)
111 numElems += iter->size();
112 for(map<int, double>::const_iterator iter2 = iter->begin() ; iter2 != iter->end() ; ++iter2)
114 if(!INTERP_KERNEL::epsilonEqual(iter2->second, 0.0, VOL_PREC))
120 return std::make_pair(numElems, numNonZero);
127 * Main method of the program.
128 * Intersects the meshes and outputs some information about the calculation as well as the
129 * intersection matrix on std::cout.
131 * @param argc number of arguments given to the program (should be 3, the user giving 2 mesh names)
132 * @param argv vector to the arguments as strings.
134 int main(int argc, char** argv)
136 using INTERP_TEST::PerfTestToolkit;
141 const string mesh1 = argv[1];
142 const string mesh2 = argv[2];
144 const string mesh1path = mesh1 + ".med";
145 const string mesh2path = mesh2 + ".med";
147 IntersectionMatrix m;
149 PerfTestToolkit testTools;
151 testTools.calcIntersectionMatrix(mesh1path.c_str(), mesh1.c_str(), mesh2path.c_str(), mesh2.c_str(), m);
153 testTools.dumpIntersectionMatrix(m);