Esempio n. 1
0
        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)