one_d_poisson_adapt.cc File Reference

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

Classes
class	RefineableOneDPoissonProblem< ELEMENT >

Namespaces
	ArcTanSolnForPoisson
	Namespace for exact solution for Poisson equation with "sharp step".

Functions
void	ArcTanSolnForPoisson::get_exact_u (const Vector< double > &x, Vector< double > &u)
	Exact solution as a Vector. More...
void	ArcTanSolnForPoisson::get_source (const Vector< double > &x, double &source)
	Source function required to make the solution above an exact solution. More...
int	main ()
	Driver code for 1D Poisson problem. More...

Variables
double	ArcTanSolnForPoisson::Alpha =100.0
	Parameter for steepness of "step". More...

Function Documentation

main()

int main

(

)

Driver code for 1D Poisson problem.

{
 
 //Set up the problem
 //------------------
 
 // Create the problem with 1D three-node refineable elements from the
 // RefineableLinePoissonElement family. Pass pointer to source function. 
 RefineableOneDPoissonProblem<RefineableQPoissonElement<1,3> > 
  problem(&ArcTanSolnForPoisson::get_source);
 
 // Create label for output
 // -----------------------
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 // Step number
 doc_info.number()=0;
 
 // Check if we're ready to go:
 // ---------------------------
 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);
  }
 
 // Choose a large value for the steepness of the "step"
 ArcTanSolnForPoisson::Alpha=100.0; 
 
 // Refine problem uniformly 5 times (to check automatic unrefinement)
 for(unsigned i=0;i<4;i++) { problem.refine_uniformly(); }
 
 // Solve the problem, performing up to 10 adaptive refinements
 problem.newton_solve(10);
 
 // Output the solution
 problem.doc_solution(doc_info);
 
} // End of main

References ArcTanSolnForPoisson::Alpha, ArcTanSolnForPoisson::get_source(), i, oomph::DocInfo::number(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, problem, and oomph::DocInfo::set_directory().