From 2cdf34bfaaa4a79cb61f8fcce0115f06041b269e Mon Sep 17 00:00:00 2001
From: Arne Morten Kvarving <arne.morten.kvarving@sintef.no>
Date: Mon, 13 Mar 2017 10:20:39 +0100
Subject: [PATCH] fixed: this is 3D so we have nGauss^3 integration points per
 element

---
 src/ASM/LR/ASMu3D.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/ASM/LR/ASMu3D.C b/src/ASM/LR/ASMu3D.C
index 6624e765..48ab8cde 100644
--- a/src/ASM/LR/ASMu3D.C
+++ b/src/ASM/LR/ASMu3D.C
@@ -1119,7 +1119,7 @@ bool ASMu3D::integrate (Integrand& integrand,
 
       // --- Integration loop over all Gauss points in each direction --------
 
-      fe.iGP = iEl*nGauss*nGauss; // Global integration point counter
+      fe.iGP = iEl*nGauss*nGauss*nGauss; // Global integration point counter
 
       Matrix B(p1*p2*p3, 4); // Bezier evaluation points and derivatives
       int ig = 1;