Appendix E. THOMAS ALGORITHM FOR MULTI-LAYERS OF VG MODELThomas algorithm solves a tridiagonal matrix equation: into a form: where ai's, bi's, ci's, and fi's are give constants.bi's and yi'$s can be calculated by comparing corresponding elements of the above two equations, and xi's can be expressed: This equation can be solved by assuming: and obtaining ai's, bi's, and yi's incrementally: Then xi's can be obtained by a back substitution:
Last modified: Oct 15, 1999 VG Model / Samuel Lee / VADOSE.NET |