helmholtz_point_source.cc File Reference

#include <fenv.h>
#include "math.h"
#include <complex>
#include "generic.h"
#include "helmholtz.h"
#include "meshes/triangle_mesh.h"
#include "oomph_crbond_bessel.h"

Classes
class	oomph::HelmholtzPointSourceElement< ELEMENT >
	Class to impose point source to (wrapped) Helmholtz element. More...
class	oomph::FaceGeometry< HelmholtzPointSourceElement< ELEMENT > >
class	oomph::FaceGeometry< FaceGeometry< HelmholtzPointSourceElement< ELEMENT > > >
class	HelmholtzPointSourceProblem< ELEMENT >
	Problem class. More...

Namespaces
	GlobalParameters
	Global parameters.
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).

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

Variables
double	GlobalParameters::K_squared =10.0
	Square of the wavenumber. More...
bool	GlobalParameters::DtN_BC =false
unsigned	GlobalParameters::ABC_order =3
	Flag to choose wich order to use. More...
double	GlobalParameters::Outer_radius =1.5
	Radius of outer boundary (must be a circle!) More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Solve 2D Helmholtz problem for scattering of a planar wave from a unit disk

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Define case to be run
 unsigned i_case=0;
 CommandLineArgs::specify_command_line_flag("--case",&i_case);
 
 // Parse command line
 CommandLineArgs::parse_and_assign(); 
 
 // Doc what has actually been specified on the command line
 CommandLineArgs::doc_specified_flags();
 
 // Now set flags accordingly
 switch(i_case)
  {
  case 0:
   GlobalParameters::DtN_BC=true;
   break;
   
  case 1:
   GlobalParameters::DtN_BC=false;
   GlobalParameters::ABC_order=1;
   break;
 
  case 2:
   GlobalParameters::DtN_BC=false;
   GlobalParameters::ABC_order=2;
   break;
 
  case 3:
   GlobalParameters::DtN_BC=false;
   GlobalParameters::ABC_order=3;
   break;
  }
 
 
 //Set up the problem
 //------------------
 
#ifdef ADAPTIVE
 
 //Set up the problem with 2D nine-node elements from the
 //QHelmholtzElement family.
 HelmholtzPointSourceProblem<ProjectableHelmholtzElement
                             <HelmholtzPointSourceElement
                              <THelmholtzElement<2,3> > > >
  problem;
 
#else
 
 // Set up the problem with 2D six-node elements from the
 // THelmholtzElement family. 
 HelmholtzPointSourceProblem<HelmholtzPointSourceElement
                             <THelmholtzElement<2,3> > > 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=5;
 
   // Solve the problem with 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 GlobalParameters::ABC_order, oomph::CommandLineArgs::doc_specified_flags(), GlobalParameters::DtN_BC, oomph::CommandLineArgs::parse_and_assign(), problem, oomph::DocInfo::set_directory(), oomph::CommandLineArgs::setup(), and oomph::CommandLineArgs::specify_command_line_flag().