NAG AD Library
d03pp_a1w_f (dim1_parab_remesh_fd_a1w)

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object, as created by x10aa_a1w_f.

2: npde – Integer Input

3: m – Integer Input

4: ts – Type (nagad_a1w_w_rtype) Input/Output

5: tout – Type (nagad_a1w_w_rtype) Input

6: pdedef – Subroutine External Procedure

The specification of pdedef is:

Fortran Interface copy

Subroutine pdedef ( ad_handle, npde, t, x, u, ux, nv, v, vdot, p, q, r, ires, iuser, ruser) Integer, Intent (In) :: npde, nv Integer, Intent (Inout) :: ires, iuser(*) Type (nagad_a1w_w_rtype), Intent (In) :: t, x, u(npde), ux(npde), v(nv), vdot(nv) Type (nagad_a1w_w_rtype), Intent (Inout) :: ruser(*) Type (nagad_a1w_w_rtype), Intent (Out) :: p(npde,npde), q(npde), r(npde) Type (c_ptr), Intent (Inout) :: ad_handle

C++ Header Interface copy

#include <nagad.h> void pdedef ( void *&ad_handle, const Integer &npde, const nagad_a1w_w_rtype &t, const nagad_a1w_w_rtype &x, const nagad_a1w_w_rtype u[], const nagad_a1w_w_rtype ux[], const Integer &nv, const nagad_a1w_w_rtype v[], const nagad_a1w_w_rtype vdot[], nagad_a1w_w_rtype p[], nagad_a1w_w_rtype q[], nagad_a1w_w_rtype r[], Integer &ires, Integer iuser[], nagad_a1w_w_rtype ruser[])

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object.

2: npde – Integer Input

3: t – Type (nagad_a1w_w_rtype) Input

4: x – Type (nagad_a1w_w_rtype) Input

5: u – Type (nagad_a1w_w_rtype) array Input

6: ux – Type (nagad_a1w_w_rtype) array Input

7: nv – Integer Input

8: v – Type (nagad_a1w_w_rtype) array Input

9: vdot – Type (nagad_a1w_w_rtype) array Input

10: p – Type (nagad_a1w_w_rtype) array Output

11: q – Type (nagad_a1w_w_rtype) array Output

12: r – Type (nagad_a1w_w_rtype) array Output

13: ires – Integer Input/Output

14: iuser($*$) – Integer array User Workspace

15: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

7: bndary – Subroutine External Procedure

The specification of bndary is:

Fortran Interface copy

Subroutine bndary ( ad_handle, npde, t, u, ux, nv, v, vdot, ibnd, beta, gamma, ires, iuser, ruser) Integer, Intent (In) :: npde, nv, ibnd Integer, Intent (Inout) :: ires, iuser(*) Type (nagad_a1w_w_rtype), Intent (In) :: t, u(npde), ux(npde), v(nv), vdot(nv) Type (nagad_a1w_w_rtype), Intent (Inout) :: ruser(*) Type (nagad_a1w_w_rtype), Intent (Out) :: beta(npde), gamma(npde) Type (c_ptr), Intent (Inout) :: ad_handle

C++ Header Interface copy

#include <nagad.h> void bndary ( void *&ad_handle, const Integer &npde, const nagad_a1w_w_rtype &t, const nagad_a1w_w_rtype u[], const nagad_a1w_w_rtype ux[], const Integer &nv, const nagad_a1w_w_rtype v[], const nagad_a1w_w_rtype vdot[], const Integer &ibnd, nagad_a1w_w_rtype beta[], nagad_a1w_w_rtype gamma[], Integer &ires, Integer iuser[], nagad_a1w_w_rtype ruser[])

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object.

2: npde – Integer Input

3: t – Type (nagad_a1w_w_rtype) Input

4: u – Type (nagad_a1w_w_rtype) array Input

5: ux – Type (nagad_a1w_w_rtype) array Input

6: nv – Integer Input

7: v – Type (nagad_a1w_w_rtype) array Input

8: vdot – Type (nagad_a1w_w_rtype) array Input

9: ibnd – Integer Input

