#!/usr/bin/env python import pypar import sys, math, Numeric, MLab, os, pickle, numarray from MontePython import MontePython try: ifile = sys.argv[1] except: ifile = '' # let MontePython instance make its own file d = MontePython(pypar, ifile) # 1 particle and alpha=0.5 yields a harmonic oscillator with # energy == 3/2: d.params[0] = 0.5 d.reset_params() #d.silent = True # to avoid too much noise if not os.access('checkpoint.dat', os.F_OK): d.blocks = 5 d.tau_list = [ [0.001, 100, d.blocks], ] if os.access('results_' + d.outfile, os.F_OK): os.unlink('results_' + d.outfile) else: d.blocks = d.tau_list[0][2] blocks = d.blocks loc_tmp_nrg = Numeric.empty(3, 'f')
#!/usr/bin/env python import pypar, sys, math, Numeric, MLab, os, pickle, numarray, copy from MontePython import MontePython try: ifile = sys.argv[1] except: ifile = '' # let MontePython instance make its own file d = MontePython(pypar,ifile) d.reset_params() #d.silent = True # to avoid too muvh noise d.blocks = 10 tau_list = [#[0.0015/20.,100,d.blocks], [0.0010/10.,25,d.blocks], #[0.0007/20.,200,d.blocks], #[0.0005/20.,200,d.blocks], ] if not os.access('checkpoint.dat',os.F_OK): d.tau_list = tau_list #if os.access('results_'+d.outfile,os.F_OK): # os.unlink('results_'+d.outfile) else: d.blocks = d.tau_list[0][2] blocks = d.blocks loc_tmp_nrg = Numeric.empty(3,'f')
#!/usr/bin/env python import pypar import sys, math, Numeric, MLab, os, pickle, numarray from MontePython import MontePython try: ifile = sys.argv[1] except: ifile = '' # let MontePython instance make its own file d = MontePython(pypar,ifile) # 1 particle and alpha=0.5 yields a harmonic oscillator with # energy == 3/2: d.params[0] = 0.5 d.reset_params() #d.silent = True # to avoid too much noise if not os.access('checkpoint.dat',os.F_OK): d.blocks = 5 d.tau_list = [[0.001,100,d.blocks],] if os.access('results_'+d.outfile,os.F_OK): os.unlink('results_'+d.outfile) else: d.blocks = d.tau_list[0][2] blocks = d.blocks loc_tmp_nrg = Numeric.empty(3,'f') tmp_nrg = Numeric.empty(3*d.numproc,'f') tmp_nrg.shape = (d.numproc,3)
#!/usr/bin/env python import pypar, sys, math, Numeric, MLab, os, pickle, numarray, copy from MontePython import MontePython try: ifile = sys.argv[1] except: ifile = '' # let MontePython instance make its own file d = MontePython(pypar, ifile) d.reset_params() #d.silent = True # to avoid too muvh noise d.blocks = 10 tau_list = [ #[0.0015/20.,100,d.blocks], [0.0010 / 10., 25, d.blocks], #[0.0007/20.,200,d.blocks], #[0.0005/20.,200,d.blocks], ] if not os.access('checkpoint.dat', os.F_OK): d.tau_list = tau_list #if os.access('results_'+d.outfile,os.F_OK): # os.unlink('results_'+d.outfile) else: d.blocks = d.tau_list[0][2] blocks = d.blocks loc_tmp_nrg = Numeric.empty(3, 'f')