anchor.diff:         ENDIF
       ENDIF
   
  +    WRITE(*,*)' SpAMM SpAMM SpAMM SpAMM SpAMM SpAMM '
       CALL Multiply(P,P,P2)             ! The only multiplication is a square
  +    WRITE(*,*)' SpAMM SpAMM SpAMM SpAMM SpAMM SpAMM '
   
       dP2=MATMUL(dP,dP)
       dP2b=Quad2Dense(P2)
  
positive.code:     trc=Zero
      DO I=1,SIZE(A,1)
         trc=trc+A(I,I)
         !       WRITE(*,*)I,' TRACE ',A(I,I)
      ENDDO
    END FUNCTION Trace_dns
  
  
negative.code:     REAL(DOUBLE) :: Norm, CR1, CR2, TrP, TrP2,    TrdP,TrdP2
      INTEGER      :: I,  II,J
  
      INTEGER, PARAMETER:: NN=1024
      REAL(DOUBLE),DIMENSION(NN,NN) :: dP,dP2,dTmp1, dTmp2,dP2b
      !-------------------------------------------------------------------------------
  
      dP=Quad2Dense(P)
  
      IF (I.EQ.1) THEN
  