anchor.diff:        SH = SHAPE(M)
         IF (ABS(M(ROW,COL)).GE.TOL) THEN
           RATIO = M(ROW,COL) / M(COL,COL)
  -        DO 110 I = SCOL, SH(2)
  +        DO I = SCOL, SH(2)
           IF (I.EQ.COL) THEN
  -!
  -!  By definition this location must be zero.  Instead of doing the
  -!  calculation, just set the value to zero.
  -!
  -          M(ROW,I) = 0D0
  +          M(ROW,I) = 0D0 ! Zero by definition.  
  +                         ! don't calculate, just zero.
           ELSE
             M(ROW,I) = M(ROW,I) - RATIO * M(COL,I)
           ENDIF
  -  110  CONTINUE
  +       ENDDO
         ELSE
           M(ROW,COL) = 0D0
         END IF
  -  999 RETURN
     END SUBROUTINE MGEROW
   !------------------------------------------------------------------------------+
   ! MBSUB: Matrix Back SUBstitution
  
positive.code:     INTEGER(KIND=4),INTENT(IN) :: COL
      INTEGER(KIND=4)            :: SH(2),J
  !
        SH = SHAPE(M)
        DO 110 J = COL+1, SH(1)
    110 CALL MGEROW(M,J,COL,COL+1)
    999 RETURN
    END SUBROUTINE MGECOL
  !------------------------------------------------------------------------------+
  ! MGJEROW: Matrix Gauss-Jordan Eliminate ROW
  
negative.code:     INTEGER(KIND=4)            :: SH(2),J
  !
      SH = SHAPE(M)
      DO J = 1, SH(1)
        IF (J.NE.COL) CALL MGEROW(M,J,COL,1)
      ENDDO
    END SUBROUTINE MGJECOL
  !------------------------------------------------------------------------------+
  ! MGECOL: Matrix Gauss Eliminate COLumn
  