@@ -0,0 +1,1439 @@
+# coding: utf-8%0A%0A%22%22%22 Check the status of frequency mapping. %22%22%22%0A%0Afrom __future__ import division, print_function%0A%0A__author__ = %22adrn %3Cadrn@astro.columbia.edu%3E%22%0A%0Afrom streammorphology.ensemble import read_allkld%0A%0Adef main(path, nkld):%0A    # read allfreqs into structured array%0A    d = read_allkld(path, nkld=nkld)%0A%0A    ndone = (d%5B'status'%5D %3E 0).sum()%0A    nsuccess = (d%5B'status'%5D == 1).sum()%0A    nfail_energy = (d%5B'status'%5D == 2).sum()%0A    nfail_integrate = (d%5B'status'%5D == 3).sum()%0A    nfail_kld = (d%5B'status'%5D == 4).sum()%0A%0A    print(%22Number of orbits: %7B%7D%22.format(len(d)))%0A    print(%22Done processing: %7B%7D%22.format(ndone))%0A    print(%22Successful: %7B%7D%22.format(nsuccess))%0A    print(%22Failures: %7B%7D%22.format(ndone-nsuccess))%0A    print(%22%5Ct Energy failure: %7B%7D%22.format(nfail_energy))%0A    print(%22%5Ct Integration failure: %7B%7D%22.format(nfail_integrate))%0A    print(%22%5Ct KLD failure: %7B%7D%22.format(nfail_kld))%0A%0Aif __name__ == '__main__':%0A    from argparse import ArgumentParser%0A%0A    # Define parser object%0A    parser = ArgumentParser(description=%22%22)%0A    parser.add_argument(%22-p%22, %22--path%22, dest=%22path%22, required=True,%0A                        help=%22Path to a Numpy memmap file containing the results %22%0A                             %22of frequency mapping.%22)%0A    parser.add_argument(%22--nkld%22, dest=%22nkld%22, required=True,%0A                        help=%22Number of times the KLD was computed and saved.%22)%0A%0A    args = parser.parse_args()%0A%0A    main(args.path, nkld=args.nkld)%0A

add status script for ensemble