2018-02-08 Tyrone Rees
Version 6.3.1
   * Fix error in section 2.5.3 of documentation	

2013-03-25  Jonathan Hogg
Version 6.3.0
   * Correct maximum value for control%u
   * Fix some compiler warnings concerning type conversions.

2012-01-03  Jonathan Hogg
Version 6.2.0
   * Store info%u in the case p=0
   * Change default value of control%small=1e-20 based on experience in
     HSL_MA86 and HSL_MA97 experiements.

2011-03-24  John Reid
Version 6.1.0
Determinant found using intrinsics exponent and fraction instead of log.
(For some matrices log was dominating runtime)

2011-01-19
Version 6.0.0
   * Complex versions added. In the Hermitian case, all diagonal
     entries are forced to be real.
   * If the number of active columns exceeds nbi, look for swaps.
   * Change default umin to 1.0 (ie disabled by default).
   * Use intrinsic sign to set static pivot value.
   * If 2x2 pivot tested but not accepted, use its best 1x1 pivot.
   * Many minor changes made, particularly to the comments.
   * Bugs corrected:
      1. detpiv and detscale stored for best candidate 2x2 pivot.
      2. if s is odd, sign of perm(s+1) ignored.
      3. corrections made re relaxed threshold pivoting
      4. Correct bug around bad 2x2 pivots.
      5. Fix bug with using cntl%umin rather than umin.

2009-08-18
Version 5.1.0
   * Looks for a 2x2 pivot before looking for a 1x1 pivot.
   * Argument info of ma64_factor given intent(inout).
   * Components of objects of type MA64_info intialized to zero.

2009-07-31
Version 5.0.0
   * cntl%umin and info%u added.
   * u allowed to be in the range [0.5,1.0].
   * Test on size of la relaxed.
   * Corrected bugs re info%usmall and the stability test for 2x2 pivots.
   * Changed to -1 the value of info%usmall returned if any pivots have been
     perturbed.

2008-10-22
Version 4.0.0
   * OpenMp features, the optional argument n_threads and the control component
     p_thresh added.
   * Rearrangement of the trailing n-p columns restored when executing on more
     than one thread or a single thread with p>nb.

2008-10-09
Version 3.1.0
   * Test for infinities added (flag value -13).
   * Flag value for mod(nb,nbi)/=1 changed to -12.
   * Pivot no longer accepted if absolute value of an entry in its
     inverse >= cntl%small**(-1).

2008-08-01
Version 3.0.0
   * Inner blocks added.
   * Trailing n-p columns not rearraged - instead updated by making temporary
     rearrangements of blocks of nb/2 columns.

2008-04-30
Version 2.0.0
   * Uses long integers internally and for the arguments la and ll so that n
     can be larger.
   * Code for determinant of 2x2 pivot scaled to avoid underflow.

2007-03-16
Version 1.0.0
Original Version