diff --git a/src/electrophysiology/assemble/IElphyAssemble.cpp b/src/electrophysiology/assemble/IElphyAssemble.cpp index 055e077119aff7770a175d94b7a32a63e95e2fa1..69ad210c5e6d23edc4acb1dfa1f00c2f20b5c9e4 100644 --- a/src/electrophysiology/assemble/IElphyAssemble.cpp +++ b/src/electrophysiology/assemble/IElphyAssemble.cpp @@ -162,7 +162,7 @@ void IElphyAssemble::PlotExcitation() const { void IElphyAssemble::PlotActivation() const { if (problem.UpdateActivationTime == false) return; - //writeActivationTimeFile(); + writeActivationTimeFile(); auto &plot = mpp::plot("ActivationTime"); plot.AddData("Activation", *activationTime); plot.PlotFile("ActivationTime"); @@ -185,14 +185,15 @@ void IElphyAssemble::writeActivationTimeFile() const { std::cout << "failed to open " << filename << '\n'; } else { for (int i = 0; i < PPM->Size(); i++) { - for (row r = (*activationTime).rows(); r != (*activationTime).rows_end(); ++r) { - Point p; - double actiTime; - B.Receive(i) >> p >> actiTime ; - file_out << p << " " <<actiTime << endl; - } + while(B.Receive(i).size()<B.Receive(i).Size()) { + Point p; + double actiTime; + B.Receive(i) >> p >> actiTime; + file_out << p << " " << actiTime << endl; + } } } + file_out.close(); } }