Type Mismatch in DirectSolve.f90

Clearly defined bug reports and their fixes
Post Reply
sarz
Posts: 20
Joined: 06 Dec 2013, 16:53
Antispam: Yes

Type Mismatch in DirectSolve.f90

Post by sarz »

Hi there,

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
In the subroutine ComplexBandSolver x is of Real Type. Later the subroutines SolveComplexBandLapack and SolveComplexSBandLapack (placed in SolveBand.f90 and SolveSBand.f90) are called but expect x to be of Type Complex.
What do you think?

Best regards, Stefan
Post Reply