I think I discovered an error in the src-file DirectSolve.f90 (Elmer Rev: 6509M). At least my compiler (ifort version 14.0.1) throws an error message when I have warning mesages on.
Code: Select all
! DirectSolve.f90:
SUBROUTINE ComplexBandSolver( A,x,b, Free_fact )
...
REAL(KIND=dp) :: x(*),b(*)
...
CALL SolveComplexBandLapack( N,1,BA,x,Subband,3*Subband+1 )
...
CALL SolveComplexSBandLapack( N,1,BA,x,Subband,Subband+1 )
...
END SUBROUTINE ComplexBandSolver
! SolveBand.f90:
SUBROUTINE SolveComplexBandLapack( N,M,A,X,Subband,Band )
...
COMPLEX(KIND=dp) :: A(Band,N),X(M,N)
...
END
! SolveSBand.f90:
SUBROUTINE SolveComplexSBandLapack( N,M,A,X,Subband,Band )
...
COMPLEX(KIND=dp) :: A(Band,N),X(N,M)
...
END
What do you think?
Best regards, Stefan