2013-04-09 Jonathan Hogg Version 1.1.0. Fixed access to undefined values bugs found using nagfor. 2004-07-12 Version 1.0.0. Version numbering added. 2001-03-21 SAVE removed (parameters stored in first 6 locations of IW, so length of workspace increased). Also changed intrinsic DBLE to REAL 2000-05-09 Comments corrected. 1999-04-28 Error in flop count for static condensation variables cted. In MC53K/KD RINFO(2) = RINFO(2) + ISTAT*DBLE(JSTOP-JSTRT+1)**2 changed to NVAR = JSTOP-JSTRT+1 DO 45 J = 1,ISTAT RINFO(2) + 2*DBLE(NVAR-J)**2 45 CONTINUE 1998-12-16 IF (NELT.LE.1) THEN corrected to IF (NELT.LT.1) THEN 1998-12-16 Initialised INFO and RINFO so that they are not undefined in the event of an error message. 1998-08-14 In MC53H/HD changed bound on priority function from MAXPRT = -WEIGHT(1)*N*NSUP - WEIGHT(3)*N to MAXPRT = -WEIGHT(1)*MXGAIN - WEIGHT(3)*MXNADJ where MXGAIN holds maximum number of variables an element in the subdomain (not including the guard) can bring into the front and MXNADJ holds the maximum number of neighbours an element in the subdomain has (the guard is included as a possible neighbour) Change made to avoid integer overflow (N*NSUP can become large). 1998-06-01 Code modified to cope with the case of each element in the subdomain having NO variables in common with any other element in the subdomain (apart from the guard). In this case, we retain the original ordering. 1998-04-16 To avoid integer overflow RINFO(2) = RINFO(2) + DBLE(ISTAT* (JSTOP-JSTRT+1)**2) changed to RINFO(2) = RINFO(2) + ISTAT*DBLE(JSTOP-JSTRT+1)**2 Other operations with RINFO changed similarly. 1997-12-18 Bug found. Instead of LIWORK = 2*NZ + NELTP1 + 3 + NSUP + MAX(NELT,NSUP) we should have LIWORK = 2*NZ + NELTP1 + 3 + NSUP + MAX(NELTP1,NSUP) 1996-04-23 Changed NEED = USED + 2*NSUP + 6*NELTP1 + 2 to NEED = USED + 2*NSUP + 7*NELTP1 + 2. Change to spec. needed (increases upper bound on workspace ... term 9*NELT should be 10*NELT. 1995-10-08 code modified to include revised version of MC44. 1995-09-08 code modified to include extra information parameter INFO(8) which, on exit, holds number of variable indices ignored because they were out-of-range. 1995-07-03 code modified to include extra information parameter INFO(7) which, on exit, holds number of variables in the subdomain. 1995-03-10 Original version.