2023-03-27
Version 3.6.0
When finding free unit numbers only search up to huge(0_short)-1.

2022-11-29
Version 3.5.1
Update the spec example to avoid zero output for portability.

2012-01-03
Version 3.5.0
Change default value of cntl%small=1e-20 following experience in symmetric
case.

2010-12-22
Version 3.4.0
   * Default value of cntl%small changed to epsilon(1.0_wp).
   * Also fix small bug with overwriting unit_error.

2009-10-02
Version 3.3.0
   * Some variables required for restarted were not saved when they needed to
     be.
   * Also bug fixed when solving transpose systems with scale present.
   * Tidied up some print statements.

2009-07-30
Version 3.2.1
Some variables undefined in restart (detected by NAG compiler with
-C=undefined). These have been fixed.

2008-10-30
Version 3.2.0
   * Rewritten code for doing depth-first search of tree (was not a bug fix but
     simplifies code).
   * Removed a couple of lines of redundant code from scale77.
   * Added error return -41 (error in Fortran close)

2008-10-08
Version 3.1.0
   * Removed variable itopd (was not used since integer info. on delayed pivots
     not held in a stack).
   * Removed scalar size_fa and array ip from keep.
   * Added error return -39 (IEEE infinities)

2008-08-05
Version 3.0.2 Line lengths corrected.

2008-05-28
Version 3.0.1 Bug in ma78_scale corrected.

2008-05-12
Version 3.0.0
   * Use of long integers modified to allow code to run on 64-bit architectures.
     In particular, the frontal matrix is addressed with a long integer (and so
     too is maxstore).
   * Bug fixed in restarting.

2008-04-28
Version 2.0.0
   * Scaling offered (with option to compute matrix norm).
   * In MA78_resid, anorm replaced by anorm_bnd, and if transpose system being
     solved, computes bound on norm of A^T.
   * Bug fixed to allow ma78_input_reals to be called after a previous
     factorization.
   * Also some other minor bug fixes.

2008-03-04
Version 1.8.0
Error -16 extended to include Fortran open unsuccessful for all elements of path

2008-02-12
Version 1.7.0
Removed all recursive subroutines (since a very large number of levels of
recursion can cause code to fail)

2008-01-30
Version 1.6.0 Bug fixed in constructing tree.

2008-01-14
Version 1.5.0
   * We now keep the max front from analyse and from factor (since if factor
     called without recalling analyse, want to use the lower bound from
     analyse).
   * Also fixed bug on where ip1 is set in ma78_factor,
   * and bug in size of fa when allocated at start of factorize.
   * And in MA78_allocate, need to use reals for computing lfa_long.

2001-12-12
Version 1.4.0 Minor bug in factorize phase, info%flag = 6 set incorrectly.

2007-12-10
Version 1.3.0
hsl_of01 updated. inactive parameter for of01_write used when writing factor
data. discard parameter used when reading stack data

2007-11-20
Version 1.2.0
Bug fixed in MA78_write_real and MA78_write_integer when switching from an
array to a file.

2007-11-06
Version 1.1.0

2007-08-21
Version 1.0.0
Initial public release

2007-07-30
arrays x, b, resid are no longer assumed shape because copying could add an
overhead (warning from lf95)