precision.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_PRECISION_HH
#define DUNE_PRECISION_HH
 
#include <stdlib.h>
 
namespace Dune {
 
  template <class ctype = double>
  class FMatrixPrecision {
  public:
    static ctype pivoting_limit ()
    {
      return _pivoting;
    }
 
    static void set_pivoting_limit (ctype pivthres)
    {
      _pivoting = pivthres;
    }
 
    static ctype singular_limit ()
    {
      return _singular;
    }
 
    static void set_singular_limit (ctype singthres)
    {
      _singular = singthres;
    }
 
    static ctype absolute_limit ()
    {
      return _absolute;
    }
 
    static void set_absolute_limit (ctype absthres)
    {
      _absolute = absthres;
    }
 
  private:
    // just to demonstrate some state information
    static ctype _pivoting;
    static ctype _singular;
    static ctype _absolute;
  };
 
  template <class ctype>
  ctype FMatrixPrecision<ctype>::_pivoting = 1E-8;
  template <class ctype>
  ctype FMatrixPrecision<ctype>::_singular = 1E-14;
  template <class ctype>
  ctype FMatrixPrecision<ctype>::_absolute = 1E-80;
 
} // end namespace
 
#endif