diff --git a/types/types_mat.F90 b/types/types_mat.F90 index fb45d702d69797fc38b7ac458f672f964abae5e8..c627d136ee98c5f5b736e8cd5ea4eb02214916a4 100644 --- a/types/types_mat.F90 +++ b/types/types_mat.F90 @@ -43,7 +43,7 @@ MODULE m_types_mat REAL,ALLOCATABLE:: work(:) INTEGER,allocatable::ipiv(:) - IF ((mat%l_real.NE.vec%l_real).OR.(mat%matsize1.NE.mat%matsize2).OR.(mat%matsize1.NE.vec%matsize1)) & + IF ((mat%l_real.NEQV.vec%l_real).OR.(mat%matsize1.NE.mat%matsize2).OR.(mat%matsize1.NE.vec%matsize1)) & CALL judft_error("Invalid matices in t_mat_lproblem") IF (mat%l_real) THEN IF (ALL(ABS(mat%data_r-TRANSPOSE(mat%data_r))<1E-8)) THEN