2012-01-03  Jonathan Hogg
Version 1.5.0 Default value of cntl%small changed to 1d-20.
	(Follows similar changes to symmetric codes.)

2010-12-22
Version 1.4.0 Default value of cntl%small changed to epsilon(1.0_wp).

2008-10-09
Version 1.3.0 Check made that entries do not become infinite (added error -14) 

2008-04-23
Version 1.2.0
Modification to diagonal pivoting if threshold u <= 0.0 and static = 0.0
(added error -13).

2008-02-26
Version 1.1.0
Slight modification to rook pivoting to reduce the number of columns searched.

2007-08-21
Version 1.0.0 Original Date