anchor.diff:       &             WRITE(6,*) '  (alpha spin)'
                 IF (NSPIN.EQ.2.AND.IJSPIN.EQ.2)
        &             WRITE(6,*) '  (beta spin)'
  -              CALL WRTMAT(WORK(KLFBLK),NI,NJ,NI,NJ)
  +              CALL WRTMAT(dbl_mb(KLFBLK),NI,NJ,NI,NJ)
               END IF
   *. Block has been constructed , transfer to -complete- 
   *. symmetry blocked Fock matrix
  
positive.code:               NK = NOBPTS(KGAS,KSM)
  *. blocks of one-electron integrals and one-electron density
                CALL GETD1(dbl_mb(KLDEN),JSM,JGAS,KSM,KGAS,1)
                CALL GETH1(WORK(KLINT),ISM,IGAS,KSM,KGAS)
                IF(NTEST.GE.1000) THEN
                  WRITE(6,*) 
       &          ' 1-e ints for ISM IGAS KGAS ',ISM,IGAS,KGAS
                  CALL WRTMAT(WORK(KLINT),NI,NK,NI,NK)
                  WRITE(6,*) 
       &          ' 1-e densi for ISM JGAS KGAS ',ISM,JGAS,KGAS
                  CALL WRTMAT(dbl_mb(KLDEN),NJ,NK,NJ,NK)
                END IF
  *. And then a matrix multiply( they are pretty much in fashion 
  *. these days )
                CALL MATML7(dbl_mb(KLFBLK),WORK(KLINT),dbl_mb(KLDEN),
       &                    NI,NJ,NI,NK,NJ,NK,ONE,ONE,2)
                 IF(NTEST.GE.1000) THEN
                   WRITE(6,*) ' Updated block '
  
negative.code:                 DO M = 1, NM
                    IIOFF = KLINT + (M-1)*NKL*NI
                    IDOFF = KLDEN + (M-1)*NKL*NJ
                    CALL MATML7(dbl_mb(KLFBLK),dbl_mb(IIOFF),
       &                   dbl_mb(IDOFF),NI,NJ,NKL,NI,NKL,NJ,ONE,ONE,1)
                  END DO
                END DO
                END DO
  