#!/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')
