Commit 04ea3309 authored by Tomáš Oberhuber's avatar Tomáš Oberhuber Committed by Tomáš Oberhuber
Browse files

Subdomain overlaps for the periodic boundary conditions is set only when... 

Subdomain overlaps for the periodic boundary conditions is set only when periodic neighbor is not the same MPI process.
parent e25a7baf

src/TNL/Meshes/DistributedMeshes/SubdomainOverlapsGetter.hpp

+6 −3
Original line number Diff line number Diff line
@@ -35,18 +35,21 @@ getOverlaps( const DistributedMeshType* distributedMesh, 
   TNL_ASSERT_TRUE( distributedMesh != NULL, "" );

   

   const CoordinatesType& subdomainCoordinates = distributedMesh->getSubdomainCoordinates();

   int rank = CommunicatorType::GetRank( CommunicatorType::AllGroup );

   

   for( int i = 0; i < Dimension; i++ )

   {



      CoordinatesType neighborDirection( 0 );

      neighborDirection[ i ] = -1;

      if( subdomainCoordinates[ i ] > 0 )

         lower[ i ] = subdomainOverlapSize;

      else

      else if( distributedMesh->getPeriodicNeighbors()[ Directions::getDirection( neighborDirection ) ] != rank )

         lower[ i ] = periodicBoundariesOverlapSize[ i ];

      

      neighborDirection[ i ] = 1;

      if( subdomainCoordinates[ i ] < distributedMesh->getDomainDecomposition()[ i ] - 1 )

         upper[ i ] = subdomainOverlapSize;

      else

      else if( distributedMesh->getPeriodicNeighbors()[ Directions::getDirection( neighborDirection ) ] != rank )

         upper[ i ] = periodicBoundariesOverlapSize[ i ];

   }

}