merged with phmc

EVS.data

@@ -0,0 +1,4 @@
+ EV_min  Ev_max  n  s_low  s_max  normalisation 
+ 0.218281 3.401855 23 0.064000 1.000000 0.494939 
+ 0.138196 2.587464 23 0.053000 1.000000 0.567508 
+ 0.213994 2.519209 19 0.084000 1.000000 0.575145 

INOUT.data

@@ -0,0 +1 @@
+ 0.031000 1.000000 0.494939 

Makefile.in

@@ -29,30 +29,33 @@ COMPILE = ${CC} $(DEFS) $(INCLUDES) -o $@ ${CFLAGS}
 
 SMODULES = Hopping_Matrix_nocom Hopping_Matrix tm_operators
 
+PMODULES = init_chi_spinor_field init_chi_copy reweighting_factor_nd \
+        chebyshev_polynomial_nd Ptilde_nd  Nondegenerate_Matrix \
+	init_bispinor_field eigenvalues_bi max_eigenvalues_bi \
+	update_tm_nd hybrid_nondegenerate_update
+
 MODULES = read_input gamma hybrid_update observables start \
 	expo 2mn_integrator get_staples update_backward_gauge \
 	measure_rectangles get_rectangle_staples  \
 	test/check_geometry test/check_xchange invert_eo \
 	ext_integrator polyakov_loop getopt sighandler \
 	source_generation boundary io update_tm ranlxd  \
 	mpi_init linsolve deriv_Sb \
-	reweight_kappac xchange_deri geometry_eo linalg_eo \
-	init_bispinor_field init_moment_field init_gauge_tmp \
+	xchange_deri geometry_eo linalg_eo \
+	init_moment_field init_gauge_tmp \
 	derivative_psf xchange_field xchange_gauge \
 	init_gauge_field init_geometry_indices init_spinor_field \
-	init_dirac_halfspinor xchange_halffield stout_smear
+	init_dirac_halfspinor xchange_halffield stout_smear 
 
-### eigenvalues_bi max_eigenvalues_bi Nondegenerate_Matrix \
+PROGRAMS = hmc_tm benchmark invert gwc2ildg \
+	ildg2gwc single2double double2single reducenoise gen_sources 
 
-PROGRAMS = hmc_tm thermal_cycle_tm benchmark invert gwc2ildg ildg2gwc single2double double2single reducenoise gen_sources
-### test/test_eigenvalues
-ALLOBJ = ${MODULES} ${PROGRAMS} ${SMODULES}
+ALLOBJ = ${MODULES} ${PROGRAMS} ${SMODULES} ${PMODULES}
 SUBDIRS = linalg solver
 
 .SUFFIXES:
 
-all: Makefile all-recursive dep ${LINKLIBS} hmc_tm invert 
-### test/test_eigenvalues 
+all: Makefile all-recursive dep ${LINKLIBS} hmc_tm invert
 
 .NOTPARALLEL: 
 
@@ -69,22 +72,29 @@ ${addsuffix .o, ${MODULES}}: %.o: ${srcdir}/%.c %.d Makefile
 ${addsuffix .o, ${SMODULES}}: %.o: ${srcdir}/%.c %.d Makefile
 	${COMPILE} ${SOPTARGS} -c $< 
 
+${addsuffix .o, ${PMODULES}}: %.o: ${srcdir}/%.c %.d Makefile
+	${COMPILE} ${OPTARGS} -c $< 
+
 ${addsuffix .o, ${PROGRAMS}}: %.o: ${srcdir}/%.c %.d Makefile
 	${COMPILE} ${OPTARGS} -c $<
 
 ${PROGRAMS}: %: %.o ${addsuffix .o, ${MODULES}} ${addsuffix .o, ${SMODULES}} ${LINKLIBS}
 	 ${LINK} ${addsuffix .o, ${MODULES}} ${addsuffix .o, ${SMODULES}} $@.o $(LIBS)
 
+phmc_tm: %: %.o ${addsuffix .o, ${MODULES}} ${addsuffix .o, ${PMODULES}} ${addsuffix .o, ${SMODULES}} ${LINKLIBS}
+	 ${LINK} ${addsuffix .o, ${MODULES}} ${addsuffix .o, ${PMODULES}} ${addsuffix .o, ${SMODULES}} $@.o $(LIBS)
+
+
 dep: $(addsuffix .d,$(ALLOBJ))
 
 compile-clean: compile-clean-recursive Makefile
 	rm -f *.o *.d
 
 clean: clean-recursive Makefile
-	rm -f hmc_tm *.o *.d
-
+	rm -f hmc_tm invert phmc_tm*.o *.d
 distclean: distclean-recursive Makefile
 	rm -f hmc_tm hybrid *.o *.d *~ Makefile config.log config.status fixed_volume.h
+	rm -f phmc_tm
 	rm -f config.h
 
 .PHONY: all clean compile-clean distclean dep \