Page 190 - Basic Structured Grid Generation

Variational methods and adaptive grid generation 179

Fig. 6.7 Area-orthogonality.

2. Var1.ao.f
3. Var2.ao.f
4. SOR.ao.f
These programs solve the Euler-Lagrange equations based on the Area-Orthogonality
functional, with matrices given by eqns (6.61). The weight function ϕ is again a
constant, so that S = 0. Examples of the resulting grids for domains of various
shapes are shown in Figs 6.4–6.7.
5. Length.conjugate.f
This solves the Euler-Lagrange equations based on the Length-functional (6.40). Thus
it solves the pair of Laplace’s eqns (6.39). The numerical scheme used here is the
Conjugate Gradient method.

6.6.2 Subdirectory: Book/one.d.gds

This contains one file relevant to this chapter, as discussed in Section 6.3.3, namely

1. curve.SOR.f