next up previous
Next: Using the inequ input Up: User specified geometry and Previous: User specified geometry and

Warm start capability

Many applications involve running an equilibrium code many times while slightly varying some parameter. It is appropriate in this case to use the warm-start capability to accelerate the numerical convergence. A warm-start can be achieved by setting NDOUB to a negative value (e.g. -1) in the initialization call,

	CALL jsoInit(MN, -1, ISYM)
In this case the input values of PSIBAR0, X and Z are used as ``first guess.'' There is no mesh doubling associated with warm-starts. The sizes of X and Z arrays are therefore NTHEP*NPSI and must be pre-allocated, with

NTHEP = (ISYM-1)*(2**MN + ISYM) + 1 - 2*ISYM


NPSI = 2**MN + 1.

Mon May 10 12:32:55 EDT 1999