From 6d2439bcb244341af7947fffc48af52787425325 Mon Sep 17 00:00:00 2001
From: Stephan Hilb <stephan@ecshi.net>
Date: Fri, 25 Feb 2022 23:17:26 +0100
Subject: [PATCH] fix quadrature scaling

---
 src/operator.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/operator.jl b/src/operator.jl
index 634b2e1..15c76f7 100644
--- a/src/operator.jl
+++ b/src/operator.jl
@@ -35,7 +35,7 @@ a(op::Poisson, xloc, u, du, v, dv) = dot(du, dv)
 # degree 2 quadrature on triangle
 quadrature() = SA[1/6, 1/6, 1/6], SA[1/6 4/6 1/6; 1/6 1/6 4/6]
 # degree 1 quadrature on triangle using corner points
-quadrature_lumped() = SA[1/3, 1/3, 1/3], SA[0. 0. 1.; 0. 1. 0.]
+quadrature_lumped() = SA[1/6, 1/6, 1/6], SA[0. 0. 1.; 0. 1. 0.]
 
 
 assemble(op::Operator) = assemble(op.space,
-- 
GitLab