Commit 94915e6e by Francois Gygi

Fix gaussian normalization in SpectrumCmd.C

parent 9f4d98fd
......@@ -210,7 +210,7 @@ int SpectrumCmd::action(int argc, char **argv)
for ( int ie = 0; ie < np; ie++ )
{
const double t = ( emin + ie * de - delta_e ) / width;
sp[ispin][ie] += w * width * sqrt(M_PI) * exp(-t*t);
sp[ispin][ie] += w * ( sqrt(M_PI) / width ) * exp(-t*t);
}
// only send if not on pe 0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment