From 7af8c77a5260dc294237c91140c8fc85ae8c8825 Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Thu, 24 Apr 2025 11:37:14 -0400
Subject: [PATCH] Normalise

---
 examples/Example_Laplacian_smearing.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/examples/Example_Laplacian_smearing.cc b/examples/Example_Laplacian_smearing.cc
index 9780e8a0..5c037461 100644
--- a/examples/Example_Laplacian_smearing.cc
+++ b/examples/Example_Laplacian_smearing.cc
@@ -93,10 +93,13 @@ int main(int argc, char ** argv)
   Real coeff = (width*width) / Real(4*Iterations);
 
   chi=kronecker;
+
   //  chi = (1-p^2/2N)^N kronecker
   for(int n = 0; n < Iterations; ++n) {
     Laplacian.M(chi,psi);
     chi = chi - coeff*psi;
+    RealD n2 = norm2(chi);
+    chi = chi * (1.0/std::sqrt(n2));
   }
 
   std::cout << " Wuppertal smeared operator is chi = \n" << chi <<std::endl;