diff --git a/Solver/helmholtzTerm.h b/Solver/helmholtzTerm.h
index f79f3e5b65239e753f234bc5e682570dfe24b6c4..45f0b9961260f2fdc93570d1b39ffe08f7c624ee 100644
--- a/Solver/helmholtzTerm.h
+++ b/Solver/helmholtzTerm.h
@@ -50,9 +50,9 @@ class helmholtzTerm : public femTerm<scalar> {
   {
     MElement *e = se->getMeshElement();
     // compute integration rule
-    //const int integrationOrder = (_a) ? 2 * e->getPolynomialOrder() : 
-    //  2 * (e->getPolynomialOrder() - 1);
-        const int integrationOrder = 2 * e->getPolynomialOrder() + 1; 
+    const int integrationOrder = (_a) ? 2 * e->getPolynomialOrder() : 
+      2 * (e->getPolynomialOrder() - 1);
+   //    const int integrationOrder = 2 * e->getPolynomialOrder() + 1; 
 
     int npts; IntPt *GP;
     e->getIntegrationPoints(integrationOrder, &npts, &GP);