mirror/MUSIC
1
0
Fork 0
mirror of https://github.com/cosmo-sims/MUSIC.git synced 2024-09-19 17:03:46 +02:00
Code Releases Activity

added Gadget header flag for double vs single precision

Browse source
This commit is contained in:
Oliver Hahn 2014-09-23 10:02:32 +02:00
parent 99581697f0
commit cac22512ed
1 changed files with 12 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
plugins/output_gadget2.cc
View file

@ -16,7 +16,7 @@
#include "mg_interp.hh"
#include "mesh.hh"
const int empty_fill_bytes = 60;
const int empty_fill_bytes = 56;
template< typename T_store=float >
class gadget2_output_plugin : public output_plugin
@ -60,7 +60,8 @@ protected:
int flag_stellarage;
int flag_metals;
unsigned int npartTotalHighWord[6];
int flag_entropy_instead_u;
int flag_entropy_instead_u;
int flag_doubleprecision;
char fill[empty_fill_bytes];
}header;
@ -872,6 +873,15 @@ public:
header_.npartTotalHighWord[i] = 0;
header_.mass[i] = 0.0;
}
if( typeid(T_store)==typeid(float) )
header_.flag_doubleprecision = 0;
else if( typeid(T_store)==typeid(double) )
header_.flag_doubleprecision = 1;
else{
LOGERR("Internal error: gadget-2 output plug-in called for neither \'float\' nor \'double\'");
throw std::runtime_error("Internal error: gadget-2 output plug-in called for neither \'float\' nor \'double\'");
}
YHe_ = cf.getValueSafe<double>("cosmology","YHe",0.248);
gamma_ = cf.getValueSafe<double>("cosmology","gamma",5.0/3.0);