Add updateTerm

src/one-body-sample.cc

@@ -193,6 +193,10 @@ int main() {
           leafView, p1Basis, neumannNodes, neumannTerm);
       f += neumannTerm;
 
+      VectorType updateTerm(grid.size(grid.maxLevel(), dim));
+      stiffnessMatrix.mv(u1, updateTerm);
+      f += updateTerm;
+
       // {{{ Assemble terms for the nonlinearity
       // TODO: Random value
       std::vector<double> normalStress;