Hello All
The script below exceeds molpdf limits immediately when run this automodel job on parallel,
but run starts with a low molpdf and ends well when commenting the line "a.use_parallel_job(j)"
Any hint of how to solve will be appreciate?
Mario
from modelito import MyModel
from modeller.parallel import *
icycle=50
log.minimal() # request minimal output
env = Environ() # create a new MODELLER environment to build this model in
icycle=str(int(icycle)+1)
# directories for input atom files
env.io.atom_files_directory = ['.', 'trajs/']
# Give less weight to all soft-sphere restraints:
env.schedule_scale = physical.values(default=1.0, soft_sphere=0.7)
# Very thorough VTFM optimization:
j = Job()
for i in range(1):
j.append(local_slave())
# Comparative modeling by the AutoModel class
env.io.hetatm=True
a = MyModel(env,
alnfile = 'test.ali', # alignment filename
knowns = ('ro','odo'), # codes of the templates
sequence = 'TEST') # code of the target
a.starting_model= int(icycle) # index of the first model
a.ending_model = int(icycle)+ 0 # index of the last model
# (determines how many models to calculate)
a.library_schedule = autosched.slow
#a.restraints.write(file='restraints.rsr')
# Thorough MD optimization:
a.md_level = refine.slow
# Repeat the whole cycle 2 times and do not stop unless obj.func. > 1E6
a.use_parallel_job(j)
a.repeat_optimization = 3
a.max_molpdf = 1e6
a.make() # do the actual comparative modeling