GALPROPhttps://galprop.stanford.edu/forum/ CR from DM annihilation with a cuspy DM profile.https://galprop.stanford.edu/forum/viewtopic.php?f=7&t=79 Page 1 of 1

 Author: jhhuh [ Fri Sep 25, 2009 4:38 am ] Post subject: CR from DM annihilation with a cuspy DM profile. I used galprop v.50p for my recent paper about CR from decaying DM.While I used it, I found some bug in the code averaging DM density.( "DM_profile_av()" in "gen_DM_source.cc")The calculation with 3-d coordinate or non-cuspy profile is OK. Butthe calculation with cuspy profile like NFW in 2-d coordinate makeunacceptably large CR fluxes. It maybe the problem in averaging DMnear the origin (Jacobi factor become significant in this region.).By simply changing "int nuse" into "double ruse" and "nuse++" into "ruse+=rr",code produce appropriate result.~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ double Galprop::DM_profile_av(double r,double z,double dr,double dz,double dzz) { double DM_profile_av_=0.0; int nuse=0; double ruse=0; for (double zz=z-dz/2.; zz<=z+dz/2.; zz+=dzz) for (double rr=r-dr/2.; rr<=r+dr/2.; rr+=dr/10.) { if (rr<0.) continue; DM_profile_av_+=rr*DM_profile(rr,0,zz); ruse+=rr; } return (DM_profile_av_/ruse); }~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

 Page 1 of 1 All times are UTC - 8 hours [ DST ] Powered by phpBB® Forum Software © phpBB Grouphttps://www.phpbb.com/