TSSetForcingFunction#

Provide a function that computes a forcing term for a ODE or PDE

Synopsis#

#include "petscts.h"  
PetscErrorCode TSSetForcingFunction(TS ts, TSForcingFunction func, void *ctx)

Logically Collective

Input Parameters#

ts - the TS context obtained from TSCreate()
func - routine for evaluating the forcing function
ctx - [optional] user-defined context for private data for the function evaluation routine (may be NULL)

Calling sequence of func#

PetscErrorCode func(TS ts, PetscReal t, Vec f, void *ctx)

ts - the integrator context
t - current timestep
f - output vector
ctx - [optional] user-defined function context

Notes#

This routine is useful for testing accuracy of time integration schemes when using the Method of Manufactured Solutions to create closed-form solutions with a non-physical forcing term. It allows you to use the Method of Manufactored Solution without directly editing the definition of the problem you are solving and hence possibly introducing bugs.

This replaces the ODE F(u,u_t,t) = 0 the TS is solving with F(u,u_t,t) - func(t) = 0

This forcing function does not depend on the solution to the equations, it can only depend on spatial location, time, and possibly parameters, the parameters can be passed in the ctx variable.

For low-dimensional problems solved in serial, such as small discrete systems, TSMonitorLGError() can be used to monitor the error history.