The work array, RWORK, is declared COMPLEX, it should be REAL.
* ..
* .. Local allocatable arrays
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ REAL, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
REAL ZERO
* ..
* .. Local allocatable arrays
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ REAL, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
REAL ZERO
* ..
* .. Local allocatable arrays
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
DOUBLE PRECISION ZERO
* ..
* .. Local allocatable arrays
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
DOUBLE PRECISION ZERO