diff --git a/src/SpectrumCmd.C b/src/SpectrumCmd.C index 9f36c82..dad2dab 100644 --- a/src/SpectrumCmd.C +++ b/src/SpectrumCmd.C @@ -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