rSCINTROPARALLEL/grid.hh
7a8b45515eb6implementation/c++

grid.hh
View Options

			#ifndef GRID_HH
			#define GRID_HH

			/* -------------------------------------------------------------------------- */
			#include <vector>
			#include <array>
			#include <mpi.h>
			/* -------------------------------------------------------------------------- */

			class Grid {
			public:
			Grid(int m = 0, int n = 0);

			/// access the value [i][j] of the grid
			inline float & operator()(int i, int j) { return m_storage[i * m_n + j]; }
			inline const float & operator()(int i, int j) const {
			return m_storage[i * m_n + j];
			}

			void resize(int m, int n);

			/// set the grid to 0
			void clear();

			int m() const;
			int n() const;

			void initializeHaloCommunications(MPI_Comm & m_communicator);
			void startSynchronization();
			void waitReceives();
			void waitSends();

			private:
			int m_m, m_n;
			std::vector<float> m_storage;

			/// local neighbors
			int m_north_prank, m_south_prank;

			/// proc rank
			int m_prank;

			/// nb process
			int m_psize;

			/// request for asynchronous communications
			std::array<MPI_Request, 4> m_requests;
			};

			#endif /* GRID_HH */