anchor.diff:      LDC = M
   
       CALL DGEMM(TRANSA, TRANSB, M, N, K, ALPHA, MatA%data, LDA, MatB%data, &
  -          & LDB, BETA, MatC%data, LDC)
  -
  -    ! MatC%data = MATMUL(MatA%data,MatB%data)
  +         & LDB, BETA, MatC%data, LDC)
  +
     END SUBROUTINE MultiplyDense
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     PURE SUBROUTINE IncrementDenseMatrix(MatA,MatB,alpha_in)
  
positive.code:        CALL ConstructEmptyDenseMatrix(MatC,MatB%columns,MatA%rows)
      END IF
  
      M = MatA%rows
      N = MatB%columns
      K = MatA%columns
  
negative.code:     INTEGER :: total_counter
      TYPE(Triplet_t) :: temporary
  
      CALL ConstructEmptyDenseMatrix(dense_matrix,sparse_matrix%columns, &
           & sparse_matrix%rows)
  
      !! Loop over elements.
      dense_matrix%data = 0
  