From 2f2d0120af70ae544e6cd8b7cf51954d52e197ea Mon Sep 17 00:00:00 2001 From: "Martin D. Weinberg" Date: Wed, 6 Nov 2024 13:04:33 -0500 Subject: [PATCH] Missing removal of 2Pi factors from PolarBasis --- src/PolarBasis.cc | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/PolarBasis.cc b/src/PolarBasis.cc index 67f2e454f..876432fbf 100644 --- a/src/PolarBasis.cc +++ b/src/PolarBasis.cc @@ -401,8 +401,8 @@ void * PolarBasis::determine_coefficients_thread(void * arg) { // For biorthogonal density component and normalization // - constexpr double norm0 = 2.0*M_PI * 0.5*M_2_SQRTPI/M_SQRT2; - constexpr double norm1 = 2.0*M_PI * 0.5*M_2_SQRTPI; + constexpr double norm0 = 0.5*M_2_SQRTPI/M_SQRT2; + constexpr double norm1 = 0.5*M_2_SQRTPI; double r, r2, facL=1.0, fac1, fac2, phi, mass; double xx, yy, zz; @@ -1044,8 +1044,8 @@ void PolarBasis::multistep_update(int from, int to, Component *c, int i, int id) // For biorthogonal density component and normalization // - constexpr double norm0 = 2.0*M_PI * 0.5*M_2_SQRTPI/M_SQRT2; - constexpr double norm1 = 2.0*M_PI * 0.5*M_2_SQRTPI; + constexpr double norm0 = 0.5*M_2_SQRTPI/M_SQRT2; + constexpr double norm1 = 0.5*M_2_SQRTPI; double mass = c->Mass(i) * component->Adiabatic();