From a0e6e030e3d41cd1c20aedc12b36ffdb522fca5c Mon Sep 17 00:00:00 2001 From: Kristoffer Carlsson Date: Fri, 2 Nov 2018 13:47:57 -0400 Subject: [PATCH] improve performance in elasticity by moving out computation of X from intergration point loop --- src/problems_elasticity.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/problems_elasticity.jl b/src/problems_elasticity.jl index faff70b..fc7f6a6 100644 --- a/src/problems_elasticity.jl +++ b/src/problems_elasticity.jl @@ -130,13 +130,14 @@ function assemble!(assembly::Assembly, for element in elements u = element("displacement", time) + X = element("geometry", time) + fill!(Km, 0.0) fill!(Kg, 0.0) fill!(f_int, 0.0) fill!(f_ext, 0.0) for ip in get_integration_points(element) - X = element("geometry", time) eval_basis!(bi, X, ip) w = ip.weight*bi.detJ N = bi.N