AdaptStationary.cc 1.34 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
#include "AdaptStationary.h"
#include "Parameters.h"
#include "Estimator.h"
#include "TecPlotWriter.h"
#include "ProblemIterationInterface.h"
#include <math.h>

namespace AMDiS {

  AdaptStationary::AdaptStationary(const char *name,
				   ProblemIterationInterface *prob,
				   AdaptInfo  *info) 
    : AdaptBase(name, prob, info)
  {
    initialize();
  }

  int AdaptStationary::adapt()
  {
    FUNCNAME("AdaptStationary::adapt()");

    // initial iteration
23
24
25
26
27
    if (adaptInfo_->getSpaceIteration() == -1) {
      problemIteration_->beginIteration(adaptInfo_);
      problemIteration_->oneIteration(adaptInfo_, NO_ADAPTION);
      problemIteration_->endIteration(adaptInfo_);
      adaptInfo_->incSpaceIteration();
28
29
30
    }

    // adaption loop
31
32
33
    while (!adaptInfo_->spaceToleranceReached() && 
	   (adaptInfo_->getSpaceIteration() < adaptInfo_->getMaxSpaceIteration() || 
	    adaptInfo_->getMaxSpaceIteration() < 0) ) {
34

35
36
37
      problemIteration_->beginIteration(adaptInfo_);
      Flag adapted = problemIteration_->oneIteration(adaptInfo_, FULL_ITERATION);
      problemIteration_->endIteration(adaptInfo_);
38
39
40
41
      
      if (!adapted) 
	break;
      
42
      adaptInfo_->incSpaceIteration();
43
    }
44
45
46
47
48
49
50
51
52
53
54
55

    return 0;
  }

  void AdaptStationary::initialize()
  {
    FUNCNAME("AdaptStationary::initialize()");
  
    GET_PARAMETER(0, name_ + "->info", "%d", &info_);
  }

}