debug

Author: js1019

EXTERNAL/makefile

@@ -4,9 +4,12 @@ ifeq ($(INTEL), 1)
   OBJ = DIRECTSOL/pevsl_pardiso.o ITERSOL/chebiter.o
   INCLUDES = -I../INC -I. -I$(MKLROOT)/include
 else
-  include DIRECTSOL/makefile_mumps.in
-  OBJ = DIRECTSOL/pevsl_mumps.o ITERSOL/chebiter.o
-  INCLUDES = -I../INC -I. $(MUMPS_INC)
+  #include DIRECTSOL/makefile_mumps.in
+  #OBJ = DIRECTSOL/pevsl_mumps.o ITERSOL/chebiter.o
+  #INCLUDES = -I../INC -I. $(MUMPS_INC)
+  #include DIRECTSOL/makefile_mumps.in
+  OBJ = ITERSOL/chebiter.o
+  INCLUDES = -I../INC -I. 
 endif
 
 %.o : %.c

FORTRAN/pevsl_f90.c

@@ -282,7 +282,7 @@ void PEVSL_FORT(pevsl_zlanbounds)(uintptr_t *pevslf90, int *mlan, int *nsteps, d
   pEVSL_ParvecCreate(N, n, nfirst, pevsl->comm, &viinit);
   pEVSL_ParvecRand(&viinit);
   /*------------------- Lanczos Bounds */
-  pEVSL_ZLanTrbounds(pevsl, *mlan, *nsteps, 1e-8, &vrinit, &viinit, 1, lmin, lmax, NULL);
+  pEVSL_ZLanTrbounds(pevsl, *mlan, *nsteps, 1e-8, &vrinit, &viinit, 1, lmin, lmax, stdout);
 
   pEVSL_ParvecFree(&vrinit);
   pEVSL_ParvecFree(&viinit);

README.md

@@ -1,9 +1,5 @@
 -----------------------------------------------------------------------
-                  pEVSL (parallel EigenValue SLicing)
-
------------------------------------------------------------------------
-Modification for complex Hermitian systems
-
+                  pEVSL (parallel EigenValue SLcing)
 -----------------------------------------------------------------------
   Latest changes made on: 
 -----------------------------------------------------------------------

SRC/lantrbnd.c

@@ -461,6 +461,7 @@ int pEVSL_ZLanTrbounds(pevsl_Data *pevsl, int lanm, int maxit, double tol,
   /*-----------------------------------------------------------------------*
    * *thick restarted* Lanczos step
    *-----------------------------------------------------------------------*/
+  
   if (do_print) {
     fprintf(fstats, " LanTR for bounds: dim %d, maxits %d\n", lanm, maxit);
   }
@@ -545,10 +546,13 @@ int pEVSL_ZLanTrbounds(pevsl_Data *pevsl, int lanm, int maxit, double tol,
     /* z references the 1st columns of Z */
     pEVSL_ParvecsGetParvecShell(Zr, 0, zr);
     pEVSL_ParvecsGetParvecShell(Zi, 0, zi);
+
+    pEVSL_fprintf0(rank,fstats,"check 0 \n");
     /* B norm */
     pEVSL_ZMatvecB(pevsl, vr, vi, zr, zi);
     pEVSL_ParvecZDot(vr, vi, zr, zi, &tr, &ti);
     t = 1.0 / sqrt(tr);
+    pEVSL_fprintf0(rank,fstats,"check 1 \n");
     /* check! ti != 0.0 ? */
 
     /* z = B*v */

makefile

@@ -2,8 +2,8 @@ default:
 	( cd SRC           ; make )
 	( cd EXTERNAL      ; make )
 	( cd TESTS/Lap     ; mkdir -p OUT ; make)
-	( cd TESTS/Gen_Lap ; mkdir -p OUT ; make)
-#	( cd TESTS/MM      ; mkdir -p OUT ; make)
+#	( cd TESTS/Gen_Lap ; mkdir -p OUT ; make)
+	( cd TESTS/MM      ; mkdir -p OUT ; make)
 all:
 	( cd SRC           ; make )
 	( cd EXTERNAL      ; make )

makefile.in

@@ -1,4 +1,4 @@
-INTEL = 1
+INTEL = 0
 
 # Compilers
 ifeq ($(INTEL), 1)
@@ -12,7 +12,8 @@ LINK = $(CC)
 # Common flags
 # Debug mode
 PEVSL_DEBUG = -DPEVSL_DEBUG
-FLAGS = -DUNIX -O3 -g -Wall $(PEVSL_DEBUG) $(TIMING)
+#FLAGS = -DUNIX -O3 -g -Wall $(PEVSL_DEBUG) $(TIMING)
+FLAGS = -DUNIX -g -Wall $(PEVSL_DEBUG) $(TIMING)
 
 ##OBJ_TIME = mactime.o
 OBJ_TIME = timing.o 
@@ -25,6 +26,7 @@ ifeq ($(INTEL), 1)
 FLAGS += -DUSE_MKL
 LIBLAPACK =  -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a ${MKLROOT}/lib/intel64/libmkl_intel_thread.a ${MKLROOT}/lib/intel64/libmkl_core.a ${MKLROOT}/lib/intel64/libmkl_blacs_intelmpi_lp64.a -Wl,--end-group -liomp5 -lpthread -lm -ldl
 else
-LIBLAPACK = /home/li50/workspace/OpenBLAS/libopenblas.a -lpthread
-LIB0 = -lgfortran -lm
+#LIBLAPACK = /home/li50/workspace/OpenBLAS/libopenblas.a -lpthread
+LIBLAPACK =  -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a ${MKLROOT}/lib/intel64/libmkl_gnu_thread.a ${MKLROOT}/lib/intel64/libmkl_core.a ${MKLROOT}/lib/intel64/libmkl_blacs_intelmpi_lp64.a -Wl,--end-group -lgomp -lpthread -lm -ldl
+LIB0 = -lgfortran #-lm
 endif