encap: Fix recv/reevaluate for PFASST.

Signed-off-by: Matthew Emmett <[email protected]>

Author: memmett

include/pfasst/encap/encap_sweeper.hpp

@@ -118,7 +118,6 @@ namespace pfasst
         virtual void spread() override
         {
           for (size_t m = 1; m < this->quadrature->get_num_nodes(); m++) {
-            //            this->get_state(m)->copy(this->start_state);
             this->state[m]->copy(this->state[0]);
           }
         }
@@ -274,6 +273,9 @@ namespace pfasst
         virtual void recv(ICommunicator* comm, int tag, bool blocking) override
         {
           this->start_state->recv(comm, tag, blocking);
+          if (this->quadrature->left_is_node()) {
+            this->state[0]->copy(this->start_state);
+          }
         }
 
         virtual void broadcast(ICommunicator* comm) override

include/pfasst/encap/imex_sweeper.hpp

@@ -168,11 +168,10 @@ namespace pfasst
           time dt = this->get_controller()->get_time_step();
           if (initial_only) {
             if (this->quadrature->left_is_node()) {
-            this->f_expl_eval(this->fs_expl[0], this->start_state, t0);
-            this->f_impl_eval(this->fs_impl[0], this->start_state, t0);
-
+              this->f_expl_eval(this->fs_expl[0], this->state[0], t0);
+              this->f_impl_eval(this->fs_impl[0], this->state[0], t0);
             } else {
-
+              throw NotImplementedYet("reevaluate");
             }
           } else {
             for (size_t m = 0; m < this->quadrature->get_num_nodes(); m++) {