two_d_unsteady_heat_restarted.cc File Reference

#include "generic.h"
#include "unsteady_heat.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	UnsteadyHeatProblem< ELEMENT >
	UnsteadyHeat problem. More...

Namespaces
	ExactSolnForUnsteadyHeat
	Namespace for unforced exact solution for UnsteadyHeat equation.

Functions
void	ExactSolnForUnsteadyHeat::get_exact_u (const double &time, const Vector< double > &x, Vector< double > &u)
	Exact solution as a Vector. More...
void	ExactSolnForUnsteadyHeat::get_exact_u (const double &time, const Vector< double > &x, double &u)
	Exact solution as a scalar. More...
void	ExactSolnForUnsteadyHeat::get_source (const double &time, const Vector< double > &x, double &source)
	Source function to make it an exact solution. More...
int	main (int argc, char *argv[])

Function Documentation

main()

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

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Driver code for unsteady heat equation with option for restart from disk: Only a single command line argument is allowed. If specified it is interpreted as the name of the restart file.

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Build problem
 UnsteadyHeatProblem<QUnsteadyHeatElement<2,4> >
  problem(&ExactSolnForUnsteadyHeat::get_source);
 
 // Setup labels for output
 DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
 
 // Output number
 doc_info.number()=0;
 
 // Open a trace file
 ofstream trace_file;
 char filename[100];   
 sprintf(filename,"%s/trace.dat",doc_info.directory().c_str());
 trace_file.open(filename);
 trace_file << "VARIABLES=\"time\",\"u<SUB>FE</SUB>\","
            << "\"u<SUB>exact</SUB>\",\"norm of error\",\"norm of solution\""
            << std::endl;
 
 // Choose simulation interval
 double t_max=0.5;
 
 // Set IC
 problem.set_initial_condition();
 
 //Output initial condition
 problem.doc_solution(doc_info,trace_file);
 
 //Increment counter for solutions 
 doc_info.number()++;
 
 
 // Get (fixed) timestep dt from Time object. dt was set in 
 // set_initial_conditions(), either by assigning a 
 // value or re-cycling the one from the restart file.
 // This is necessary because the history values obtained from the
 // restart file are only correct if the same dts are used.
 double dt=problem.time_pt()->dt();
 
 // Find number of steps
 unsigned nstep = unsigned(t_max/dt);
 
 // Timestepping loop
 for (unsigned istep=0;istep<nstep;istep++)
  {
   // Take timestep
   problem.unsteady_newton_solve(dt);
   
   //Output solution
   problem.doc_solution(doc_info,trace_file);
   
   //Increment counter for solutions 
   doc_info.number()++;
  }
 
 // Close trace file
 trace_file.close();
 
 
}; // end of main

References oomph::DocInfo::directory(), MergeRestartFiles::filename, ExactSolnForUnsteadyHeat::get_source(), oomph::DocInfo::number(), problem, oomph::DocInfo::set_directory(), and Flag_definition::setup().