pml_helmholtz/scattering/unstructured_two_d_helmholtz.cc File Reference

#include <fenv.h>
#include "math.h"
#include <complex>
#include "generic.h"
#include "pml_helmholtz.h"
#include "meshes/triangle_mesh.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	PMLProblem< ELEMENT >

Namespaces
	GlobalParameters
	Global parameters.

Functions
int	main (int argc, char **argv)
	Solve 2D Helmholtz problem. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Solve 2D Helmholtz problem.

{
 //Set up the problem
 //------------------
 
#ifdef ADAPTIVE
 
 // Set up the problem with projectable 2D six-node elements from the
 // TPMLHelmholtzElement family.
 PMLProblem<ProjectablePMLHelmholtzElement
  <TPMLHelmholtzElement<2,3> > > problem;
 
 // Set up the problem with 2D ten-node elements from the
 // TPMLHelmholtzElement family.
 // PMLProblem<ProjectablePMLHelmholtzElement
 //  <TPMLHelmholtzElement<2,4> > > problem;
 
#else
 
 // Set up the problem with 2D six-node elements from the
 // TPMLHelmholtzElement family.
 PMLProblem<TPMLHelmholtzElement<2,3> >  problem;
 
 // Set up the problem with 2D ten-node elements from the
 // TPMLHelmholtzElement family.
 //   PMLProblem<TPMLHelmholtzElement<2,4> >  problem;
 
#endif
 
 // Create label for output
 //------------------------
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 
#ifdef ADAPTIVE
 
 // Max. number of adaptations
 unsigned max_adapt=1;
 
 // Solve the problem with the adaptive Newton's method, allowing
 // up to max_adapt mesh adaptations after every solve.
 problem.newton_solve(max_adapt);
 
#else
 
 // Solve the problem with Newton's method
 problem.newton_solve();
 
#endif
 
 //Output solution
 problem.doc_solution(doc_info);
 
} //end of main

References problem, and oomph::DocInfo::set_directory().