mkSolveAdapt - Specify linear algebra solvers with adaptive grids

Calling sequence:

mkSolveAdapt(solve_list,eq_list,unknown_list,old_unknown_list)

Parameters:

solve_list - List of keywords to denote different linear systems solver. One keyword for each equation in eq_list.

eq_list - List of equations to be solved.

unknown_list - The list of unknowns as discussed above.

old_unknown_list - The list of unknowns at the previous timestep, as discussed above.

Synopsis:

• Generates calls to different linear systems solvers as described in mkSolve.
• Generates a call to a subroutine performing grid refinement. The grid updates at each time step depending on values of ErrCrit (see mkAdapt).
• Introduces two parameters to indata_sample: adaptHmin specifying the smallest allowed size of a grid element and adaptTol setting threshold for grid refinement.

Examples:

> # set up different solvers to use with different equations:
> solve_list:=[iccg,iccg]:
> mkSolveAdapt(solve_list,eq_list,unknown_list,old_unknown_list);
dosolve
new rdpara.f
new indata_sample
solve