The module HSL_MP62 uses the multiple front method to solve sets of symmetric positive-definite finite-element equations that have been divided into non-overlapping subdomains. The HSL routines MA62 and MA72 are used with MPI for message passing.
TThe coefficient matrix must be of the form
where the summation is over finite elements. The element matrix is nonzero only in those rows and columns which correspond to variables in the -th element. The right-hand side(s) may optionally be in the form
where is nonzero only in those rows which correspond to variables in element .
In the multiple front method, a frontal decomposition is performed on each subdomain separately. Thus, on each subdomain, and factors are computed. Once all possible eliminations have performed within a subdomain, there remain the interface variables, which are shared by more than one subdomain. If is the remaining frontal matrix for subdomain , and is the corresponding right-hand side matrix, then the remaining problem is
(1) |
where and . By treating each as an element matrix, the interface problem (1) is also solved by the frontal method. Once (1) has been solved, back-substitution on the subdomains completes the solution.
The element data and/or the matrix factors are optionally held in direct-access files.