---

makecube.C

Go to the documentation of this file.

#include "dyn.h"

#ifdef TOOLBOX

local void  mkcube(dyn * root, int n, int u_flag)
{
    dyn  * bi;

    root->set_mass(1);
    real pmass = 1.0 / n;

    for (bi = root->get_oldest_daughter(); bi != NULL;
         bi = bi->get_younger_sister()) {
        bi->set_mass(pmass);    
        bi->set_pos(vector(randinter(-1,1),randinter(-1,1),randinter(-1,1)));
        bi->set_vel(vector(randinter(-1,1),randinter(-1,1),randinter(-1,1)));
    }

//  Transform to center-of-mass coordinates, and optionally
//  scale to standard parameters.

    root->to_com();
    putrq(root->get_log_story(), "initial_mass", 1.0);

    if (! u_flag && n > 1) {

        real kinetic, potential;

        get_top_level_energies(root, 0.0, potential, kinetic);
        scale_virial(root, -0.5, potential, kinetic);   // scales kinetic
        real energy = kinetic + potential;
        scale_energy(root, -0.25, energy);              // scales energy
        putrq(root->get_dyn_story(), "initial_total_energy", -0.25);
        putrq(root->get_log_story(), "initial_rvirial", 1.0);
    }
}

#define  SEED_STRING_LENGTH  60

main(int argc, char ** argv) {
    int  i;
    int  n;
    int  input_seed, actual_seed;
    int  c_flag = FALSE;
    int  n_flag = FALSE;
    int  s_flag = FALSE;
    int  i_flag = FALSE;
    int  l_flag = FALSE;
    int  o_flag = FALSE;
    int  u_flag = FALSE;   // if TRUE, leave scaling untouched

    char  *comment;
    char  seedlog[SEED_STRING_LENGTH];

    real  cube_size = 2;

    check_help();

    extern char *poptarg;
    int c;
    char* param_string = "c:iln:os:u";

    while ((c = pgetopt(argc, argv, param_string)) != -1)
        switch(c)
            {
            case 'c': c_flag = TRUE;
                      comment = poptarg;
                      break;
            case 'i': i_flag = TRUE;
                      break;
            case 'l': l_flag = TRUE;
                      break;
            case 'n': n_flag = TRUE;
                      n = atoi(poptarg);
                      break;
            case 'o': o_flag = TRUE;
                      break;
            case 's': s_flag = TRUE;
                      input_seed = atoi(poptarg);
                      break;
            case 'u': u_flag = TRUE;
                      break;
            case '?': params_to_usage(cerr, argv[0], param_string);
                      get_help();
                      exit(1);
            }            
    
    if (n_flag == FALSE) {
        cerr << "makecube: must specify the number # of";
        cerr << " particles with -n#\n";
        exit(1);
    }
    
    if (n < 1) {
        cerr << "mkcube: n < 1 not allowed\n";
        exit(1);
    }

    dyn *b, *by, *bo;
    b = new dyn();
    bo = new dyn();
    if (i_flag) 
        bo->set_label(1);
    b->set_oldest_daughter(bo);
    bo->set_parent(b);
    if (l_flag)
        putrq(b->get_dyn_story(), "cube_size", cube_size);

    for (i = 1; i < n; i++) {
        by = new dyn();
        if (i_flag) by->set_label(i+1);
        by->set_parent(b);
        bo->set_younger_sister(by);
        by->set_elder_sister(bo);
        bo = by;
    }

    if (c_flag) b->log_comment(comment);

    b->log_history(argc, argv);

    if (s_flag == FALSE) input_seed = 0;
    actual_seed = srandinter(input_seed);

    if (o_flag) cerr << "mkcube: random seed = " << actual_seed << endl;

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

    if (n > 0) mkcube(b, n, u_flag);

    put_node(cout, *b);
}

#endif

/* end of: mkcube.c */

---