mpi/distribution/eigenproblem/harmonic.cc File Reference

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

Classes
class	ComplexLess< T >
class	HarmonicEquations
class	QHarmonicElement< NNODE_1D >
class	HarmonicProblem< ELEMENT, EIGEN_SOLVER >
	1D Harmonic problem in unit interval. More...

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

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver for 1D Poisson problem.

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

{
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 // Set up the problem: 
 unsigned n_element=100; //Number of elements
 
#ifdef OOMPH_HAS_TRILINOS
//Solve with Anasazi (non distributed)
 clock_t t_start1 = clock();
 {
  HarmonicProblem<QHarmonicElement<3>,ANASAZI> problem(n_element);
  
  problem.solve(1);
 }
 clock_t t_end1 = clock();
 
 
 //Solve with Anasazi (distributed)
 clock_t t_start2 = clock();
 {
 
  HarmonicProblem<QHarmonicElement<3>,ANASAZI> problem(n_element);
  //In the validation case there should be two processors
  if(problem.communicator_pt()->nproc()==2)
   {
    //Setup a dummy partition
    unsigned n_element = problem.mesh_pt()->nelement();
    Vector<unsigned> element_partition(n_element);
    for(unsigned e=0;e<n_element/2;e++) {element_partition[e]=0;}
    for(unsigned e=n_element/2;e<n_element;e++) {element_partition[e]=1;}
    problem.distribute(element_partition);
   }
  else
   {
    problem.distribute();
   }
  problem.solve(2);
 }
 clock_t t_end2 = clock();
 
 oomph_info << "ANASAZI TIME (non-distributed): " 
            << (double)(t_end1 - t_start1)/CLOCKS_PER_SEC
            << std::endl;
 
 
 oomph_info << "ANASAZI TIME (distributed): " << 
  (double)(t_end2 - t_start2)/CLOCKS_PER_SEC
            << std::endl;
#endif
 
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::finalize();
#endif
 
} // end of main

References e(), oomph::oomph_info, and problem.