2004-07-12 Version 1.0.0 Version numbering added. 2001-04 MF49 replaced by MC59 (to make code threadsafe) We also now exploit zeros within the front. This can reduce factor size (but can increase times because of extra searching). 2000-02 Changes needed because of changes to ME42. In particular, so that user's data is not altered by ME43B/BD, we take a copy of each row of A before it is used. The real workspace required must be increased by length of row with most entries. 1997-09 Changes made to the code because for equation entry to ME42, the symbolic factorization accurately predicts the max. frontsize and storage required by the factors. This is because for equation entry to the frontal solver, pivots are NEVER delayed (unless matrix is singular). For equation entry, the pivoting stategy is partial pivoting. This makes the control parameter CNTL(2) redundant. 1993-04-05 Original version of ME43.