From 1505205dcd29ee055ea62131b321ae3132e760e5 Mon Sep 17 00:00:00 2001
From: Jukka Aho <ahojukka5@gmail.com>
Date: Mon, 14 Aug 2017 20:50:21 +0300
Subject: [PATCH] use FEMBasis / BasisInfo structure

Calculate basic element info using BasisInfo
---
 src/problems_elasticity.jl | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

diff --git a/src/problems_elasticity.jl b/src/problems_elasticity.jl
index 766b6a1..61c2778 100644
--- a/src/problems_elasticity.jl
+++ b/src/problems_elasticity.jl
@@ -140,12 +140,15 @@ function assemble!{El<:Elasticity3DVolumeElements}(assembly::Assembly,
         Km = zeros(ndofs, ndofs)
         Kg = zeros(ndofs, ndofs)
         f = zeros(ndofs)
+        bi = BasisInfo(El)
 
         for ip in get_integration_points(element)
-            detJ = element(ip, time, Val{:detJ})
+            X = element("geometry", time)
+            eval_basis!(bi, X, ip)
+            detJ = bi.detJ
+            N = bi.N
+            dN = bi.grad
             w = ip.weight*detJ
-            N = element(ip, time)
-            dN = element(ip, time, Val{:Grad})
 
             # kinematics; calculate deformation gradient and strain