dmap->ComputeStructuredCoordinates( buf, ijk, pcoords );
double radius =
dmap->GetScalarComponentAsDouble( ijk[ 0 ], ijk[ 1 ], ijk[ 2 ], 0 );
- double rad_dis = _TDist( 1e-5, radius * 0.9 )( gen );
+ double rad_dis = _TDist( 1e-2, radius * 0.9 )( gen );
double the_dis = _TDist( 0, 3.14159265359 )( gen );
double phi_dis = _TDist( 0, 6.28318530718 )( gen );
buf[ 0 ] += rad_dis * std::sin( the_dis ) * std::cos( phi_dis );