10: beta – Type (nagad_a1w_w_rtype) array Output

11: gamma – Type (nagad_a1w_w_rtype) array Output

12: ires – Integer Input/Output

13: iuser($*$) – Integer array User Workspace

14: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

8: uvinit – Subroutine External Procedure

The specification of uvinit is:

Fortran Interface copy

Subroutine uvinit ( ad_handle, npde, npts, nxi, x, xi, u, nv, v, iuser, ruser) Integer, Intent (In) :: npde, npts, nxi, nv Integer, Intent (Inout) :: iuser(*) Type (nagad_a1w_w_rtype), Intent (In) :: x(npts), xi(nxi) Type (nagad_a1w_w_rtype), Intent (Inout) :: ruser(*) Type (nagad_a1w_w_rtype), Intent (Out) :: u(npde,npts), v(nv) Type (c_ptr), Intent (Inout) :: ad_handle

C++ Header Interface copy

#include <nagad.h> void uvinit ( void *&ad_handle, const Integer &npde, const Integer &npts, const Integer &nxi, const nagad_a1w_w_rtype x[], const nagad_a1w_w_rtype xi[], nagad_a1w_w_rtype u[], const Integer &nv, nagad_a1w_w_rtype v[], Integer iuser[], nagad_a1w_w_rtype ruser[])

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object.

2: npde – Integer Input

3: npts – Integer Input

4: nxi – Integer Input

5: x – Type (nagad_a1w_w_rtype) array Input

6: xi – Type (nagad_a1w_w_rtype) array Input

7: u – Type (nagad_a1w_w_rtype) array Output

8: nv – Integer Input

9: v – Type (nagad_a1w_w_rtype) array Output

10: iuser($*$) – Integer array User Workspace

11: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

9: u(neqn) – Type (nagad_a1w_w_rtype) array Input/Output

10: npts – Integer Input

11: x(npts) – Type (nagad_a1w_w_rtype) array Input/Output

12: nv – Integer Input

13: odedef – Subroutine External Procedure

The specification of odedef is:

Fortran Interface copy

Subroutine odedef ( ad_handle, npde, t, nv, v, vdot, nxi, xi, ucp, ucpx, rcp, ucpt, ucptx, f, ires, iuser, ruser) Integer, Intent (In) :: npde, nv, nxi Integer, Intent (Inout) :: ires, iuser(*) Type (nagad_a1w_w_rtype), Intent (In) :: t, v(nv), vdot(nv), xi(nxi), ucp(npde,nxi), ucpx(npde,nxi), rcp(npde,nxi), ucpt(npde,nxi), ucptx(npde,nxi) Type (nagad_a1w_w_rtype), Intent (Inout) :: ruser(*) Type (nagad_a1w_w_rtype), Intent (Out) :: f(nv) Type (c_ptr), Intent (Inout) :: ad_handle

C++ Header Interface copy

#include <nagad.h> void odedef ( void *&ad_handle, const Integer &npde, const nagad_a1w_w_rtype &t, const Integer &nv, const nagad_a1w_w_rtype v[], const nagad_a1w_w_rtype vdot[], const Integer &nxi, const nagad_a1w_w_rtype xi[], const nagad_a1w_w_rtype ucp[], const nagad_a1w_w_rtype ucpx[], const nagad_a1w_w_rtype rcp[], const nagad_a1w_w_rtype ucpt[], const nagad_a1w_w_rtype ucptx[], nagad_a1w_w_rtype f[], Integer &ires, Integer iuser[], nagad_a1w_w_rtype ruser[])

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object.

2: npde – Integer Input

3: t – Type (nagad_a1w_w_rtype) Input

4: nv – Integer Input

5: v – Type (nagad_a1w_w_rtype) array Input

6: vdot – Type (nagad_a1w_w_rtype) array Input

7: nxi – Integer Input

8: xi – Type (nagad_a1w_w_rtype) array Input

9: ucp – Type (nagad_a1w_w_rtype) array Input

10: ucpx – Type (nagad_a1w_w_rtype) array Input

11: rcp – Type (nagad_a1w_w_rtype) array Input

