1 // Copyright (C) 2007-2014 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, or (at your option) any later version.
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 // Author : Anthony Geay (CEA/DEN)
21 #include "InterpolationOptions.hxx"
22 #include "InterpKernelException.hxx"
26 const double INTERP_KERNEL::InterpolationOptions::DFT_MEDIAN_PLANE=0.5;
28 const double INTERP_KERNEL::InterpolationOptions::DFT_SURF3D_ADJ_EPS=1.e-4;
30 const double INTERP_KERNEL::InterpolationOptions::DFT_MAX_DIST_3DSURF_INTERSECT=-1.;
32 const double INTERP_KERNEL::InterpolationOptions::DFT_MIN_DOT_BTW_3DSURF_INTERSECT=-1.;
34 const char INTERP_KERNEL::InterpolationOptions::PRECISION_STR[]="Precision";
36 const char INTERP_KERNEL::InterpolationOptions::MEDIANE_PLANE_STR[]="MedianPlane";
38 const char INTERP_KERNEL::InterpolationOptions::BOUNDING_BOX_ADJ_STR[]="BoundingBoxAdjustment";
40 const char INTERP_KERNEL::InterpolationOptions::BOUNDING_BOX_ADJ_ABS_STR[]="BoundingBoxAdjustmentAbs";
42 const char INTERP_KERNEL::InterpolationOptions::MAX_DISTANCE_3DSURF_INSECT_STR[]="MaxDistance3DSurfIntersect";
44 const char INTERP_KERNEL::InterpolationOptions::MIN_DOT_BTW_3DSURF_INSECT_STR[]="MinDotBetween3DSurfIntersect";
46 const char INTERP_KERNEL::InterpolationOptions::PRINT_LEV_STR[]="PrintLevel";
48 const char INTERP_KERNEL::InterpolationOptions::DO_ROTATE_STR[]="DoRotate";
50 const char INTERP_KERNEL::InterpolationOptions::ORIENTATION_STR[]="Orientation";
52 const char INTERP_KERNEL::InterpolationOptions::MEASURE_ABS_STR[]="MeasureAbs";
54 const char INTERP_KERNEL::InterpolationOptions::INTERSEC_TYPE_STR[]="IntersectionType";
56 const char INTERP_KERNEL::InterpolationOptions::SPLITTING_POLICY_STR[]="SplittingPolicy";
58 const char INTERP_KERNEL::InterpolationOptions::TRIANGULATION_INTERSECT2D_STR[]="Triangulation";
60 const char INTERP_KERNEL::InterpolationOptions::CONVEX_INTERSECT2D_STR[]="Convex";
62 const char INTERP_KERNEL::InterpolationOptions::GEOMETRIC_INTERSECT2D_STR[]="Geometric2D";
64 const char INTERP_KERNEL::InterpolationOptions::POINTLOCATOR_INTERSECT_STR[]="PointLocator";
66 const char INTERP_KERNEL::InterpolationOptions::BARYCENTRIC_INTERSECT_STR[]="Barycentric";
68 const char INTERP_KERNEL::InterpolationOptions::BARYCENTRICGEO2D_INTERSECT_STR[]="BarycentricGeo2D";
70 const char INTERP_KERNEL::InterpolationOptions::PLANAR_SPLIT_FACE_5_STR[]="PLANAR_FACE_5";
72 const char INTERP_KERNEL::InterpolationOptions::PLANAR_SPLIT_FACE_6_STR[]="PLANAR_FACE_6";
74 const char INTERP_KERNEL::InterpolationOptions::GENERAL_SPLIT_24_STR[]="GENERAL_24";
76 const char INTERP_KERNEL::InterpolationOptions::GENERAL_SPLIT_48_STR[]="GENERAL_48";
78 void INTERP_KERNEL::InterpolationOptions::init()
81 _intersection_type=Triangulation;
83 _median_plane=DFT_MEDIAN_PLANE;
85 _bounding_box_adjustment=DFT_SURF3D_ADJ_EPS;
86 _bounding_box_adjustment_abs=0.;
87 _max_distance_for_3Dsurf_intersect=DFT_MAX_DIST_3DSURF_INTERSECT;
88 _min_dot_btw_3Dsurf_intersect=DFT_MIN_DOT_BTW_3DSURF_INTERSECT;
91 _splitting_policy=PLANAR_FACE_5;
94 std::string INTERP_KERNEL::InterpolationOptions::getIntersectionTypeRepr() const
96 if(_intersection_type==INTERP_KERNEL::Triangulation)
97 return std::string(TRIANGULATION_INTERSECT2D_STR);
98 else if(_intersection_type==INTERP_KERNEL::Convex)
99 return std::string(CONVEX_INTERSECT2D_STR);
100 else if(_intersection_type==INTERP_KERNEL::Geometric2D)
101 return std::string(GEOMETRIC_INTERSECT2D_STR);
102 else if(_intersection_type==INTERP_KERNEL::PointLocator)
103 return std::string(POINTLOCATOR_INTERSECT_STR);
104 else if(_intersection_type==INTERP_KERNEL::Barycentric)
105 return std::string(BARYCENTRIC_INTERSECT_STR);
106 else if(_intersection_type==INTERP_KERNEL::BarycentricGeo2D)
107 return std::string(BARYCENTRICGEO2D_INTERSECT_STR);
109 return std::string("UNKNOWN_INTERSECT_TYPE");
112 bool INTERP_KERNEL::InterpolationOptions::setOptionDouble(const std::string& key, double value)
114 if(key==PRECISION_STR)
119 else if(key==MEDIANE_PLANE_STR)
121 setMedianPlane(value);
124 else if(key==BOUNDING_BOX_ADJ_STR)
126 setBoundingBoxAdjustment(value);
129 else if(key==BOUNDING_BOX_ADJ_ABS_STR)
131 setBoundingBoxAdjustmentAbs(value);
134 else if(key==MAX_DISTANCE_3DSURF_INSECT_STR)
136 setMaxDistance3DSurfIntersect(value);
139 else if(key==MIN_DOT_BTW_3DSURF_INSECT_STR)
141 setMinDotBtwPlane3DSurfIntersect(value);
148 bool INTERP_KERNEL::InterpolationOptions::setOptionInt(const std::string& key, int value)
150 if(key==PRINT_LEV_STR)
152 setPrintLevel(value);
155 else if(key==DO_ROTATE_STR)
157 setDoRotate(value != 0);
160 else if(key==ORIENTATION_STR)
162 setOrientation(value);
165 else if(key==MEASURE_ABS_STR)
167 bool valBool=(value!=0);
168 setMeasureAbsStatus(valBool);
175 bool INTERP_KERNEL::InterpolationOptions::setOptionString(const std::string& key, const std::string& value)
177 if(key==INTERSEC_TYPE_STR)
179 if(value==TRIANGULATION_INTERSECT2D_STR)
181 setIntersectionType(INTERP_KERNEL::Triangulation);
184 else if(value==CONVEX_INTERSECT2D_STR)
186 setIntersectionType(INTERP_KERNEL::Convex);
189 else if(value==GEOMETRIC_INTERSECT2D_STR)
191 setIntersectionType(INTERP_KERNEL::Geometric2D);
194 else if(value==POINTLOCATOR_INTERSECT_STR)
196 setIntersectionType(INTERP_KERNEL::PointLocator);
199 else if(value==BARYCENTRIC_INTERSECT_STR)
201 setIntersectionType(INTERP_KERNEL::Barycentric);
204 else if(value==BARYCENTRICGEO2D_INTERSECT_STR)
206 setIntersectionType(INTERP_KERNEL::BarycentricGeo2D);
210 else if(key==SPLITTING_POLICY_STR)
212 if(value==PLANAR_SPLIT_FACE_5_STR)
214 setSplittingPolicy(INTERP_KERNEL::PLANAR_FACE_5);
217 else if(value==PLANAR_SPLIT_FACE_6_STR)
219 setSplittingPolicy(INTERP_KERNEL::PLANAR_FACE_6);
222 else if(value==GENERAL_SPLIT_24_STR)
224 setSplittingPolicy(INTERP_KERNEL::GENERAL_24);
227 else if(value==GENERAL_SPLIT_48_STR)
229 setSplittingPolicy(INTERP_KERNEL::GENERAL_48);
238 std::string INTERP_KERNEL::InterpolationOptions::getSplittingPolicyRepr() const
240 if(_splitting_policy==INTERP_KERNEL::PLANAR_FACE_5)
241 return std::string(PLANAR_SPLIT_FACE_5_STR);
242 else if(_splitting_policy==INTERP_KERNEL::PLANAR_FACE_6)
243 return std::string(PLANAR_SPLIT_FACE_6_STR);
244 else if(_splitting_policy==INTERP_KERNEL::GENERAL_24)
245 return std::string(GENERAL_SPLIT_24_STR);
246 else if(_splitting_policy==INTERP_KERNEL::GENERAL_48)
247 return std::string(GENERAL_SPLIT_48_STR);
249 return std::string("UNKNOWN_SPLITTING_POLICY");
252 std::string INTERP_KERNEL::InterpolationOptions::filterInterpolationMethod(const std::string& meth) const
254 return std::string(meth);
257 bool INTERP_KERNEL::InterpolationOptions::setInterpolationOptions(long print_level,
258 std::string intersection_type,
262 double bounding_box_adjustment,
263 double bounding_box_adjustment_abs,
264 double max_distance_for_3Dsurf_intersect,
267 std::string splitting_policy)
269 _print_level=print_level;
270 _precision=precision;
271 _median_plane=median_plane;
272 _do_rotate=do_rotate;
273 _bounding_box_adjustment=bounding_box_adjustment;
274 _bounding_box_adjustment_abs=bounding_box_adjustment_abs;
275 _max_distance_for_3Dsurf_intersect=max_distance_for_3Dsurf_intersect;
276 _orientation=orientation;
277 _measure_abs=measure_abs;
278 return(setOptionString(INTERSEC_TYPE_STR,intersection_type) && setOptionString(SPLITTING_POLICY_STR,splitting_policy));
281 std::string INTERP_KERNEL::InterpolationOptions::printOptions() const
283 std::ostringstream oss; oss.precision(15); oss << "Interpolation Options ******" << std::endl;
284 oss << "Print level : " << _print_level << std::endl;
285 oss << "Intersection type : " << getIntersectionTypeRepr() << std::endl;
286 oss << "Precision : " << _precision << std::endl;
287 oss << "Median plane : " << _median_plane << std::endl;
288 oss << "Do Rotate status : " << std::boolalpha << _do_rotate << std::endl;
289 oss << "Bounding box adj : " << _bounding_box_adjustment << std::endl;
290 oss << "Bounding box adj abs : " << _bounding_box_adjustment_abs << std::endl;
291 oss << "Max distance for 3DSurf intersect : " << _max_distance_for_3Dsurf_intersect << std::endl;
292 oss << "Min dot between plane for 3DSurf intersect : " << _min_dot_btw_3Dsurf_intersect << std::endl;
293 oss << "Orientation : " << _orientation << std::endl;
294 oss << "Measure abs : " << _measure_abs << std::endl;
295 oss << "Splitting policy : " << getSplittingPolicyRepr() << std::endl;
296 oss << "****************************" << std::endl;
300 void INTERP_KERNEL::InterpolationOptions::CheckAndSplitInterpolationMethod(const std::string& method, std::string& srcMeth, std::string& trgMeth)
302 const int NB_OF_METH_MANAGED=4;
303 const char *METH_MANAGED[NB_OF_METH_MANAGED]={"P0P0","P0P1","P1P0","P1P1"};
305 for(int i=0;i<NB_OF_METH_MANAGED && !found;i++)
306 found=(method==METH_MANAGED[i]);
309 std::string msg("The interpolation method : \'"); msg+=method; msg+="\' not managed by INTERP_KERNEL interpolators ! Supported are \"P0P0\", \"P0P1\", \"P1P0\" and \"P1P1\".";
310 throw INTERP_KERNEL::Exception(msg.c_str());
312 srcMeth=method.substr(0,2);
313 trgMeth=method.substr(2);