From fb1691be80adfff153df72cd82471de68700c582 Mon Sep 17 00:00:00 2001 From: Jonathon Misiewicz Date: Wed, 9 Oct 2024 12:15:09 -0400 Subject: [PATCH] Hotfixes --- src/madness/mra/gfit.h | 2 +- src/madness/mra/operator.h | 4 +++- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/src/madness/mra/gfit.h b/src/madness/mra/gfit.h index c457a4af4f7..24782c1d002 100644 --- a/src/madness/mra/gfit.h +++ b/src/madness/mra/gfit.h @@ -522,7 +522,7 @@ class GFit { // in the moment list // // cj = - const long nmom = 1; + const long nmom = 0; if (nmom > 0) { Tensor q(4), qg(4); double range = sqrt(-log(1e-6)/expnt[nmom-1]); diff --git a/src/madness/mra/operator.h b/src/madness/mra/operator.h index f9cafbca971..62a2c84c15e 100644 --- a/src/madness/mra/operator.h +++ b/src/madness/mra/operator.h @@ -245,8 +245,10 @@ namespace madness { Tensor coeff=fit.coeffs(); Tensor expnt=fit.exponents(); + // WARNING! More fine-grained control over the last argument is needed. + // This is a hotfix. if (bc(0,0) == BC_PERIODIC) { - fit.truncate_periodic_expansion(coeff, expnt, cell_width.max(), false); + fit.truncate_periodic_expansion(coeff, expnt, cell_width.max(), true); } return std::make_pair(coeff,expnt);