FixedPointIteration.h

Go to the documentation of this file.

// Copyright (C) 2003 Johan Hoffman and Anders Logg.
// Licensed under the GNU GPL Version 2.

#ifndef __FIXED_POINT_ITERATION_H
#define __FIXED_POINT_ITERATION_H

#include <dolfin/constants.h>
#include <dolfin/Iteration.h>

namespace dolfin
{
  class Solution;
  class RHS;
  class Adaptivity;
  class TimeSlab;
  class Element;
  class ElementGroup;
  class ElementGroupList;
  

  class FixedPointIteration
  {
  public:

    FixedPointIteration(Solution& u, RHS& f, Adaptivity& adaptivity);
    
    ~FixedPointIteration();

    bool iterate(TimeSlab& timeslab);

    bool iterate(ElementGroupList& list);

    real iterate(ElementGroup& group);
    
    real iterate(Element& element);
    
    void reset(ElementGroupList& list);

    void reset(ElementGroup& group);

    void reset(Element& element);

    void stabilization(real& alpha, unsigned int& m) const;

    void report() const;

  private:

    // Start iteration on group list
    void start(ElementGroupList& list);

    // Start iteration on element group
    void start(ElementGroup& group);

    // Start iteration on element
    void start(Element& element);

    // End iteration on group list
    void end(ElementGroupList& list);

    // End iteration on element group
    void end(ElementGroup& group);

    // End iteration on element
    void end(Element& element);

    // Update group list
    void update(ElementGroupList& list, Iteration::Increments& d);
    
    // Update element group
    void update(ElementGroup& group, Iteration::Increments& d);

    // Update element
    void update(Element& element, Iteration::Increments& d);

    // Stabilize group list
    void stabilize(ElementGroupList& list, const Iteration::Increments& d, 
                   unsigned int n);
    
    // Stabilize element group
    void stabilize(ElementGroup& group, const Iteration::Increments& r,
                   unsigned int n);

    // Stabilize element
    void stabilize(Element& element, const Iteration::Increments& d,
                   unsigned int n);

    // Check convergence for group list
    bool converged(ElementGroupList& list, const Iteration::Increments& d,
                   unsigned int n);
    
    // Check convergence for element group
    bool converged(ElementGroup& group, const Iteration::Increments& d,
                   unsigned int n);

    // Check convergence for element
    bool converged(Element& element, const Iteration::Increments& d,
                   unsigned int n);

    // Check divergence for group list
    bool diverged(ElementGroupList& list, const Iteration::Increments& d,
                  unsigned int n, Iteration::State& newstate);
    
    // Check divergence for element group
    bool diverged(ElementGroup& group, const Iteration::Increments& d,
                  unsigned int n, Iteration::State& newstate);
    
    // Check divergence for element
    bool diverged(Element& element, const Iteration::Increments& d,
                  unsigned int n, Iteration::State& newstate);
    
    void init(ElementGroup& group);
    
    // Update initial data for element
    void init(Element& element);

    // Change state
    void changeState(Iteration::State newstate);

    // Write debug info for group list iteration
    void debug(ElementGroupList& list, const Iteration::Increments& d,
               unsigned int n);
    
    // Write debug info for element group iteration
    void debug(ElementGroup& group, const Iteration::Increments& d,
               unsigned int n);
    
    // Write debug info for element iteration
    void debug(Element& element, const Iteration::Increments& d,
               unsigned int n);
    
    //--- Data for fixed point iteration ---

    // Solution
    Solution& u;

    // Right-hand side f
    RHS& f;

    // Maximum number of iterations
    unsigned int maxiter;

    // Maximum allowed divergence
    real maxdiv;

    // Maximum allowed convergence
    real maxconv;

    // Tolerance for discrete residual
    real tol;

    // Check if we should debug iterations
    bool debug_iter;

    // Current state
    Iteration* state;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.