// Construct from components lineDivide::lineDivide(const curvedEdge& bc, const label n, const scalar xratio) : points_(n + 1), divisions_(n + 1), noPoints_(n) { scalar np(n); scalar lambda(0.0); if (xratio == 1.0) { scalar y(1.0/np); for (label i=0; i<=noPoints_; i++) { lambda = scalar(i)/np; points_[i] = bc.position(lambda); divisions_[i] = y*i; } } else { points_[0] = bc.position(0.0); divisions_[0] = 0.0; scalar xrpower = 1.0; for (label i=1; i<=noPoints_; i++) { lambda = (1.0 - pow(xratio, i))/(1.0 - pow(xratio, np)); points_[i] = bc.position(lambda); divisions_[i] = lambda; xrpower *= xratio; } } }
Foam::lineDivide::lineDivide ( const curvedEdge& cedge, const label ndiv, const scalar xratio ) : points_(ndiv + 1), divisions_(ndiv + 1) { divisions_[0] = 0.0; divisions_[ndiv] = 1.0; // calculate the spacing if (xratio == 1.0) { for (label i=1; i < ndiv; i++) { divisions_[i] = scalar(i)/ndiv; } } else { for (label i=1; i < ndiv; i++) { divisions_[i] = (1.0 - pow(xratio, i))/(1.0 - pow(xratio, ndiv)); } } // calculate the points for (label i=0; i <= ndiv; i++) { points_[i] = cedge.position(divisions_[i]); } }