12: ucpt – Type (nagad_a1w_w_rtype) array Input

13: ucptx – Type (nagad_a1w_w_rtype) array Input

14: f – Type (nagad_a1w_w_rtype) array Output

15: ires – Integer Input/Output

16: iuser($*$) – Integer array User Workspace

17: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

14: nxi – Integer Input

15: xi(nxi) – Type (nagad_a1w_w_rtype) array Input

16: neqn – Integer Input

17: rtol($*$) – Type (nagad_a1w_w_rtype) array Input

18: atol($*$) – Type (nagad_a1w_w_rtype) array Input

19: itol – Integer Input

20: norm – character Input

21: laopt – character Input

22: algopt($30$) – Type (nagad_a1w_w_rtype) array Input

23: remesh – logical Input

24: nxfix – Integer Input

25: xfix(nxfix) – Type (nagad_a1w_w_rtype) array Input

26: nrmesh – Integer Input

27: dxmesh – Type (nagad_a1w_w_rtype) Input

28: trmesh – Type (nagad_a1w_w_rtype) Input

29: ipminf – Integer Input

30: xratio – Type (nagad_a1w_w_rtype) Input

31: con – Type (nagad_a1w_w_rtype) Input

32: monitf – Subroutine External Procedure

The specification of monitf is:

Fortran Interface copy

Subroutine monitf ( ad_handle, t, npts, npde, x, u, r, fmon, iuser, ruser) Integer, Intent (In) :: npts, npde Integer, Intent (Inout) :: iuser(*) Type (nagad_a1w_w_rtype), Intent (In) :: t, x(npts), u(npde,npts), r(npde,npts) Type (nagad_a1w_w_rtype), Intent (Inout) :: ruser(*) Type (nagad_a1w_w_rtype), Intent (Out) :: fmon(npts) Type (c_ptr), Intent (Inout) :: ad_handle

C++ Header Interface copy

#include <nagad.h> void monitf ( void *&ad_handle, const nagad_a1w_w_rtype &t, const Integer &npts, const Integer &npde, const nagad_a1w_w_rtype x[], const nagad_a1w_w_rtype u[], const nagad_a1w_w_rtype r[], nagad_a1w_w_rtype fmon[], Integer iuser[], nagad_a1w_w_rtype ruser[])

1: ad_handle – Type (c_ptr) Input/Output

On entry: a handle to the AD configuration data object.

2: t – Type (nagad_a1w_w_rtype) Input

3: npts – Integer Input

4: npde – Integer Input

5: x – Type (nagad_a1w_w_rtype) array Input

6: u – Type (nagad_a1w_w_rtype) array Input

7: r – Type (nagad_a1w_w_rtype) array Input

8: fmon – Type (nagad_a1w_w_rtype) array Output

9: iuser($*$) – Integer array User Workspace

10: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

33: rsave(lrsave) – Type (nagad_a1w_w_rtype) array Communication Array

34: lrsave – Integer Input

35: isave(lisave) – Integer array Communication Array

36: lisave – Integer Input

37: itask – Integer Input

38: itrace – Integer Input

39: ind – Integer Input/Output

40: cwsav($10$) – character array Communication Array

If $\mathbf{ind}=0$, cwsav, lwsav, iwsav and rwsav need not be set on entry.

41: lwsav($100$) – logical array Communication Array

If $\mathbf{ind}=0$, cwsav, lwsav, iwsav and rwsav need not be set on entry.

42: iwsav($505$) – Integer array Communication Array

If $\mathbf{ind}=0$, cwsav, lwsav, iwsav and rwsav need not be set on entry.

43: rwsav($1100$) – Type (nagad_a1w_w_rtype) array Communication Array

If $\mathbf{ind}=0$, cwsav, lwsav, iwsav and rwsav need not be set on entry.

44: iuser($*$) – Integer array User Workspace

User workspace.

45: ruser($*$) – Type (nagad_a1w_w_rtype) array User Workspace

User workspace.

46: ifail – Integer Input/Output

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

10.1 Adjoint mode (a1w)

10.2 Tangent mode (t1w)

10.3 Passive mode (p0w)