A profile solver in c for finite element equations

Abstract

This paper presents an efficient, pointer based profile solver with standard matrix indexing. Constrained equations Ax = b where x contains known and unknown values are solved and the full vectors x and b are obtained. Pseudo-code algorithms are formulated for a row oriented form of the LDL T factorization and implemented directly as a C code. The solver is implemented in C because of the close relation between two-dimensional arrays and pointers which makes it possible to write a clear and efficient code.