---

makeheavystar.C

Go to the documentation of this file.

// Kimberly Engle, November 1997

#include "node.h"

#ifdef TOOLBOX

// double_mass: Double the mass of a star to approximate a binary.
//                Only the mass is modified at this stage.

local void double_mass(node* original, real mass_ratio)
{
    // Set new masses.

    original->set_mass(mass_ratio*original->get_mass());
    //    putiq(original->get_log_story(), "mass_doubled", 1);
}

local void mkheavystar(node* b, real fraction_doubled, real lower_limit,
                       real mass_ratio)
{
 
    real sum = 0;
    b->set_mass(0);

    PRI(6); putrq(b->get_log_story(),
                  "fraction_of_masses_doubled", fraction_doubled);

    for_all_daughters(node, b, bi) {
        sum += fraction_doubled;
        if (sum >= 1) {
            sum -= 1;

            double_mass(bi, mass_ratio);
            putrq(bi->get_log_story(), "mass_doubled", 1);
        }
        else putiq(bi->get_log_story(), "mass_doubled", 0);

        b->inc_mass(bi->get_mass());
    }
}

void main(int argc, char ** argv)
{
    real fraction_doubled = 0.1;
    real lower_limit = 1;
    int random_seed = 0;
    real mass_ratio = 2;
    //    char seedlog[64];

    check_help();

    extern char *poptarg;
    int c;
    char* param_string = "f:l:q:s:";

    while ((c = pgetopt(argc, argv, param_string)) != -1)
        switch(c) {

            case 'f': fraction_doubled = atof(poptarg);
                      break;
            case 'l': lower_limit = atof(poptarg);
                      break;
            case 'q': mass_ratio = atof(poptarg);
                      break;
            case 's': random_seed = atoi(poptarg);
                      break;
            case '?': params_to_usage(cerr, argv[0], param_string);
                      get_help();
                      exit(1);
        }

    if (fraction_doubled < 0 || fraction_doubled > 1)
        err_exit("mkheavystar: Illegal doubling fraction");
    //    if (lower_limit <= 0 || lower_limit > 1)
    //  err_exit("mkheavystar: Illegal mass ratio limit");

    node* b;
    b = get_node(cin);

    b->log_history(argc, argv);

    //    int actual_seed = srandinter(random_seed);

    //    sprintf(seedlog,
    //            "       random number generator seed = %d",actual_seed);
    //    b->log_comment(seedlog);

    mkheavystar(b, fraction_doubled, lower_limit, mass_ratio);

    put_node(cout, *b);

}
#endif

---