refactoring

mluq/src/pdesolvers/EllipticPDESolver.hpp

@@ -31,8 +31,10 @@ protected:
 public:
   EllipticPDESolver(IStochasticEllipticAssemble *assemble,
                     const string &quantity, const string &costMeasure) :
-    PDESolver(), assemble(assemble), quantity(quantity), costMeasure(costMeasure),
-    newton(std::make_unique<Newton>()) {
+    PDESolver(), assemble(assemble), quantity(quantity), costMeasure(costMeasure) {
+
+    newton = std::make_unique<Newton>();
+
     if (verbose > 0) {
       mout.PrintInfo("EllipticPDESolver", verbose,
                      PrintInfoEntry("Quantity", quantity),

mluq/src/pdesolvers/TransportPDESolver.hpp

@@ -28,13 +28,15 @@ protected:
 public:
   TransportPDESolver(IStochasticLinearTransportAssemble *assemble,
                      const string &quantity, const string &costMeasure) :
-    PDESolver(), assemble(assemble), quantity(quantity), costMeasure(costMeasure),
-    timeInt(TimeIntegratorCreator(IMPLICIT_MIDPOINT).
+    PDESolver(), assemble(assemble), quantity(quantity), costMeasure(costMeasure) {
+
+    timeInt = TimeIntegratorCreator(IMPLICIT_MIDPOINT).
       WithLinearSolver(new GMRES(GetPC("PointBlockJacobi_dG"))).
       WithLinearSolver(new GMRES(GetPC("PointBlockJacobi_dG"))).
       WithLinearSolver(new GMRES(GetPC("PointBlockJacobi_dG"))).
       WithLinearSolver(new GMRES(GetPC("PointBlockJacobi_dG"))).
-      CreateUnique()) {
+      CreateUnique();
+
     if (verbose)
       mout.PrintInfo("TransportPDESolver", verbose,
                      PrintInfoEntry("Quantity", quantity),