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

tipsy output can be made in double precision (format = tipsy_double) even if MUSIC is run in single precision. This improves accuracy anyway since the single precision displacement vectors are converted to double before being converted to particle positions

Browse source
This commit is contained in:
Oliver Hahn 2015-06-14 19:43:24 +02:00
parent 381562c8b8
commit 86cd42f41f
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
plugins/output_tipsy.cc
View file

@ -843,7 +843,7 @@ public:
gh.cell_pos(ilevel, i, j, k, xx); gh.cell_pos(ilevel, i, j, k, xx);
//xx[coord] = fmod( (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) + 1.0, 1.0 ) - 0.5; //xx[coord] = fmod( (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) + 1.0, 1.0 ) - 0.5;
xx[coord] = (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) - 0.5; xx[coord] = (xx[coord]+(T_store)(*gh.get_grid(ilevel))(i,j,k)) - 0.5;
if( temp_data.size() < block_buf_size_ ) if( temp_data.size() < block_buf_size_ )
temp_data.push_back( xx[coord] ); temp_data.push_back( xx[coord] );
@ -1044,7 +1044,7 @@ public:
xx[coord] += 0.5*h; xx[coord] += 0.5*h;
//xx[coord] = fmod( (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) + 1.0, 1.0 ) - 0.5; //xx[coord] = fmod( (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) + 1.0, 1.0 ) - 0.5;
xx[coord] = (xx[coord]+(*gh.get_grid(ilevel))(i,j,k)) - 0.5; xx[coord] = (xx[coord]+(T_store)(*gh.get_grid(ilevel))(i,j,k)) - 0.5;
if( temp_data.size() < block_buf_size_ ) if( temp_data.size() < block_buf_size_ )
temp_data.push_back( xx[coord] ); temp_data.push_back( xx[coord] );
@ -1105,7 +1105,7 @@ int tipsy_output_plugin<double>::xdr_dump( XDR *xdrs, double*p )
namespace{ namespace{
output_plugin_creator_concrete< tipsy_output_plugin<float> > creator1("tipsy"); output_plugin_creator_concrete< tipsy_output_plugin<float> > creator1("tipsy");
#ifndef SINGLE_PRECISION //#ifndef SINGLE_PRECISION
output_plugin_creator_concrete< tipsy_output_plugin<double> > creator2("tipsy_double"); output_plugin_creator_concrete< tipsy_output_plugin<double> > creator2("tipsy_double");
#endif //#endif
} }