old_for_doc.cc File Reference

#include "generic.h"
#include "poisson.h"
#include "meshes/fish_mesh.h"
#include "circle_as_generalised_element.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...
template<class ELEMENT >
void	demo_fish_poisson (const string &directory_name)
template<class ELEMENT >
void	demo_elastic_fish_poisson (const string &directory_name)
int	main (int argc, char *argv[])

Function Documentation

demo_elastic_fish_poisson()

template<class ELEMENT >

void demo_elastic_fish_poisson

(

const string &

directory_name

)

Demonstrate how to solve coupled "elastic" 2D Poisson problem in deformable fish-shaped domain with mesh adaptation.

{
 
  //Set up the problem with "elastic" fish back
 bool fix_position=false;
 RefineableFishPoissonProblem<ELEMENT> problem(fix_position,directory_name);
  
 // Change/doc targets for mesh adaptation
 if (CommandLineArgs::Argc>1) 
  {
   problem.fish_mesh_pt()->max_permitted_error()=0.05;
   problem.fish_mesh_pt()->min_permitted_error()=0.005;
  }
 problem.fish_mesh_pt()->doc_adaptivity_targets(cout);
  
  
 // Do some uniform mesh refinement first
 problem.refine_uniformly();
 problem.refine_uniformly();
 
 // Initial guess for "load" on fish back
 problem.load()=0.0;
 
 // Solve/doc fully coupled problem
 unsigned max_solve=2; 
 problem.newton_solve(max_solve);
 problem.doc_solution();
 
}

References oomph::CommandLineArgs::Argc, and problem.

demo_fish_poisson()

template<class ELEMENT >

void demo_fish_poisson

(

const string &

directory_name

)

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Demonstrate how to solve 2D Poisson problem in deformable fish-shaped domain with mesh adaptation.

{
 
 // Set up the problem with prescribed displacement of fish back
 bool fix_position=true;
 RefineableFishPoissonProblem<ELEMENT> problem(fix_position,directory_name);  
 
 
 // Change/doc targets for mesh adaptation
 if (CommandLineArgs::Argc>1) 
  {
   problem.fish_mesh_pt()->max_permitted_error()=0.05;
   problem.fish_mesh_pt()->min_permitted_error()=0.005;
  }
 else
  {
   problem.fish_mesh_pt()->max_permitted_error()=0.005;
   problem.fish_mesh_pt()->min_permitted_error()=0.0005;
  }
 problem.fish_mesh_pt()->doc_adaptivity_targets(cout);
    
 
 // Do some uniform mesh refinement first
 //--------------------------------------
 problem.refine_uniformly();
 problem.refine_uniformly(); 
 
 // Initial value for the vertical displacement of the fish's back
 problem.y_c()=-0.3;
 
 // Loop for different fish shapes
 //-------------------------------
 
 // Number of steps
 unsigned nstep=5;
 if (CommandLineArgs::Argc>1) nstep=1;
 
 // Increment in displacement
 double dyc=0.6/double(nstep-1);
 
 // Loop over different fish shapes
 for (unsigned istep=0;istep<nstep;istep++)
  {    
   // Solve/doc
   unsigned max_solve=2; 
   problem.newton_solve(max_solve);
   problem.doc_solution();
   
   //Increment counter for solutions 
   problem.doc_info().number()++;    
   
   // Change vertical displacement
   problem.y_c()+=dyc;
  } 
 
}

References oomph::CommandLineArgs::Argc, and problem.

main()

int main	(	int argc	,
		char *	argv[]
	)

Driver for "elastic" fish poisson solver with adaptation. If there are any command line arguments, we regard this as a validation run and reduce the targets for the mesh adaptation.

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Shorthand for element type
 typedef MacroElementNodeUpdateElement<RefineableQPoissonElement<2,3> > ELEMENT;
 
 // Compute solution of Poisson equation in various domains -- prescribed
 // boundary shape.
 demo_fish_poisson<ELEMENT>("RESLT");
 
 // Compute coupled "elastic" coupled solution directly
 demo_elastic_fish_poisson<ELEMENT>("RESLT_coupled"); 
 
}

References Flag_definition::setup().