diff --git a/src/ASM/SplineField2D.C b/src/ASM/SplineField2D.C
index 9ff4ede4..cf5f0579 100644
--- a/src/ASM/SplineField2D.C
+++ b/src/ASM/SplineField2D.C
@@ -118,6 +118,7 @@ bool SplineField2D::gradFE (const FiniteElement& fe, Vector& grad) const
     }
     dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
 
+    ip.clear();
     ASMs2D::scatterInd(basis->numCoefs_u(),basis->numCoefs_v(),
 		       basis->order_u(),basis->order_v(),
 		       spline.left_idx,ip);
diff --git a/src/ASM/SplineField3D.C b/src/ASM/SplineField3D.C
index 02c1c42f..9c0d9700 100644
--- a/src/ASM/SplineField3D.C
+++ b/src/ASM/SplineField3D.C
@@ -123,6 +123,7 @@ bool SplineField3D::gradFE (const FiniteElement& fe, Vector& grad) const
     }
     dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
 
+    ip.clear();
     ASMs3D::scatterInd(basis->numCoefs(0),basis->numCoefs(1),basis->numCoefs(2),
 		       basis->order(0),basis->order(1),basis->order(2),
 		       spline.left_idx,ip);
diff --git a/src/ASM/SplineFields2D.C b/src/ASM/SplineFields2D.C
index 30535bbe..b63fa105 100644
--- a/src/ASM/SplineFields2D.C
+++ b/src/ASM/SplineFields2D.C
@@ -125,6 +125,7 @@ bool SplineFields2D::gradFE (const FiniteElement& fe, Matrix& grad) const
     }
     dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
 
+    ip.clear();
     ASMs2D::scatterInd(basis->numCoefs_u(),basis->numCoefs_v(),
 		       basis->order_u(),basis->order_v(),
 		       spline.left_idx,ip);
diff --git a/src/ASM/SplineFields3D.C b/src/ASM/SplineFields3D.C
index 1d17824b..c17cae64 100644
--- a/src/ASM/SplineFields3D.C
+++ b/src/ASM/SplineFields3D.C
@@ -130,6 +130,7 @@ bool SplineFields3D::gradFE(const FiniteElement& fe, Matrix& grad) const
     }
     dNdX.multiply(dNdu,Jac); // dNdX = dNdu * Jac
 
+    ip.clear();
     ASMs3D::scatterInd(basis->numCoefs(0),basis->numCoefs(1),basis->numCoefs(2),
 		       basis->order(0),basis->order(1),basis->order(2),
 		       spline.left_idx,ip);