diff --git a/BEM_run.m b/BEM_run.m
index cebfefe41925bf4fae44cc61f726e43c3959f967..5f131cd84bb187371bf5d1ee29d9a2a4d2457aed 100644
--- a/BEM_run.m
+++ b/BEM_run.m
@@ -96,7 +96,8 @@ for i = 1:Nconv
     vLR = vL + vR;
     
     % Half-wave rectify and square
-    vLR = HalfSquare(vLR);
+    %vLR = HalfSquare(vLR);
+    vLR = max(vLR,0).^2;
     
     % Reshape
     vLR = permute(vLR, [2 3 1]);