demo_drivers/mpi/distribution/hanging_node_reconciliation/adaptive_driven_cavity.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableDrivenCavityProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char **argv)
	Driver for RefineableDrivenCavity test problem. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver for RefineableDrivenCavity test problem.

{
 
#ifdef OOMPH_HAS_MPI
 
 // Initialise MPI
 MPI_Helpers::init(argc,argv);
 
#endif
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Set output directory
 DocInfo doc_info;
 doc_info.set_directory("RESLT");
 
 // Set max. number of black-box adaptation
 unsigned max_adapt=2;
 
 // Solve problem with Taylor Hood elements
 //---------------------------------------
 {
  //Build problem
  RefineableDrivenCavityProblem<RefineableQTaylorHoodElement<2> > problem;
 
  
#ifdef OOMPH_HAS_MPI
  
  // DocInfo object specifies directory in which we document
  // the distribution
  DocInfo mesh_doc_info;
  mesh_doc_info.set_directory("RESLT_TH_MESH");
  
  // Create storage for pre-determined partitioning
  const unsigned n_element=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_element);
  
  // Eqidistribute mesh to create non-synchronised hanging node
  // where needed
  for (unsigned e=0;e<8;e++)
   {
    element_partition[e]=0;
    element_partition[e+8]=1;
   }
 
  // Now perform the distribution and document it
  bool report_stats=true;
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  
#endif
  
  // solve with adaptation
  problem.newton_solve(max_adapt);
  
  //Output solution
  problem.doc_solution(doc_info);
  
#ifdef OOMPH_HAS_MPI
  
  mesh_doc_info.number()=1;
  problem.mesh_pt()->doc_mesh_distribution(mesh_doc_info);
#endif
 
 
} // end of Taylor Hood elements
 
 
 // Solve problem with Crouzeix Raviart elements
 //--------------------------------------------
 {
  // Build problem
  RefineableDrivenCavityProblem<RefineableQCrouzeixRaviartElement<2> > problem;
  
#ifdef OOMPH_HAS_MPI
  
  DocInfo mesh_doc_info;
  mesh_doc_info.set_directory("RESLT_CR_MESH");
  mesh_doc_info.number()=0;
  
  
  // Eqidistribute mesh to create non-synchronised hanging node
  // where needed
  unsigned n_el=problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_el);
  for (unsigned e=0;e<8;e++)
   {
    element_partition[e]=0;
    element_partition[e+8]=1;
   }
  
  // Now perform the distribution
  bool report_stats=true;
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  
#endif
  
  // Re-solve with adaptation
  problem.newton_solve(max_adapt);
  
  // change doc_info number
  doc_info.number()=1;
  
  //Output solution
  problem.doc_solution(doc_info);
  
#ifdef OOMPH_HAS_MPI
  mesh_doc_info.number()=1;
  problem.mesh_pt()->doc_mesh_distribution(mesh_doc_info);
#endif
 
 
 } // end of Crouzeix Raviart elements
 
 
// Finalise MPI
#ifdef OOMPH_HAS_MPI
 
 MPI_Helpers::finalize();
 
#endif
 
} // end_of_main

References e(), oomph::DocInfo::number(), problem, oomph::DocInfo::set_directory(), and Flag_definition::setup().