Commit 9f44ea1d authored by Matouš Fencl's avatar Matouš Fencl Committed by Tomáš Oberhuber
Browse files

Change of int to IndexType and preparations for OpenMPI.

parent 552e90c4

src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem.h

+2 −0
Original line number Diff line number Diff line
@@ -25,6 +25,8 @@ 
#include <solvers/pde/tnlLinearSystemAssembler.h>

#include <functions/tnlMeshFunction.h>



#include <TNL/Meshes/DistributedMeshes/DistributedGridIO.h>



template< typename Mesh,

		    typename DifferentialOperator,

		    typename BoundaryCondition,

src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem_impl.h

+16 −4
Original line number Diff line number Diff line
@@ -123,11 +123,23 @@ setInitialCondition( const Config::ParameterContainer& parameters, 
{

  this->bindDofs( mesh, dofs );

  const String& initialConditionFile = parameters.getParameter< String >( "initial-condition" );

  if(CommunicatorType::isDistributed())

  {

    std::cout<<"Nodes Distribution: " << uPointer->getMesh().getDistributedMesh()->printProcessDistr() << std::endl;

    if(distributedIOType==Meshes::DistributedMeshes::MpiIO)

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

    if(distributedIOType==Meshes::DistributedMeshes::LocalCopy)

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

    uPointer->template synchronize<CommunicatorType>();

  }

  else

  {

    if( ! this->solution.boundLoad( initialConditionFile ) )

    {

      std::cerr << "I am not able to load the initial condition from the file " << initialConditionFile << "." <<std::endl;

      return false;

    }

  }

   return true;

}