Commit 85f0dae7 authored by Vít Hanousek's avatar Vít Hanousek
Browse files

Allow existance of non setup synchronizer object. MeshFunction creates grid... 

Allow existance of non setup synchronizer object. MeshFunction creates grid synchronizer automagicli with binding.
(only if distr grid was chnged)
parent b7fab41b

src/TNL/Functions/MeshFunction.h

+15 −0
Original line number Diff line number Diff line
@@ -13,6 +13,8 @@ 
#include <TNL/Functions/MeshFunctionGnuplotWriter.h>

#include <TNL/Functions/MeshFunctionVTKWriter.h>

#include <TNL/SharedPointer.h>

#include <TNL/Meshes/DistributedGrid.h>

#include <TNL/Meshes/DistributedGridSynchronizer.h>



#pragma once
@@ -156,13 +158,26 @@ class MeshFunction : 
 

      using Object::boundLoad;



#ifdef USE_MPI

      void synchronize(void);    

#endif

 

   protected:



#ifdef USE_MPI

      Meshes::DistributedGridSynchronizer<Meshes::DistributedGrid<MeshType>,ThisType> synchronizer;    

#endif

      

      MeshPointer meshPointer;

      

      VectorType data;

 

      template< typename, typename > friend class MeshFunctionEvaluator;



#ifdef USE_MPI

   private:

      void SetupSynchronizer(Meshes::DistributedGrid<Mesh> *distrgrid);

#endif   

};



} // namespace Functions

src/TNL/Functions/MeshFunction_impl.h

+67 −8
Original line number Diff line number Diff line
@@ -34,8 +34,11 @@ template< typename Mesh, 
          typename Real >

MeshFunction< Mesh, MeshEntityDimension, Real >::

MeshFunction( const MeshPointer& meshPointer )

