Finite Element Techniques for Fluid Flow

86 THE FINITE ELEMENT TECHNIQUE

The reader can verify by comparing this listing with that for the

full matrix that we are carrying out the same types of operations but

taking into account the symmetry and banded character of the

matrix and the different location of the column elements, which are

now stored in a 'skew' form. The vector C is simply a working space.

A check statement like IF(N.LT.I) G0T0 30 stops the program

from working outside the storage space for the real matrix, i.e. in

the 'dummy' area [see equation (c) of Example 2.3].

(b) SECOND CASE

For this case we have the classical eigenvalue problem*

tf*U = XJi*U (2.55)

where JT an ...