From 0d4b08f576160ef4dc259253408b541c0aa23f8a Mon Sep 17 00:00:00 2001 From: Jukka Aho Date: Sat, 5 Nov 2016 18:27:39 +0200 Subject: [PATCH] fixed bug in integration rule of tet10 --- src/integrate.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/integrate.jl b/src/integrate.jl index 9663e92..8cd8ecd 100644 --- a/src/integrate.jl +++ b/src/integrate.jl @@ -205,7 +205,7 @@ function get_integration_points(element::TetrahedralElement, ::Type{Val{4}}) w2 = (2665.0 - 14.0*sqrt(15.0))/226800.0 w3 = (2665.0 + 14.0*sqrt(15.0))/226800.0 w4 = 5.0/567.0 - weights = [w1, w2, w2, w2, w2, w3, w3, w3, w3, w4, w4, w4, w4, w4] + weights = [w1, w2, w2, w2, w2, w3, w3, w3, w3, w4, w4, w4, w4, w4, w4] points = Vector{Float64}[ [a, a, a], [b1, b1, b1],