---

sys_stats.C

Go to the documentation of this file.

       //=======================================================//    _\|/_
      //  __  _____           ___                    ___       //      /|\ ~
     //  /      |      ^     |   \  |         ^     |   \     //          _\|/_
    //   \__    |     / \    |___/  |        / \    |___/    //            /|\ ~
   //       \   |    /___\   |  \   |       /___\   |   \   // _\|/_
  //     ___/   |   /     \  |   \  |____  /     \  |___/  //   /|\ ~
 //                                                       //            _\|/_
//=======================================================//              /|\ ~

// The help information below is copied directly from the dyn version.
// The only difference between this program and the dyn version is that
// possible GRAPE and dstar extensions are used here.


#include "hdyn.h"

#ifdef TOOLBOX

#include "../evolve/kira_grape_include.C"       // GRAPE-specific functions:
                                                // Invoke by setting USE_GRAPE
                                                // at compilation time.

// Main code follows dyn/util/sys_stats.C version.

main(int argc, char **argv)
{
    check_help();

    bool binaries, long_binary_output, B_flag, verbose, out, n_sq, calc_e;
    int which_lagr;

    // (Defaults are set in dyn/util/sys_stats.C)

    if (!parse_sys_stats_main(argc, argv,
                              which_lagr,
                              binaries, long_binary_output, B_flag,
                              calc_e, n_sq, out, verbose)) {
        get_help();
        exit(1);
    }

    if (!n_sq) {

#ifndef USE_GRAPE
        calc_e = false;
#endif

    }

    // Loop over input until no more data remain.

    hdyn *b;
    int i = 0;

    while (b = get_hdyn(cin)) {

        check_addstar(b);
        check_set_external(b, true);    // true ==> verbose output
        cerr << endl;

        if (B_flag || check_kira_flag(b, "kira_evolve_binaries"))
            b->set_use_dstar(true);

        if (i++ > 0) cerr << endl;

        // Simply invoke the dyn function with appropriate hdyn
        // extensions, as in kira_log.C.

        sys_stats(b,
                  0.5,                  // energy cutoff
                  verbose,
                  binaries,
                  long_binary_output,
                  which_lagr,
                  true,                 // print_time
                  calc_e,
                  n_sq,
                  kira_calculate_energies,
                  print_dstar_params,
                  print_dstar_stats);

        if (out) {
            b->log_history(argc, argv);
            put_node(cout, *b);
        }

        rmtree(b);      // causes core dump if B_flag is enabled...
    }
}

#endif

---