test/parallel-fem.cpp

#include "DAPTA.hpp"

#include <set>
#include <fstream>
#include <iostream>
#include <iomanip>
#undef SEEK_SET
#undef SEEK_CUR
#undef SEEK_END
#include "mpi.h"
using namespace DAPTA;

typedef Mesh<Tri3MeshTraits> TriMesh;
typedef MeshParalleliser<Tri3MeshTraits> TriMeshParalleliser;

int main(int argc, char **argv) {
   // Initialize MPI
   MPI_Init(&argc, &argv);
   MPI_Status status;
   
   // Load in the complete mesh
   TriMesh *mesh = new TriMesh();
   mesh->ClearMesh();
   std::ifstream meshfile;
   meshfile.open(argv[1]);
   if(!meshfile) {
      std::cout << "ERROR READING FILE!" << std::endl;
      exit(1);
   }
   mesh->ReadFile(meshfile);
   meshfile.close();
   
   TriMeshParalleliser *mesh_paralleliser = new TriMeshParalleliser(MPI_COMM_WORLD, mesh);
   mesh_paralleliser->initialiseCommunication();
   
   // ********************
   // * BEGIN FEM CODE
   // ********************
   
   int myRank, numProcs;
   MPI_Comm_rank(MPI_COMM_WORLD, &myRank);
   MPI_Comm_size(MPI_COMM_WORLD, &numProcs);
   
   int timestep=1;
   
   while(true) {
      if(myRank == 0) {
         std::cout << "STEP: " << timestep << std::endl;
      }
      
      for(TriMesh::ElementMapIterator e_iter=mesh->ElementMapBegin(); e_iter!=mesh->ElementMapEnd(); ++e_iter) {
      }
      
      timestep++;
      sleep(1);
   }
   
   // ********************
   // * END FEM CODE
   // ********************
   
   // Shut Down MPI
   MPI_Finalize();
   
   return 0;
}

---