mesh_from_triangle_navier_stokes.cc File Reference

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

Classes
class	FlowPastBoxProblem< ELEMENT >
	Flow past a box in a channel. More...

Namespaces
	Global_Physical_Variables
	Global variables.

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

Function Documentation

main()

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

Driver for FlowPastBox test problem.

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////

{
 
// Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 
 // Check number of command line arguments: Need exactly two.
 if (argc!=4)
  {
   std::string error_message =
    "Wrong number of command line arguments.\n";
   error_message +=
    "Must specify the following file names  \n";
   error_message += 
    "filename.node then filename.ele then filename.poly\n";
   
   throw OomphLibError(error_message,
                       OOMPH_CURRENT_FUNCTION,
                       OOMPH_EXCEPTION_LOCATION);
  }
 
 
 // Convert arguments to strings that specify the input file names
 string node_file_name(argv[1]);
 string element_file_name(argv[2]);
 string poly_file_name(argv[3]);
 
 
 // Set up doc info
 // ---------------
 
 // Label for output
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 // Step number
 doc_info.number()=0;
 
#ifdef ADAPT
 const unsigned max_adapt = 3;
#endif
 
// Doing TTaylorHoodElements
 {
#ifdef ADAPT
  // Build the problem with TTaylorHoodElements
  FlowPastBoxProblem<ProjectableTaylorHoodElement<TTaylorHoodElement<2> > >
   problem(node_file_name, element_file_name, poly_file_name);
#else
  // Build the problem with TTaylorHoodElements
  FlowPastBoxProblem<TTaylorHoodElement<2> > 
   problem(node_file_name, element_file_name, poly_file_name);
#endif
  // Output boundaries 
  problem.mesh_pt()->output_boundaries("RESLT/boundaries.dat");
 
  // Outpt the solution
  problem.doc_solution(doc_info);
  
  // Step number
  doc_info.number()++;
  
#ifdef ADAPT
  // Solve the problem
  problem.newton_solve(max_adapt);
#else
  // Solve the problem
  problem.newton_solve();
#endif
  
  // Outpt the solution
  problem.doc_solution(doc_info);
  
  // Step number
  doc_info.number()++;
  
 } // end of QTaylorHoodElements
 
 
} // end_of_main

References oomph::DocInfo::number(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, problem, oomph::DocInfo::set_directory(), Flag_definition::setup(), and oomph::Global_string_for_annotation::string().