root/old_code/tumor_growth/heat.cc

DEFINITIONS

1. main

   1 #include <config.h>
   2 
   3 // iostream includes
   4 #include <iostream>
   5 
   6 // include GeometryGrid
   7 #include <dune/grid/geometrygrid.hh>
   8 
   9 // include grid part
  10 #include <dune/fem/gridpart/adaptiveleafgridpart.hh>
  11 
  12 // include discrete function space
  13 //#include <dune/fem/space/lagrangespace.hh> veraltet
  14 #include <dune/fem/space/lagrange.hh>
  15 
  16 // include discrete function
  17 #include <dune/fem/function/adaptivefunction.hh>
  18 
  19 // include solvers
  20 //#include <dune/fem/solver/inverseoperators.hh> veraltet
  21 #include <dune/fem/solver/cginverseoperator.hh>
  22 #include <dune/fem/solver/oemsolver.hh>
  23 
  24 // include Lagrange interpolation
  25 #include <dune/fem/operator/lagrangeinterpolation.hh>
  26 
  27 // include norms
  28 #include <dune/fem/misc/l2norm.hh>
  29 #include <dune/fem/misc/h1norm.hh>
  30 
  31 // include parameter handling
  32 #include <dune/fem/io/parameter.hh>
  33 
  34 // include output
  35 #include <dune/fem/io/file/dataoutput.hh>
  36 
  37 // local includes
  38 #include "deformation.hh"
  39 #include "elliptic.hh"
  40 #include "heat.hh"
  41 #include "rhs.hh"
  42 // local CAPG includes
  43 #include "dune_typedef_heat.hpp"
  44 // this header only contains lots of typedef's
  45 #include "dune_bdf.hpp"
  46 #include "dune_heat_algorithm.hpp"
  47 
  48 // main
  49 // ----
  50 
  51 int main ( int argc, char **argv )
  52 try
  53 {
  54   // initialize MPI, if necessary
  55   Dune::Fem::MPIManager::initialize(argc, argv);
  56 
  57   // append all given parameters
  58   Dune::Fem::Parameter::append(argc, argv);
  59   for( int i = 1; i < argc; ++i )
  60     Dune::Fem::Parameter::append(argv[i]);
  61   Dune::Fem::Parameter::append("../data/parameter");
  62 
  63   // BDF::ie_heat_algorithm();
  64   // BDF::bdf2_heat_algorithm();  
  65   // BDF::nonlinear_algorithm();
  66   // nonlinear_evolution();
  67   tumor_growth_code();
  68 
  69   return 0;
  70 }
  71 catch( const Dune::Exception& e )
  72 {
  73   std::cerr << "Error: " << e << std::endl;
  74   return 1;
  75 }