spr: save the reduction result

duhhhhh

Author: cwsmith

src/MeshField_SPR_ErrorEstimator.hpp

@@ -24,7 +24,7 @@ class SInt : public MeshField::Integrator {
   void reset() { r = 0; }
   MeshField::Real r;
   Omega_h::CommPtr comm;
-  void parallelReduce() { comm->allreduce(r, OMEGA_H_SUM); }
+  void parallelReduce() { r = comm->allreduce(r, OMEGA_H_SUM); }
 };
 
 template <typename ShapeField> class Estimation {