mpi/distribution/bifurcation_tracking/track_pitch.cc File Reference

#include <iostream>
#include <fstream>
#include <cstdio>
#include "generic.h"
#include "axisym_navier_stokes.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	CylinderMesh< ELEMENT >
class	RotatingProblem< ELEMENT >

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

{
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 
 
 //Construct the problem, you can use either element types
 RotatingProblem<AxisymmetricQCrouzeixRaviartElement> problem(20,20,10,10);
 
 //Set up a dummy partition
 
#ifdef OOMPH_HAS_MPI
 //Set up a dummy partition
 unsigned n_element = problem.mesh_pt()->nelement();
 Vector<unsigned> element_partition(n_element);
 for(unsigned e=0;e<n_element/2;e++) {element_partition[e]=0;}
 for(unsigned e=n_element/2;e<n_element;e++) {element_partition[e]=1;}
 
 //DocInfo mesh_doc_info;
 //bool report_stats=true;
 //mesh_doc_info.set_directory("RESLT_MESH");
 problem.distribute(element_partition);//,mesh_doc_info,report_stats);
 //problem.check_halo_schemes(mesh_doc_info);
 //problem.distribute();
#endif
 
 
 //Solve the problem :)
 problem.solve_system();
 
 
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::finalize();
#endif
}

References e(), and problem.