petsc-3.6.1 2015-08-06

TSSetForcingFunction

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

Synopsis

#include "petscts.h"  
PetscErrorCode  TSSetForcingFunction(TS ts,PetscErrorCode (*f)(TS,PetscReal,Vec,void*),void *ctx)

Logically Collective on TS

Input Parameters

	ts	- the TS context obtained from TSCreate()
	f	- 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

    func (TS ts,PetscReal t,Vec u,void *ctx);

	t	- current timestep
	u	- 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.

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

Keywords