Base the energyNorm on the stiffnessMatrix

src/one-body-sample.cc

@@ -135,7 +135,8 @@ int main() {
     BlockGSStep<OperatorType, VectorType> blockGSStep(stiffnessMatrix, u, f);
     blockGSStep.ignoreNodes_ = &ignoreNodes;
 
-    EnergyNorm<OperatorType, VectorType> energyNorm(blockGSStep);
+    // FIXME: Does this make any sense?
+    EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix);
     LoopSolver<VectorType> solver(&blockGSStep, solver_maxIterations,
                                   solver_tolerance, &energyNorm, Solver::FULL);