MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Combining InterpolatingFunctions

  • To: mathgroup at
  • Subject: [mg107186] Re: Combining InterpolatingFunctions
  • From: Roland Franzius <roland.franzius at>
  • Date: Fri, 5 Feb 2010 03:20:53 -0500 (EST)
  • References: <hk8nqs$8ne$>

Simon Pearce schrieb:
> Hi MathGroup,
> I have two sets of InterpolatingFunctions coming from two separate
> NDSolve's. One of them is defined over the region [0,rc] and the other
> over the region [rc,2]. I would like Mathematica to automatically choose
> the correct one when I use a replacement rule. If I could tell it never
> to extrapolate this would be perfect, though I don't seem to be able to.
> I've tried using FunctionInterpolation, but in order to keep my error
> terms down I had to increase the InterpolationPoints to 1000, which
> increases the calculation time from approximately .5sec to 1.5sec.
> Can anyone suggest an efficient way of combining InterpolatingFunctions
> without re-interpolating them? Or turning the extrapolation off!

The interpolating functions have standard data structures, eg

a=NDSolve[{y'[x]== Cos[x],y[0]==1},y,{x,0,2pi} ]//First

b=NDSolve[{y'[x]== -Sin[x],y[2Pi]==(y/.a)[2Pi]},y,{x,2pi,4PI}]//First

a[[1,2,1]]  etc

The domains and value tables of a and b can be combined in an obvious 
way using Part to construct the joined List of domains and value tables.

The List, then, may be converted Applying InterpolatingFunction.


Roland Franzius

  • Prev by Date: Re: Bug? Analytical integration of cosines gets the sign wrong
  • Next by Date: Re: DeleteDuplicates is too slow?
  • Previous by thread: Re: Combining InterpolatingFunctions
  • Next by thread: Re: Combining InterpolatingFunctions