2023-03-08 Version 1.1.1 Updated spec examples and associated documentation to avoid zeros in the solution vector. 2005-02-22 Version 1.1.0 FD05 dependence changed to FD15. 2004-07-12 Version 1.0.0 Version numbering added. 2002-12-06 The test for convergence of iterative refinement changed to avoid any problem with comparisons of numbers held in registers. 2001-08-09 Use of MF41 changed to MF71 1998-07-21 Redundant variables ONE removed. IFLAG in MA48DD made of length 7. Temporary variable NP introduced to avoid PFORT warnings. Third argument in calls to MA48DD changed to LA-NE. Bug corrected in print statements at end of ME48BD and start of ME48CD. Resetting of IPTRD and IRN after call to ME50BD moved to immediately after the call so that they will be reset correctly in the event of an error return from ME50BD. 1994-12-12 Calls of MC13D and MC21A changed to MC13DD and MC21AD 1994-06-14 Minor bugs in printing corrected. Minor bug rank calculation corrected. 1993-12-10 Minor bugs in printing corrected. 1993-10-18 Minor bug in ME48BD when ICNTL(3)=3 corrected. 1993-09-24 Some IVDEP comments added for speed on the Cray, some minor bugs fixed, defaults changed to BLAS 3 with block size 32 and no amalgamation of blocks of the block triangular form. 1993-05-20 Original version of ME48.