Skip to content

Commit

Permalink
Progress to resolving the bug with complex testers
Browse files Browse the repository at this point in the history
  • Loading branch information
Daniel Bielich committed Feb 9, 2023
1 parent 21059dc commit c65a7a7
Show file tree
Hide file tree
Showing 5 changed files with 599 additions and 298 deletions.
7 changes: 5 additions & 2 deletions SRC/zgedmdq.f90
Original file line number Diff line number Diff line change
Expand Up @@ -8,9 +8,9 @@ SUBROUTINE ZGEDMDQ( JOBS, JOBZ, JOBR, JOBQ, JOBT, JOBF, &
!.....
! USE PRECISION, ONLY: WP => DP ! double precision
! IMPLICIT NONE
use iso_fortran_env, only: real32
use iso_fortran_env, only: real64
IMPLICIT NONE
integer, parameter :: WP = real32
integer, parameter :: WP = real64

!.....
! Scalar arguments
Expand Down Expand Up @@ -616,6 +616,8 @@ SUBROUTINE ZGEDMDQ( JOBS, JOBZ, JOBR, JOBQ, JOBT, JOBF, &
ZWORK(1) = MLWORK
ZWORK(2) = OLWORK
WORK(1) = MLRWRK
!WRITE(*,*) 'ZWORK(1) = ', ZWORK(1)
!WRITE(*,*) 'ZWORK(2) = ', ZWORK(2)
RETURN
END IF
!.....
Expand Down Expand Up @@ -645,6 +647,7 @@ SUBROUTINE ZGEDMDQ( JOBS, JOBZ, JOBR, JOBQ, JOBT, JOBF, &
EIGS, Z, LDZ, RES, B, LDB, V, LDV, &
S, LDS, ZWORK(MINMN+1), LZWORK-MINMN, &
WORK, LWORK, IWORK, LIWORK, INFO1 )
!S, LDS, ZWORK(MINMN+1), LZWORK-MINMN, &
IF ( INFO1 == 2 .OR. INFO1 == 3 ) THEN
! Return with error code. See ZGEDMD for details.
INFO = INFO1
Expand Down
10 changes: 9 additions & 1 deletion TESTING/EIG/Makefile
Original file line number Diff line number Diff line change
Expand Up @@ -83,6 +83,8 @@ CEIGTST = cchkee.o \
csgt01.o cslect.o csyl01.o\
cstt21.o cstt22.o cunt01.o cunt03.o

CDMDEIGTST = cchkdmd.o

DZIGTST = dlafts.o dlahd2.o dlasum.o dlatb9.o dstech.o dstect.o \
dsvdch.o dsvdct.o dsxt1.o

Expand Down Expand Up @@ -129,13 +131,16 @@ all: single complex double complex16

.PHONY: single complex double complex16
single: xeigtsts xdmdeigtsts
complex: xeigtstc
complex: xeigtstc xdmdeigtstc
double: xeigtstd xdmdeigtstd
complex16: xeigtstz xdmdeigtstz

xdmdeigtsts: $(SDMDEIGTST) $(TMGLIB) $(LAPACKLIB) $(BLASLIB)
$(FC) $(FFLAGS) $(LDFLAGS) -o $@ $^

xdmdeigtstc: $(CDMDEIGTST) $(TMGLIB) $(LAPACKLIB) $(BLASLIB)
$(FC) $(FFLAGS) $(LDFLAGS) -o $@ $^

xdmdeigtstd: $(DDMDEIGTST) $(TMGLIB) $(LAPACKLIB) $(BLASLIB)
$(FC) $(FFLAGS) $(LDFLAGS) -o $@ $^

Expand All @@ -155,6 +160,7 @@ xeigtstz: $(ZEIGTST) $(DZIGTST) $(AEIGTST) $(TMGLIB) $(LAPACKLIB) $(BLASLIB)
$(FC) $(FFLAGS) $(LDFLAGS) -o $@ $^

$(SDMDEIGTST): $(FRC)
$(CDMDEIGTST): $(FRC)
$(DDMDEIGTST): $(FRC)
$(ZDMDEIGTST): $(FRC)
$(AEIGTST): $(FRC)
Expand Down Expand Up @@ -185,6 +191,8 @@ zchkee.o: zchkee.F
$(FC) $(FFLAGS_DRV) -c -o $@ $<
schkdmd.o: schkdmd.f90
$(FC) $(FFLAGS_DRV) -c -o $@ $<
cchkdmd.o: cchkdmd.f90
$(FC) $(FFLAGS_DRV) -c -o $@ $<
dchkdmd.o: dchkdmd.f90
$(FC) $(FFLAGS_DRV) -c -o $@ $<
zchkdmd.o: zchkdmd.f90
Expand Down
Loading

0 comments on commit c65a7a7

Please sign in to comment.