one_d_spectral.cc File Reference

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

Classes
class	OneDPoissonProblem< ELEMENT >
	1D Poisson problem in unit interval. More...

Namespaces
	FishSolnOneDPoisson
	Namespace for fish-shaped solution of 1D Poisson equation.

Functions
void	FishSolnOneDPoisson::get_exact_u (const Vector< double > &x, Vector< double > &u)
	Exact, fish-shaped solution as a 1D vector. More...
void	FishSolnOneDPoisson::source_function (const Vector< double > &x, double &source)
	Source function required to make the fish shape an exact solution. More...
template<unsigned SPECTRAL_ORDER>
void	construct_and_solve_problem (const unsigned &i)
int	main ()
	Driver for 1D Poisson problem. More...

Function Documentation

construct_and_solve_problem()

template<unsigned SPECTRAL_ORDER>

void construct_and_solve_problem

(

const unsigned &

i

)

{
 // Set up the problem: 
 // Solve a 1D Poisson problem using a source function that generates
 // a fish shaped exact solution
 unsigned n_element=10; //Number of elements
 
 OneDPoissonProblem<QSpectralPoissonElement<1,SPECTRAL_ORDER> > //Element type as template parameter
  problem(n_element,FishSolnOneDPoisson::source_function);
 
 // Check whether the problem can be solved
 cout << "\n\n\nProblem self-test ";
 if (problem.self_test()==0)  
  {
   cout << "passed: Problem can be solved." << std::endl;
  }
 else 
  {
   throw OomphLibError("Self test failed",
                       OOMPH_CURRENT_FUNCTION,
                       OOMPH_EXCEPTION_LOCATION);
  }
 
 // Set the sign of the source function:
 cout << "\n\n\nSolving with Spectral Order :\n" << SPECTRAL_ORDER << std::endl;
 FishSolnOneDPoisson::Sign=-1;
 
 // Solve the problem with this Sign
 problem.newton_solve();
 
 //Output solution for this case (label output files with "i")
 problem.doc_solution(i);
}

References i, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, problem, FishSolnOneDPoisson::Sign, and FishSolnOneDPoisson::source_function().

main()

int main

(

)

Driver for 1D Poisson problem.

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

{
 construct_and_solve_problem<2>(0);
 construct_and_solve_problem<3>(1);
 construct_and_solve_problem<4>(2);
 construct_and_solve_problem<5>(3);
 construct_and_solve_problem<6>(4);
 construct_and_solve_problem<7>(5);
 construct_and_solve_problem<8>(6);
} // end of main