Commit 3dcb180a authored by Vít Hanousek's avatar Vít Hanousek
Browse files

Fix explicit updater synchronization and distributed Heat equation initial condition load 

This version should compute same results if distributed or not distributed.
parent 72beb5f8

src/TNL/Problems/HeatEquationProblem_impl.h

+2 −0
Original line number Diff line number Diff line
@@ -149,6 +149,7 @@ setInitialCondition( const Config::ParameterContainer& parameters, 
   if(CommunicatorType::isDistributed())

    {

        Meshes::DistributedMeshes::DistributedGridIO<MeshFunctionType,Meshes::DistributedMeshes::MpiIO> ::load(initialConditionFile, *uPointer );

        uPointer->template Synchronize<CommunicatorType>();

    }

    else

    {
@@ -253,6 +254,7 @@ getExplicitUpdate( const RealType& time, 
   this->explicitUpdater.setBoundaryConditions( this->boundaryConditionPointer ),

   this->explicitUpdater.setRightHandSide( this->rightHandSidePointer ),

   this->explicitUpdater.template update< typename Mesh::Cell, CommType >( time, tau, meshPointer, this->uPointer, fuPointer );



}



template< typename Mesh,

src/TNL/Solvers/PDE/ExplicitUpdater.h

+1 −1
Original line number Diff line number Diff line
@@ -150,7 +150,7 @@ class ExplicitUpdater 
                                             userDataPointer );



         if(CommunicatorType::isDistributed())

            uPointer->template Synchronize<CommunicatorType>();

            fuPointer->template Synchronize<CommunicatorType>();



      }