complex_harmonic.cc File Reference

#include "generic.h"
#include "meshes/one_d_mesh.h"

Classes
class	ComplexLess< T >
class	ComplexHarmonicEquations
class	QComplexHarmonicElement< NNODE_1D >
class	ComplexHarmonicProblem< ELEMENT, EIGEN_SOLVER >
	1D ComplexHarmonic problem in unit interval. More...

Namespaces
	EigenproblemShift
	Namespace for the shift applied to the eigenproblem.

Functions
int	main (int argc, char **argv)
	Driver for 1D Poisson problem. More...

Variables
double	EigenproblemShift::Mu = 6.5

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver for 1D Poisson problem.

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

{
//Want to test Trilinos if we have it, so we must initialise MPI
//if we have compiled with it
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 
 // Set up the problem: 
 unsigned n_element=100; //Number of elements
 
 clock_t t_start1 = clock();
 //Solve with ARPACK
 {
  ComplexHarmonicProblem<QComplexHarmonicElement<3>,ARPACK> 
   problem(n_element);
  
  std::cout << "Matrix size " << problem.ndof() << std::endl;
  
  problem.solve(1);
  }
 clock_t t_end1 = clock();
 
 clock_t t_start2 = clock();
 //Solve with LAPACK_QZ
 {
  ComplexHarmonicProblem<QComplexHarmonicElement<3>,LAPACK_QZ> 
   problem(n_element);
  
  problem.solve(2);
 }
 clock_t t_end2 = clock();
 
#ifdef OOMPH_HAS_TRILINOS
 clock_t t_start3 = clock();
//Solve with Anasazi
 {
  ComplexHarmonicProblem<QComplexHarmonicElement<3>,ANASAZI> problem(n_element);
  problem.solve(3);
 }
 clock_t t_end3 = clock();
#endif
 
 std::cout << "ARPACK TIME: " << (double)(t_end1 - t_start1)/CLOCKS_PER_SEC
           << std::endl;
 
 std::cout << "LAPACK TIME: " << (double)(t_end2 - t_start2)/CLOCKS_PER_SEC
           << std::endl;
 
#ifdef OOMPH_HAS_TRILINOS
  std::cout << "ANASAZI TIME: " << (double)(t_end3 - t_start3)/CLOCKS_PER_SEC
           << std::endl;
#endif
 
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::finalize();
#endif
 
} // end of main

References problem.