36#include "FixedVertexBasedDivisionRule.hpp"
38template <
unsigned SPACE_DIM>
41 if (fabs(norm_2(rDivisionVector) - 1.0) > DBL_EPSILON)
43 EXCEPTION(
"Input argument must be a unit vector");
45 mDivisionVector = rDivisionVector;
48template <
unsigned SPACE_DIM>
53 return mDivisionVector;
56template <
unsigned SPACE_DIM>
59 return mDivisionVector;
#define EXCEPTION(message)
#define EXPORT_TEMPLATE_CLASS_SAME_DIMS(CLASS)
FixedVertexBasedDivisionRule(c_vector< double, SPACE_DIM > &rDivisionVector)
virtual c_vector< double, SPACE_DIM > CalculateCellDivisionVector(CellPtr pParentCell, VertexBasedCellPopulation< SPACE_DIM > &rCellPopulation)
const c_vector< double, SPACE_DIM > & rGetDivisionVector() const