def test_repnums(self):
if TuneRex is None:
self.skipTest("test requires the rpy2 module")
TuneRex.debug = True
rn = self.read_replicanums_file()
changed, params = TuneRex.tune_params_flux(rn, self.temps)
print params
def test_repnums(self):
if TuneRex is None:
self.skipTest("test requires the rpy2 module")
TuneRex.debug=True
rn=self.read_replicanums_file()
changed,params=TuneRex.tune_params_flux(rn,self.temps)
print params
def test_flux(self):
if TuneRex is None:
self.skipTest("test requires the rpy2 module")
# self.replica.tune_data['dumb_scale']=0.5
# for i in range(1000):
# self.replica.replica_exchange()
# self.replica.write_rex_stats()
del self.replica.tune_data['dumb_scale']
# self.replica.tune_data['CvMethod']='constant'
# self.replica.tune_data['goodMethod']='step'
# self.replica.tune_data['badMethod']='step'
self.replica.tune_data['rate'] = 500
self.replica.tune_data['method'] = 'flux'
#self.replica.tune_temps = False
for i in range(5000):
self.replica.replica_exchange()
self.replica.write_rex_stats()
# if os.path.exists('pup'):
# os.rename('pup','pup_%d' % i)
self.replica.tune_data['dumb_scale'] = 0.1
replicanums = self.read_replicanums_file()
indicators = TuneRex.compute_indicators(replicanums, subs=1, start=0)
print array([sum(ind) / float(len(ind)) for ind in indicators])
def test_flux(self):
if TuneRex is None:
self.skipTest("test requires the rpy2 module")
#self.replica.tune_data['dumb_scale']=0.5
#for i in range(1000):
# self.replica.replica_exchange()
# self.replica.write_rex_stats()
del self.replica.tune_data['dumb_scale']
#self.replica.tune_data['CvMethod']='constant'
#self.replica.tune_data['goodMethod']='step'
#self.replica.tune_data['badMethod']='step'
self.replica.tune_data['rate']=500
self.replica.tune_data['method']='flux'
#self.replica.tune_temps = False
for i in range(5000):
self.replica.replica_exchange()
self.replica.write_rex_stats()
#if os.path.exists('pup'):
# os.rename('pup','pup_%d' % i)
self.replica.tune_data['dumb_scale']=0.1
replicanums=self.read_replicanums_file()
indicators = TuneRex.compute_indicators(replicanums, subs=1,start=0)
print array([sum(ind)/float(len(ind)) for ind in indicators])
