|
From: Luigi B. <lui...@gm...> - 2010-02-09 16:47:16
|
On Tue, 2010-02-09 at 09:40 +0000, na...@us... wrote:
> Log Message:
> -----------
> added MixedInterpolation
Question:
> + Real primitive(Real x) const {
> + if (x<*(this->xBegin2_))
> + return interpolation1_.primitive(x, true);
> + return interpolation2_.primitive(x, true);
> + }
shouldn't the second one be:
return interpolation1_.primitive(*xBegin2_, true) +
interpolation2_.primitive(x, true);
to account for the first piece of the integral?
Luigi
--
I am extraordinarily patient, provided I get my own way in the end.
-- Margaret Thatcher
|