two_d_linear_wave.cc File Reference

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

Classes
class	LinearWaveProblem< ELEMENT, TIMESTEPPER >
	LinearWave problem in rectanglular domain. More...

Namespaces
	TanhSolnForLinearWave

Functions
double	TanhSolnForLinearWave::exact_u (const double &time, const Vector< double > &x)
	Exact solution. More...
double	TanhSolnForLinearWave::exact_dudt (const double &time, const Vector< double > &x)
	1st time-deriv of exact solution More...
double	TanhSolnForLinearWave::exact_d2udt2 (const double &time, const Vector< double > &x)
	2nd time-deriv of exact solution More...
void	TanhSolnForLinearWave::get_exact_u (const double &time, const Vector< double > &x, Vector< double > &u)
	Exact solution as a vector. More...
void	TanhSolnForLinearWave::get_source (const double &time, const Vector< double > &x, double &source)
	Source function to make it an exact solution. More...
int	main (int argc, char *argv[])
	Demonstrate how to solve LinearWave problem. More...

Variables
double	TanhSolnForLinearWave::Alpha
	Parameter for steepness of step. More...
double	TanhSolnForLinearWave::Phi
	Orientation of step wave. More...

Function Documentation

main()

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

Demonstrate how to solve LinearWave problem.

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

{
 
 // Store command line arguments: If a command line argument is specied
 // we regard this as validation run.
 CommandLineArgs::setup(argc,argv);
 
 // Number of elements in x direction
 unsigned n_x=10;
 
 // Number of elements in y direction
 unsigned n_y=20;
 
 // Impulsive start?
 bool impulsive_start;
 
 // Run with impulsive start
 // ------------------------
 {
  impulsive_start=true;
 
  // Build problem
  LinearWaveProblem<QLinearWaveElement<2,3>, Newmark<1> >
   problem(n_x,n_y,impulsive_start,&TanhSolnForLinearWave::get_source);
  
  // Run it
  problem.unsteady_run();
 }
 
 // Run with "smooth" start
 // -----------------------
 {
  impulsive_start=false;
 
  // Build problem
  LinearWaveProblem<QLinearWaveElement<2,3>, Newmark<1> >
   problem(n_x,n_y,impulsive_start,&TanhSolnForLinearWave::get_source);
  
  // Run it
  problem.unsteady_run();
 }
 
 
}; // end of main

References TanhSolnForLinearWave::get_source(), problem, and Flag_definition::setup().