Real
PiecewiseConstant::integral() const
{
const unsigned len = functionSize();
Real sum = 0;
unsigned offset = 0;
if (_direction == RIGHT)
offset = 1;
for (unsigned i = 0; i < len - 1; ++i)
sum += range(i + offset) * (domain(i + 1) - domain(i));
return _scale_factor * sum;
}
Real
PiecewiseConstant::value(Real t, const Point & p)
{
Real func_value(0);
Real x = t;
if (_has_axis)
{
x = p(_axis);
}
unsigned i = 1;
const unsigned len = functionSize();
const Real toler = 1e-14;
// endpoint cases
if ( (_direction == LEFT && x < (1+toler) * domain(0)) ||
(_direction == RIGHT && x < (1-toler) * domain(0)) )
{
func_value = range(0);
i = len;
}
if ( (_direction == LEFT && x > (1+toler) * domain(len-1)) ||
(_direction == RIGHT && x > (1-toler) * domain(len-1)) )
{
func_value = range(len-1);
i = len;
}
for ( ; i < len; ++i )
{
if ( _direction == LEFT &&
x < (1+toler) * domain(i) )
{
func_value = range(i-1);
break;
}
else if ( (_direction == RIGHT &&
x < (1-toler) * domain(i)) )
{
func_value = range(i);
break;
}
}
return _scale_factor * func_value;
}
Real
PiecewiseConstant::average() const
{
return integral() / (domain(functionSize() - 1) - domain(0));
}
