1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML3.2 EN">
<HTML>
<HEAD>
<META NAME="GENERATOR" CONTENT="DOCTEXT">
<TITLE>TSARKIMEXRegister</TITLE>
</HEAD>
<BODY BGCOLOR="FFFFFF">
<A NAME="TSARKIMEXRegister"><H1>TSARKIMEXRegister</H1></A>
register an ARK IMEX scheme by providing the entries in the Butcher tableau and optionally embedded approximations and interpolation
<H3><FONT COLOR="#CC3333">Synopsis</FONT></H3>
<PRE>
#include "petscts.h"
PetscErrorCode TSARKIMEXRegister(const TSARKIMEXType name,PetscInt order,PetscInt s,
const PetscReal At[],const PetscReal bt[],const PetscReal ct[],
const PetscReal A[],const PetscReal b[],const PetscReal c[],
PetscInt pinterp,const PetscReal binterpt[],const PetscReal binterp[])
</PRE>
Not Collective, but the same schemes should be registered on all processes on which they will be used
<P>
<H3><FONT COLOR="#CC3333">Input Parameters</FONT></H3>
<TABLE border="0" cellpadding="0" cellspacing="0">
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>name </B></TD><TD>- identifier for method
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>order </B></TD><TD>- approximation order of method
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>s </B></TD><TD>- number of stages, this is the dimension of the matrices below
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>At </B></TD><TD>- Butcher table of stage coefficients for stiff part (dimension s*s, row-major)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>bt </B></TD><TD>- Butcher table for completing the stiff part of the step (dimension s; <A HREF="../Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A> to use the last row of At)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>ct </B></TD><TD>- Abscissa of each stiff stage (dimension s, <A HREF="../Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A> to use row sums of At)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>A </B></TD><TD>- Non-stiff stage coefficients (dimension s*s, row-major)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>b </B></TD><TD>- Non-stiff step completion table (dimension s; <A HREF="../Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A> to use last row of At)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>c </B></TD><TD>- Non-stiff abscissa (dimension s; <A HREF="../Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A> to use row sums of A)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>pinterp </B></TD><TD>- Order of the interpolation scheme, equal to the number of columns of binterpt and binterp
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>binterpt </B></TD><TD>- Coefficients of the interpolation formula for the stiff part (dimension s*pinterp)
</TD></TR>
<TR><TD WIDTH=40></TD><TD ALIGN=LEFT VALIGN=TOP><B>binterp </B></TD><TD>- Coefficients of the interpolation formula for the non-stiff part (dimension s*pinterp; <A HREF="../Sys/PETSC_NULL.html#PETSC_NULL">PETSC_NULL</A> to reuse binterpt)
</TD></TR></TABLE>
<P>
<H3><FONT COLOR="#CC3333">Notes</FONT></H3>
Several ARK IMEX methods are provided, this function is only needed to create new methods.
<P>
<P>
<H3><FONT COLOR="#CC3333">Keywords</FONT></H3>
<A HREF="../TS/TS.html#TS">TS</A>, register
<BR>
<P>
<H3><FONT COLOR="#CC3333">See Also</FONT></H3>
<A HREF="../TS/TSARKIMEX.html#TSARKIMEX">TSARKIMEX</A>
<BR><P><B><P><B><FONT COLOR="#CC3333">Level:</FONT></B>advanced
<BR><FONT COLOR="#CC3333">Location:</FONT></B><A HREF="../../../src/ts/impls/arkimex/arkimex.c.html#TSARKIMEXRegister">src/ts/impls/arkimex/arkimex.c</A>
<BR><A HREF="./index.html">Index of all TS routines</A>
<BR><A HREF="../../index.html">Table of Contents for all manual pages</A>
<BR><A HREF="../singleindex.html">Index of all manual pages</A>
<P><H3><FONT COLOR="#CC3333">Examples</FONT></H3>
<A HREF="../../../src/ts/examples/tutorials/ex16.c.html">src/ts/examples/tutorials/ex16.c.html</A><BR>
</BODY></HTML>
|
