InitialData::InitialData(const string& radius_file,
const string& density_file,
const string& temperature_file,
const string& velocity_file):
radius_list(load_txt(radius_file)),
radius_mid(mid_array(radius_list)),
density_list(decapitate(load_txt(density_file))),
temperature_list(decapitate(load_txt(temperature_file))),
velocity_list(load_txt(velocity_file)),
tracers_list(get_composition_data()) {}
map<string,vector<double> > get_composition_data(void)
{
map<string,vector<double> > res;
const map<string,pair<double,double> > atomic_properties = generate_atomic_properties();
for(map<string,pair<double,double> >::const_iterator it =
atomic_properties.begin();
it != atomic_properties.end(); ++it)
res[it->first] = decapitate(load_txt(string("tracer_")+it->first+".txt"));
return res;
}
int scene::load(char *p_dir, char *p_file, SCENE_TYPE type)
{
int ret = 0;
switch(type) {
case SCENE_B2M:
ret = load_b2m(p_dir,p_file);
break;
case SCENE_3DS:
ret = load_3ds(p_dir,p_file);
break;
case SCENE_TXT:
ret = load_txt(p_dir,p_file);
break;
}
load_insert();
return(ret);
}
