diff --git a/src/operator.jl b/src/operator.jl index ea428e8ea6098d47ef335867dff539f83b139fa4..c8adc04c3e086f49c84c9837469c2defc6a709ae 100644 --- a/src/operator.jl +++ b/src/operator.jl @@ -34,6 +34,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] +#quadrature() = SA[1/3, 1/3, 1/3], SA[0. 0. 1.; 0. 1. 0.] assemble(op::Operator) = assemble(op.space,