Introducing sources in ODE?
Barry Smith
bsmith at mcs.anl.gov
Tue Aug 18 21:52:00 CDT 2009
I'm sorry, we never thought to add this support to the "linear
case". You could use the nonlinear TS interface to solve your problem,
but that is a bit overkill since it is only linear in U.
What particular integrator are you interested in using? Perhaps you
could add it yourself once you understand the source code in TS. Or
perhaps someone else has time to try to add this support?
Me, I'd just try doing it with the nonlinear interface where the
entire Au(t) + f(t) is the nonlinear RHS you provide and A is the
Jacobian you provide.
Barry
On Aug 17, 2009, at 9:58 AM, zhifeng sheng wrote:
> Dear all
>
> I would like to use Time-steppers - ODE integrators - TS to solve my
> time domain problem...
>
> However, the equation has the form Au(t) + f(t) = u'(t) ... and by
> looking into the example code, I could not find a way to introduce
> the source term f(t)... Does anyone know how to do that? If you do,
> please help with this.
>
> Thanks and best regards
> Zhifeng
More information about the petsc-users
mailing list