diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index 15ed967965faebf3d0cff17c39e1f02d2b2970a8..eaf97abfcda7b59d52243b4825f2f6e707413138 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -738,7 +738,7 @@ getStrain(const std::vector<Configuration>& sol,
             FieldVector<double,blocksize> localStrain = getStrain(sol, it, quad[pt].position());
             
             // Sum it all up
-            strain.axpy(weight, localStrain);
+            strain[elementIdx].axpy(weight, localStrain);
 
         }