fish_poisson_node_update.cc File Reference

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

Classes
class	RefineableFishPoissonProblem< ELEMENT >

Namespaces
	ConstSourceForPoisson
	Namespace for const source term in Poisson equation.

Functions
void	ConstSourceForPoisson::get_source (const Vector< double > &x, double &source)
	Const source function. More...
int	main ()

Function Documentation

main()

int main

(

)

Demonstrate how to solve 2D Poisson problem in fish-shaped domain with black-box mesh adaptation and domain updates in response to changes in the domain shape.

{
 
 //Set up the problem with 9 node refineable Poisson elements
 RefineableFishPoissonProblem<RefineableQPoissonElement<2,3> > problem;
 
 // Setup labels for output
 //------------------------
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT"); 
 
 // Adjust the domain shape by changing the width of the fish
 //----------------------------------------------------------
 unsigned nstep=3;
 for (unsigned i=0;i<nstep;i++)
  {
   // Get pointer to GeomObject that defines the position of the
   // fish's back:
   GeomObject* fish_back_pt=problem.mesh_pt()->fish_back_pt();
 
   // Recast to pointer to Circle object to get access to the member function
   // that sets the y-position of the Circle's centre and decrease its
   // value, making the fish narrower
   dynamic_cast<Circle*>(fish_back_pt)->y_c()-=0.1;
   
   // Update the domain shape in response to the changes in its
   // boundary
   problem.mesh_pt()->node_update();
 
   // Solve the problem, allowing for up to two levels of refinement
   problem.newton_solve(2);
   
   //Output solution
   problem.doc_solution(doc_info);
   
   //Increment counter for solutions 
   doc_info.number()++;
  }
 
} // end of main

References i, oomph::DocInfo::number(), problem, and oomph::DocInfo::set_directory().