On 07/09/2010 05:35 PM, DE wrote:
just a hint for developers think about the so called "dense output" methods e.g. you want to get solution on a grid of 100 points but in fact a stepper does not always need to do 100 steps to produce the desired figures a frequent case may be that a stepper makes say 42 steps to reach the end point maintaining given accuracy and produces the solution on the grid points by some special means (a variant of interpolation) there is an item on this topic in the book "numerical recipes 3rd ed." 17.2.2 Dense Output (p.915)
Thanks for this hint. I think just that you mentioned this we will include some steppers with dense output, they are in fact very nice. We have tought in an earlier stage of this library about including them. But we droped that idea, because there exist only a few steppers which provide a output function and this concept can not be applied to all steppers in general. Karsten