Wgap set incorrectly in subset case

SebastianAchilles · Jun 24, 2013 · c00b8ab · c00b8ab
1 parent 612809a
commit c00b8ab
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/SRC/mrrr_val.c b/SRC/mrrr_val.c
@@ -563,6 +563,7 @@ int mrrr_val(char *jobz, char *range, int nthreads, in_t *Dstruct,
   	W[jee] = W[je];
   	Wcopy[jee] = Wcopy[je];
   	Werr[jee] = Werr[je];
+	Wgap[jee] = Wgap[je];
   	Windex[jee] = Windex[je];
   	iblock[jee] = iblock[je];
   	jee++;

diff --git a/SRC/process_singleton.c b/SRC/process_singleton.c
@@ -137,7 +137,7 @@ int PMR_process_s_task(singleton_t *sng, int tid, counter_t *num_left,
     right   = Wshifted[i] + Werr[i];
     i_local = Windex[i];
     r       = 0;
-    
+
     /* compute left and right gap */
     if (i == bl_begin)
       lgap = DBL_EPSILON * fmax(fabs(left), fabs(right));

diff --git a/make.inc b/make.inc
@@ -4,8 +4,8 @@ CC = gcc
 FC = gfortran
 
 # Compiler flags
-CFLAGS = -O3 -pthread 
-FFLAGS = -O3 -funderscoring
+CFLAGS = -O3 -g -pthread 
+FFLAGS = -O3 -g -funderscoring
 
 # Archiver and flags used when building the archive
 AR = /usr/bin/ar