meme_em['niters'] = start.niters
meme_em['em_time'] = start.em_time
meme_em['cons'] = start.cons_after_em
meme_em['nsites'] = start.nsites
meme_em['sig'] = start.sig
logging.info(
'%s: cons0=%20s; nsites0=%3d; niters=%4d; elapsed=%7.1fs; per iteration=%6.2fs; cons=%20s; nsites0=%3d; sig=%e',
data_set, start.cons0, start.nsites0, start.niters, start.em_time, start.em_time /
start.niters, start.cons, start.nsites, start.sig
)
meme_em.append()
h5file.close() # Close (and flush) the HDF5 file
if '__main__' == __name__:
#
# parse options
#
parser = OptionParser()
test_data.add_options(parser)
stempy.add_options(parser)
options, args = parser.parse_args()
usage = 'USAGE: %s <options> <h5 file>' % sys.argv[0]
if len(args) != 1:
raise RuntimeError(usage)
if len(args) < 1:
raise RuntimeError(usage)
filename = args[0]
run(options)