anchor.diff:         CALL ConstructEmptyDenseMatrix(MatC,MatB%columns,MatA%rows)
       END IF
   
  +    !! Setup Lapack
       M = MatA%rows
       N = MatB%columns
       K = MatA%columns
  
positive.code:     INTEGER, INTENT(IN) :: rows
      INTEGER, INTENT(IN) :: columns
  
      CALL DestructDenseMatrix(this)
      ALLOCATE(this%data(rows,columns))
  
      this%rows = rows
      this%columns = columns
  
    END SUBROUTINE ConstructEmptyDenseMatrix
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
negative.code:     LDC = M
  
      CALL DGEMM(TRANSA, TRANSB, M, N, K, ALPHA, MatA%data, LDA, MatB%data, &
           & LDB, BETA, MatC%data, LDC)
  
    END SUBROUTINE MultiplyDense
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    PURE SUBROUTINE IncrementDenseMatrix(MatA,MatB,alpha_in)
  