Commit dac34ebc authored by Tomáš Oberhuber's avatar Tomáš Oberhuber
Browse files

Refactoring the explicit solver updater. 

Adding OpenMP support.
parent ed7a9e18

.gitignore

+1 −0
Original line number Diff line number Diff line
@@ -41,3 +41,4 @@ 
# /tools/src/

/tools/src/Makefile.in

/Testing

/CMakeLists.txt.user

CMakeLists.txt

+1 −0
Original line number Diff line number Diff line
@@ -99,6 +99,7 @@ endif( NOT WITH_CUDA STREQUAL "no" ) 
#

find_package( OpenMP ) 

if( OPENMP_FOUND )

   AddCompilerFlag( "-DHAVE_OPENMP -fopenmp" )

# TODO: finish this

endif()

examples/heat-equation/heatEquationSolver_impl.h

+0 −1
Original line number Diff line number Diff line
@@ -159,7 +159,6 @@ void heatEquationSolver< Mesh,Diffusion,BoundaryCondition,RightHandSide > 


   this->bindDofs( mesh, _u );

   explicitUpdater.template update< Mesh::Dimensions >( time,

                                                        tau,

                                                        mesh,

                                                        this->differentialOperator,

                                                        this->boundaryCondition,

src/implementation/mesh/tnlTraversal_Grid1D_impl.h

+3 −0
Original line number Diff line number Diff line
@@ -63,6 +63,9 @@ processEntities( const GridType& grid, 
    */

   CoordinatesType coordinates;

   const IndexType& xSize = grid.getDimensions().x();

#ifdef HAVE_OPENMP

//#pragma omp parallel for

#endif

   for( coordinates.x() = 0; coordinates.x() <= xSize; coordinates.x() ++ )

   {

      const IndexType index = grid.getVertexIndex( coordinates );

src/implementation/mesh/tnlTraversal_Grid2D_impl.h

+3 −0
Original line number Diff line number Diff line
@@ -37,6 +37,9 @@ processEntities( const GridType& grid, 
   CoordinatesType coordinates;

   const IndexType& xSize = grid.getDimensions().x();

   const IndexType& ySize = grid.getDimensions().y();

#ifdef HAVE_OPENMP

//#pragma omp parallel for

#endif

   for( coordinates.y() = 0; coordinates.y() < ySize; coordinates.y() ++ )

      for( coordinates.x() = 0; coordinates.x() < xSize; coordinates.x() ++ )

      {