: meshPointer( meshPointer )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.setSize( getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(), 

               std::cerr << "this->data.getSize() = " << this->data.getSize() << std::endl
@@ -47,8 +50,11 @@ template< typename Mesh, 
          typename Real >

MeshFunction< Mesh, MeshEntityDimension, Real >::

MeshFunction( const ThisType& meshFunction )

: meshPointer( meshFunction.meshPointer )

{

#ifdef USE_MPI

    SetupSynchronizer(meshFunction.meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshFunction.meshPointer;

   this->data.bind( meshFunction.getData() );

}
@@ -60,8 +66,12 @@ MeshFunction< Mesh, MeshEntityDimension, Real >:: 
MeshFunction( const MeshPointer& meshPointer,

              Vector& data,

              const IndexType& offset )

: meshPointer( meshPointer )

//: meshPointer( meshPointer )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.bind( data, offset, getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(), 

               std::cerr << "this->data.getSize() = " << this->data.getSize() << std::endl
@@ -77,8 +87,12 @@ MeshFunction< Mesh, MeshEntityDimension, Real >:: 
MeshFunction( const MeshPointer& meshPointer,

              SharedPointer< Vector >& data,

              const IndexType& offset )

: meshPointer( meshPointer )

//: meshPointer( meshPointer )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.bind( *data, offset, getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(), 

               std::cerr << "this->data.getSize() = " << this->data.getSize() << std::endl
@@ -176,7 +190,10 @@ void 
MeshFunction< Mesh, MeshEntityDimension, Real >::

bind( ThisType& meshFunction )

{

   this->meshPointer = meshFunction.getMeshPointer();

#ifdef USE_MPI

    SetupSynchronizer(meshFunction.meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshFunction.meshPointer;

   this->data.bind( meshFunction.getData() );

}
@@ -190,6 +207,9 @@ bind( const MeshPointer& meshPointer, 
      const Vector& data,

      const IndexType& offset )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.bind( data, offset, getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(),
@@ -207,6 +227,9 @@ bind( const MeshPointer& meshPointer, 
      const SharedPointer< Vector >& data,

      const IndexType& offset )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.bind( *data, offset, getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(),
@@ -222,6 +245,9 @@ void 
MeshFunction< Mesh, MeshEntityDimension, Real >::

setMesh( const MeshPointer& meshPointer )

{

#ifdef USE_MPI

    SetupSynchronizer(meshPointer->GetDistGrid());

#endif

   this->meshPointer=meshPointer;

   this->data.setSize( getMesh().template getEntitiesCount< typename Mesh::template EntityType< MeshEntityDimension > >() );

   TNL_ASSERT( this->data.getSize() == this->getMesh().template getEntitiesCount< typename MeshType::template EntityType< MeshEntityDimension > >(),
@@ -514,6 +540,39 @@ write( const String& fileName, 
   return true;

}



#ifdef USE_MPI

template< typename Mesh,

          int MeshEntityDimension,

          typename Real >

void

MeshFunction< Mesh, MeshEntityDimension, Real >:: 

synchronize(void)

{

    //když je synchronizer korektně natavený

    auto distrgrid = this->getMesh().GetDistGrid();

    if(distrgrid != NULL && distrgrid->isMPIUsed())

    {

        this->synchronizer.Synchronize(*this);

    }

}



template< typename Mesh,

          int MeshEntityDimension,

          typename Real >

void

MeshFunction< Mesh, MeshEntityDimension, Real >:: 

SetupSynchronizer(Meshes::DistributedGrid<Mesh> *distrgrid)

{

    if(distrgrid)//pokud síť kterou se snaží uživatel nastavit je distruibuovaná

    {

        if(this->getMesh().GetDistGrid()!=distrgrid)//pokud má nová síť jinou distribuovanou síť než předchozí

        {

            this->synchronizer.SetDistributedGrid(distrgrid);

        }

    }

}

#endif

 

} // namespace Functions

} // namespace TNL

src/TNL/Meshes/DistributedGrid.h

+115 −66
Original line number Diff line number Diff line
@@ -78,12 +78,12 @@ class DistributedGrid <GridType,1> 


    private : 

        

        GridType GlobalGrid;

        PointType localorigin;

        CoordinatesType localbegin;

        CoordinatesType localsize;

        CoordinatesType localgridsize;

        CoordinatesType overlap;

        PointType spaceSteps;

        

        

        IndexType Dimensions;
@@ -97,18 +97,32 @@ class DistributedGrid <GridType,1> 
        int left;

        int right;



        bool isSet;

        

     

   public:

       DistributedGrid()

       {

            isSet=false;

       };



       //compute everithing 

       DistributedGrid(GridType globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {      

           SetGlobalGrid(globalGrid,overlap,distribution);      

       };



       void SetGlobalGrid(GridType globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {

           isSet=true;



           this->overlap=overlap;

           

           left=-1;

           right=-1;

           

           Dimensions= GridType::getMeshDimension();

           GlobalGrid=globalGrid;

           spaceSteps=globalGrid.getSpaceSteps();

           //Detect MPI and number of process

           mpiInUse=false;

           if(MPI::Is_initialized())
@@ -126,18 +140,18 @@ class DistributedGrid <GridType,1> 
           {

               //Without MPI



               std::cout <<"MEZ MPI"<<std::endl;

               std::cout <<"BEZ MPI"<<std::endl;

               rank=0;

               localorigin=GlobalGrid.getOrigin();

               localsize=GlobalGrid.getDimensions();

               localgridsize=GlobalGrid.getDimensions();

               localorigin=globalGrid.getOrigin();

               localsize=globalGrid.getDimensions();

               localgridsize=globalGrid.getDimensions();

               

               localbegin=CoordinatesType(0);

               return;

           }

           else

           {

               //With MPI

               

               //nearnodes

               if(rank!=0)

                   left=rank-1;
@@ -145,26 +159,26 @@ class DistributedGrid <GridType,1> 
                   right=rank+1;

                  

               //compute local mesh size               

               numberoflarger=GlobalGrid.getDimensions().x()%nproc;

               numberoflarger=globalGrid.getDimensions().x()%nproc;

                 

               localsize.x()=(GlobalGrid.getDimensions().x()/nproc);               

               localsize.x()=(globalGrid.getDimensions().x()/nproc);               

               if(numberoflarger>rank)

                    localsize.x()+=1;                      

                                  

               if(numberoflarger>rank)

                   localorigin.x()=GlobalGrid.getOrigin().x()

                                +(rank*localsize.x()-overlap.x())*GlobalGrid.getSpaceSteps().x();

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(rank*localsize.x()-overlap.x())*globalGrid.getSpaceSteps().x();

               else

                   localorigin.x()=GlobalGrid.getOrigin().x()

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(numberoflarger*(localsize.x()+1)+(rank-numberoflarger)*localsize.x()-overlap.x())

                                *GlobalGrid.getSpaceSteps().x();

                                *globalGrid.getSpaceSteps().x();

              

              localbegin=overlap;

               

               //vlevo neni prekryv

               if(left==-1)

               {

                   localorigin.x()+=overlap.x()*GlobalGrid.getSpaceSteps().x();

                   localorigin.x()+=overlap.x()*globalGrid.getSpaceSteps().x();

                   localbegin.x()=0;

               }
@@ -180,10 +194,11 @@ class DistributedGrid <GridType,1> 
       

       void SetupGrid( GridType& grid)

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by SetupGrid");

           grid.setOrigin(localorigin);

           grid.setDimensions(localgridsize);

           //compute local proporions by sideefect

           grid.setSpaceSteps(GlobalGrid.getSpaceSteps());

           grid.setSpaceSteps(spaceSteps);

           grid.SetDistGrid(this);

       };
@@ -204,11 +219,13 @@ class DistributedGrid <GridType,1> 
       

       int getLeft()

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by getLeft");

           return this->left;

       };

       

       int getRight()

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by getRight");

           return this->right;

       };
@@ -251,7 +268,7 @@ class DistributedGrid <GridType,2> 
    

    private : 

        

        GridType GlobalGrid;

        PointType spaceSteps;

        PointType localorigin;

        CoordinatesType localsize;//velikost gridu zpracovavane danym uzlem bez prekryvu

        CoordinatesType localbegin;//souradnice zacatku zpracovavane vypoctove oblasi
@@ -271,17 +288,31 @@ class DistributedGrid <GridType,2> 
        

        int neighbors[8];



        bool isSet;

     

   public:

       DistributedGrid()

       {

            isSet=false;

       };



       //compute everithing 

       DistributedGrid(GridType globalGrid,CoordinatesType overlap,int *distribution=NULL)

       DistributedGrid(GridType &globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {      

           SetGlobalGrid(globalGrid,overlap,distribution);      

       };



       void SetGlobalGrid(GridType &globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {

           isSet=true;



           this->overlap=overlap;

           

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

                neighbors[i]=-1;

           

           Dimensions= GridType::getMeshDimension();

           GlobalGrid=globalGrid;

           spaceSteps=globalGrid.getSpaceSteps();

           //Detect MPI and number of process

           mpiInUse=false;
@@ -305,9 +336,9 @@ class DistributedGrid <GridType,2> 
               myproccoord[1]=0;

               procsdistr[0]=1;

               procsdistr[1]=1;

               localorigin=GlobalGrid.getOrigin();

               localgridsize=GlobalGrid.getDimensions();

               localsize=GlobalGrid.getDimensions();

               localorigin=globalGrid.getOrigin();

               localgridsize=globalGrid.getDimensions();

               localsize=globalGrid.getDimensions();

               localbegin.x()=0;

               localbegin.y()=0;
@@ -332,11 +363,11 @@ class DistributedGrid <GridType,2> 
               myproccoord[1]=rank/procsdistr[0];        



               //compute local mesh size              

               numberoflarger[0]=GlobalGrid.getDimensions().x()%procsdistr[0];

               numberoflarger[1]=GlobalGrid.getDimensions().y()%procsdistr[1];

               numberoflarger[0]=globalGrid.getDimensions().x()%procsdistr[0];

               numberoflarger[1]=globalGrid.getDimensions().y()%procsdistr[1];



               localsize.x()=(GlobalGrid.getDimensions().x()/procsdistr[0]);

               localsize.y()=(GlobalGrid.getDimensions().y()/procsdistr[1]);

               localsize.x()=(globalGrid.getDimensions().x()/procsdistr[0]);

               localsize.y()=(globalGrid.getDimensions().y()/procsdistr[1]);



               if(numberoflarger[0]>myproccoord[0])

                    localsize.x()+=1;
@@ -344,20 +375,20 @@ class DistributedGrid <GridType,2> 
                   localsize.y()+=1;



               if(numberoflarger[0]>myproccoord[0])

                   localorigin.x()=GlobalGrid.getOrigin().x()

                                +(myproccoord[0]*localsize.x()-overlap.x())*GlobalGrid.getSpaceSteps().x();

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(myproccoord[0]*localsize.x()-overlap.x())*globalGrid.getSpaceSteps().x();

               else

                   localorigin.x()=GlobalGrid.getOrigin().x()

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(numberoflarger[0]*(localsize.x()+1)+(myproccoord[0]-numberoflarger[0])*localsize.x()-overlap.x())

                                *GlobalGrid.getSpaceSteps().x();

                                *globalGrid.getSpaceSteps().x();

               

               if(numberoflarger[1]>myproccoord[1])

                   localorigin.y()=GlobalGrid.getOrigin().y()

                                +(myproccoord[1]*localsize.y()-overlap.y())*GlobalGrid.getSpaceSteps().y();

                   localorigin.y()=globalGrid.getOrigin().y()

                                +(myproccoord[1]*localsize.y()-overlap.y())*globalGrid.getSpaceSteps().y();

               else

                   localorigin.y()=GlobalGrid.getOrigin().y()

                   localorigin.y()=globalGrid.getOrigin().y()

                                +(numberoflarger[1]*(localsize.y()+1)+(myproccoord[1]-numberoflarger[1])*localsize.y()-overlap.y())

                                *GlobalGrid.getSpaceSteps().y();

                                *globalGrid.getSpaceSteps().y();



               //nearnodes

               if(myproccoord[0]>0)
@@ -381,13 +412,13 @@ class DistributedGrid <GridType,2> 


               if(neighbors[Left]==-1)

               {

                    localorigin.x()+=overlap.x()*GlobalGrid.getSpaceSteps().x();

                    localorigin.x()+=overlap.x()*globalGrid.getSpaceSteps().x();

                    localbegin.x()=0;

               }



               if(neighbors[Up]==-1)

               {

                   localorigin.y()+=overlap.y()*GlobalGrid.getSpaceSteps().y();

                   localorigin.y()+=overlap.y()*globalGrid.getSpaceSteps().y();

                   localbegin.y()=0;

               }
@@ -410,10 +441,11 @@ class DistributedGrid <GridType,2> 
       

       void SetupGrid( GridType& grid)

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by SetupGrid");

           grid.setOrigin(localorigin);

           grid.setDimensions(localgridsize);

           //compute local proporions by sideefect

           grid.setSpaceSteps(GlobalGrid.getSpaceSteps());

           grid.setSpaceSteps(spaceSteps);

           grid.SetDistGrid(this);

       };
@@ -439,6 +471,7 @@ class DistributedGrid <GridType,2> 
       

       int * getNeighbors()

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by getNeighbors");

           return this->neighbors;

       }
@@ -488,7 +521,7 @@ class DistributedGrid <GridType,3> 


    private : 

        

        GridType GlobalGrid;

        PointType spaceSteps;

        PointType localorigin;

        CoordinatesType localsize;

        CoordinatesType localgridsize;
@@ -508,10 +541,24 @@ class DistributedGrid <GridType,3> 
        

        int neighbors[26];



        bool isSet;

     

   public:

       DistributedGrid()

       {

            isSet=false;

       };



       //compute everithing 

       DistributedGrid(GridType globalGrid,CoordinatesType overlap,int *distribution=NULL)

       DistributedGrid(GridType &globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {      

           SetGlobalGrid(globalGrid,overlap,distribution);      

       };



       void SetGlobalGrid(GridType &globalGrid, CoordinatesType overlap, int *distribution=NULL)

       {



           isSet=true;           



           this->overlap=overlap;
@@ -519,7 +566,7 @@ class DistributedGrid <GridType,3> 
                neighbors[i]=-1;

           

           Dimensions= GridType::getMeshDimension();

           GlobalGrid=globalGrid;

           spaceSteps=globalGrid.getSpaceSteps();

           //Detect MPI and number of process

           mpiInUse=false;
@@ -547,8 +594,8 @@ class DistributedGrid <GridType,3> 
               procsdistr[1]=1;

               procsdistr[2]=1;               

               

               localorigin=GlobalGrid.getOrigin();

               localsize=GlobalGrid.getDimensions();

               localorigin=globalGrid.getOrigin();

               localsize=globalGrid.getDimensions();

               localgridsize=localsize;

               return;

           }
@@ -574,13 +621,13 @@ class DistributedGrid <GridType,3> 
               myproccoord[0]=(rank%(procsdistr[0]*procsdistr[1]))%procsdistr[0];



               //compute local mesh size               

               numberoflarger[0]=GlobalGrid.getDimensions().x()%procsdistr[0];

               numberoflarger[1]=GlobalGrid.getDimensions().y()%procsdistr[1];

               numberoflarger[2]=GlobalGrid.getDimensions().z()%procsdistr[2];

               numberoflarger[0]=globalGrid.getDimensions().x()%procsdistr[0];

               numberoflarger[1]=globalGrid.getDimensions().y()%procsdistr[1];

               numberoflarger[2]=globalGrid.getDimensions().z()%procsdistr[2];

                 

               localsize.x()=(GlobalGrid.getDimensions().x()/procsdistr[0]);

               localsize.y()=(GlobalGrid.getDimensions().y()/procsdistr[1]);

               localsize.z()=(GlobalGrid.getDimensions().z()/procsdistr[2]);

               localsize.x()=(globalGrid.getDimensions().x()/procsdistr[0]);

               localsize.y()=(globalGrid.getDimensions().y()/procsdistr[1]);

               localsize.z()=(globalGrid.getDimensions().z()/procsdistr[2]);

               

               if(numberoflarger[0]>myproccoord[0])

                    localsize.x()+=1;
@@ -590,28 +637,28 @@ class DistributedGrid <GridType,3> 
                   localsize.z()+=1;

                                  

               if(numberoflarger[0]>myproccoord[0])

                   localorigin.x()=GlobalGrid.getOrigin().x()

                                +(myproccoord[0]*localsize.x()-overlap.x())*GlobalGrid.getSpaceSteps().x();

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(myproccoord[0]*localsize.x()-overlap.x())*globalGrid.getSpaceSteps().x();

               else

                   localorigin.x()=GlobalGrid.getOrigin().x()

                   localorigin.x()=globalGrid.getOrigin().x()

                                +(numberoflarger[0]*(localsize.x()+1)+(myproccoord[0]-numberoflarger[0])*localsize.x()-overlap.x())

                                *GlobalGrid.getSpaceSteps().x();

                                *globalGrid.getSpaceSteps().x();

               

               if(numberoflarger[1]>myproccoord[1])

                   localorigin.y()=GlobalGrid.getOrigin().y()

                                +(myproccoord[1]*localsize.y()-overlap.y())*GlobalGrid.getSpaceSteps().y();

                   localorigin.y()=globalGrid.getOrigin().y()

                                +(myproccoord[1]*localsize.y()-overlap.y())*globalGrid.getSpaceSteps().y();

               else

                   localorigin.y()=GlobalGrid.getOrigin().y()

                   localorigin.y()=globalGrid.getOrigin().y()

                                +(numberoflarger[1]*(localsize.y()+1)+(myproccoord[1]-numberoflarger[1])*localsize.y()-overlap.y())

                                *GlobalGrid.getSpaceSteps().y();

                                *globalGrid.getSpaceSteps().y();



               if(numberoflarger[2]>myproccoord[2])

                   localorigin.z()=GlobalGrid.getOrigin().z()

                                +(myproccoord[2]*localsize.z()-overlap.z())*GlobalGrid.getSpaceSteps().z();

                   localorigin.z()=globalGrid.getOrigin().z()

                                +(myproccoord[2]*localsize.z()-overlap.z())*globalGrid.getSpaceSteps().z();

               else

                   localorigin.z()=GlobalGrid.getOrigin().z()

                   localorigin.z()=globalGrid.getOrigin().z()

                                +(numberoflarger[2]*(localsize.z()+1)+(myproccoord[2]-numberoflarger[2])*localsize.z()-overlap.z())

                                *GlobalGrid.getSpaceSteps().z();

                                *globalGrid.getSpaceSteps().z();

               

               //nearnodes

               //X Y Z
@@ -678,17 +725,17 @@ class DistributedGrid <GridType,3> 
               

               if(neighbors[West]==-1)

               {

                   localorigin.x()+=overlap.x()*GlobalGrid.getSpaceSteps().x();

                   localorigin.x()+=overlap.x()*globalGrid.getSpaceSteps().x();

                   localbegin.x()=0;

               }

               if(neighbors[Nord]==-1)

               {

                   localorigin.y()+=overlap.y()*GlobalGrid.getSpaceSteps().y();

                   localorigin.y()+=overlap.y()*globalGrid.getSpaceSteps().y();

                   localbegin.y()=0;

               }

               if(neighbors[Bottom]==-1)

               {

                   localorigin.z()+=overlap.z()*GlobalGrid.getSpaceSteps().z();

                   localorigin.z()+=overlap.z()*globalGrid.getSpaceSteps().z();

                   localbegin.z()=0;

               }
@@ -714,10 +761,11 @@ class DistributedGrid <GridType,3> 
       

       void SetupGrid( GridType& grid)

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by SetupGrid");

           grid.setOrigin(localorigin);

           grid.setDimensions(localgridsize);

           //compute local proporions by sideefect

           grid.setSpaceSteps(GlobalGrid.getSpaceSteps());

           grid.setSpaceSteps(spaceSteps);

           grid.SetDistGrid(this);

       };
@@ -743,6 +791,7 @@ class DistributedGrid <GridType,3> 
       

       int * getNeighbors()

       {

           TNL_ASSERT_TRUE(isSet,"DistributedGrid is not set, but used by getNeighbors");

           return this->neighbors;

       }

src/TNL/Meshes/DistributedGridSynchronizer.h

+120 −42

File changed.

Preview size limit exceeded, changes collapsed.

src/UnitTests/Mpi/DistributedGridIOTest.cpp

+4 −3
Original line number Diff line number Diff line
@@ -330,7 +330,7 @@ class TestDistributedGridIO{ 
        CoordinatesType overlap;

        overlap.setValue(1);

        DistributedGridType distrgrid(globalGrid,overlap, distr);

        DistributedGridSynchronizer<DistributedGridType,MeshFunctionType> synchronizer(&distrgrid);

        //DistributedGridSynchronizer<DistributedGridType,MeshFunctionType> synchronizer(&distrgrid);



        //Crete "distributedgrid driven" grid filed by load

        SharedPointer<MeshType> loadGridptr;
@@ -346,7 +346,7 @@ class TestDistributedGridIO{ 
        DistributedGridIO<MeshFunctionType> ::load(file, *loadMeshFunctionptr );

        file.close();



        synchronizer.Synchronize(*loadMeshFunctionptr);//need synchronization for overlaps to be filled corectly in loadDof

        loadMeshFunctionptr->synchronize();//need synchronization for overlaps to be filled corectly in loadDof





        //Crete "distributedgrid driven" grid filed by evaluated linear function
@@ -359,7 +359,8 @@ class TestDistributedGridIO{ 
        meshFunctionptr->bind(gridptr,dof);

        

        linearFunctionEvaluator.evaluateAllEntities(meshFunctionptr , linearFunctionPtr);

        synchronizer.Synchronize(*meshFunctionptr);//need synchronization for overlaps to be filled corectly in dof

        //synchronizer.Synchronize(*meshFunctionptr);//need synchronization for overlaps to be filled corectly in dof

        meshFunctionptr->synchronize();



        for(int i=0;i<localDof.getSize();i++)

        {