Update test_semilinear1d_nonconstant_bounds.py

Author: milzj

tests/algorithms/test_semilinear1d_nonconstant_bounds.py

@@ -110,7 +110,8 @@ def solve_problem(n, n_ref,  u_init=None, maxiter=1000, gtol=1e-15, ftol=-np.inf
         u_vec = u.vector()[:]
         lb_vec = project(lb, U).vector()[:]
         ub_vec = project(ub, U).vector()[:]
-        u_vec = np.clip(u_vec, lb_vec, ub_vec)
+        w_vec = np.clip(u_vec, -beta, beta)
+        u_vec = np.clip(u_vec-w_vec, lb_vec, ub_vec)
         u.vector()[:] = u_vec
 
     y = Function(V)