matrix[11]=-p0[0]*matrix[8]-p0[1]*matrix[9]-p0[2]*matrix[10];
return true;
}
-
- template<int SPACEDIM>
- void ComputeTriangleHeight(const double *PA, const double *PB, const double *PC, double *res)
- {
- double AB = getDistanceBtw2Pts<SPACEDIM>(PA,PB);
- double BC = getDistanceBtw2Pts<SPACEDIM>(PB,PC);
- double CA = getDistanceBtw2Pts<SPACEDIM>(PC,PA);
- double perim( (AB+BC+CA)*0.5 );
- double num( 2*sqrt(perim*(perim-AB)*(perim-BC)*(perim-CA)) );
- res[0] = num/AB; res[1] = num/BC; res[2] = num/CA;
- }
}
#include "VolSurfUser.hxx"
#include "VolSurfFormulae.hxx"
#include "InterpolationUtils.hxx"
+#include "VectorUtils.hxx"
#include <algorithm>
return computeBarycenter<ConnType,numPolConn,1>(type,connec,lgth,coords,res);
throw INTERP_KERNEL::Exception("Invalid spaceDim specified for compute barycenter : must be 1, 2 or 3");
}
+
+ template<int SPACEDIM>
+ void ComputeTriangleHeight(const double *PA, const double *PB, const double *PC, double *res)
+ {
+ double AB = getDistanceBtw2Pts<SPACEDIM>(PA,PB);
+ double BC = getDistanceBtw2Pts<SPACEDIM>(PB,PC);
+ double CA = getDistanceBtw2Pts<SPACEDIM>(PC,PA);
+ double perim( (AB+BC+CA)*0.5 );
+ double num( 2*sqrt(perim*(perim-AB)*(perim-BC)*(perim-CA)) );
+ res[0] = num/AB; res[1] = num/BC; res[2] = num/CA;
+ }
}
#endif
#include "DiameterCalculator.hxx"
#include "OrientationInverter.hxx"
#include "InterpKernelAutoPtr.hxx"
-#include "VolSurfUser.cxx"
+#include "VolSurfUser.txx"
using namespace MEDCoupling;
this->checkAllocated();
constexpr char MSG1[] = "DataArrayInt::locateComponentId : single component array expected";
valToSearchIntoTuples->checkNbOfComps(1,MSG1); tupleIdHint->checkNbOfComps(1,MSG1);
- auto nbOfCompo( this->getNumberOfComponents() ),thisNbTuples( this->getNumberOfTuples() );
- auto nbOfTuples( valToSearchIntoTuples->getNumberOfTuples() );
+ std::size_t nbOfCompo( this->getNumberOfComponents() );
+ mcIdType thisNbTuples( this->getNumberOfTuples() );
+ mcIdType nbOfTuples( valToSearchIntoTuples->getNumberOfTuples() );
tupleIdHint->checkNbOfTuples(nbOfTuples,"Number of tuples of input arrays must be the same.");
const T *cPtr(this->begin()),*valSearchPt(valToSearchIntoTuples->begin());
const mcIdType *tHintPtr(tupleIdHint->begin());