next up previous contents
Next: Residual vector : Up: Matrix G and vector Previous: Matrix G and vector

Declarations:

In the above algorithm, G is an (it+2) tex2html_wrap_inline572 (it+2) matrix. Thus, at the top of your main program, you should declare space for G as a one-based, (Itmx+2) tex2html_wrap_inline572 (Itmx+2) dmatrix and also declare space for u as a one-based dvector of length Itmx+2. Finally, be good and free these variables at the end of main().


Tomas Arias
Mon Apr 2 13:24:52 EDT 2001