void DriftTransport::operator() ( const std::array<double,2> &x , std::array<double,2> &dxdt , const double /* t */ )
{
// FIXME: avoid temporal creation overhead
Array<double> e_field((std::size_t) 2); // temp wrap for e. field
double e_field_mod;
Point eval_point(x[0],x[1],0.0);
(*_d_f_grad)(e_field, eval_point);
e_field_mod = sqrt(e_field[0]*e_field[0] + e_field[1]*e_field[1]);
dxdt[0] = _sign*_mu.obtain_mobility(e_field_mod) * e_field[0];
dxdt[1] = _sign*_mu.obtain_mobility(e_field_mod) * e_field[1];
}
void CatmullRomCurveEvaluator::evaluateCurve(const std::vector<Point>& controlPoints,
std::vector<Point>& evaluatedPoints,
const float& animationLength,
const bool& beWrap) const
{
evaluatedPoints.clear();
vector<Point> controlPointsCopy;
if (beWrap)
{
Point start_p1 = Point((controlPoints.end() - 2)->x - animationLength,
(controlPoints.end() - 2)->y);
Point start_p2 = Point((controlPoints.end() - 1)->x - animationLength,
(controlPoints.end() - 1)->y);
Point end_p1 = Point((controlPoints.begin())->x + animationLength,
(controlPoints.begin())->y);
Point end_p2 = Point((controlPoints.begin() + 1)->x + animationLength,
(controlPoints.begin() + 1)->y);
controlPointsCopy.push_back(start_p1);
controlPointsCopy.push_back(start_p2);
controlPointsCopy.insert(controlPointsCopy.end(), controlPoints.begin(), controlPoints.end());
controlPointsCopy.push_back(end_p1);
controlPointsCopy.push_back(end_p2);
}
else
{
controlPointsCopy.push_back(Point(0, controlPoints.front().y));
controlPointsCopy.insert(controlPointsCopy.end(), controlPoints.begin(), controlPoints.end());
controlPointsCopy.push_back(Point(animationLength, controlPoints.back().y));
}
const Mat4d basis = Mat4d(
-1, 3, -3, 1,
2, -5, 4, -1,
-1, 0, 1, 0,
0, 2, 0, 0) / 2.0;
for (size_t cnt = 0; cnt + 3 < controlPointsCopy.size(); ++cnt)
{
const Vec4d param_x(controlPointsCopy[cnt].x, controlPointsCopy[cnt + 1].x,
controlPointsCopy[cnt + 2].x, controlPointsCopy[cnt + 3].x);
const Vec4d param_y(controlPointsCopy[cnt].y, controlPointsCopy[cnt + 1].y,
controlPointsCopy[cnt + 2].y, controlPointsCopy[cnt + 3].y);
for (int i = 0; i < SEGMENT; ++i)
{
const double t = i / (double) SEGMENT;
const Vec4d param_time(t*t*t, t*t, t, 1);
Point eval_point(param_time * basis * param_x, param_time*basis*param_y);
// avoid wave curve occurs
if (evaluatedPoints.empty() || eval_point.x > evaluatedPoints.back().x)
{
evaluatedPoints.push_back(eval_point);
}
}
}
if (!beWrap)
{
// avoid slope interpolate at the end
if (controlPoints.back().x > evaluatedPoints.back().x)
{
evaluatedPoints.push_back(controlPoints.back());
}
evaluatedPoints.push_back(Point(0, controlPoints.front().y));
evaluatedPoints.push_back(Point(animationLength, controlPoints.back().y));
}
}
