An Object Finite Element Library

Estimator.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free
   Software Foundation; Version 2 of the License.

   This program is distributed in the hope that it will be useful, but WITHOUT
   ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
   FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
   details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

              Definition of class 'Estimator' for error post-estimation

  ==============================================================================*/


#ifndef __ESTIMATOR_H
#define __ESTIMATOR_H

#include "OFELI_Config.h"
#include "Mesh.h"
#include "NodeVect.h"
#include "ElementVect.h"
#include "SideVect.h"


namespace OFELI {

#ifndef DOXYGEN_SHOULD_SKIP_THIS
struct ErrorInEstimator {
   void Message(const char *file, size_t line, int code, int p1=0);
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */


//template class NodeVect<double>;
//template class ElementVect<double>;
//template class NodeVect<Point<double> >;

class Estimator
{

public :

    Estimator() {
#ifdef _OFELI_DEBUG
      std::clog << "An instance of class Estimator is constructed.\n";
      std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    Estimator(Mesh &m);

    ~Estimator() {
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Estimator is destructed.\n";
       std::clog << "File: " << __FILE__ << ", Line: " << __LINE__ << endl;
#endif
    }

    void setError(const NodeVect<double> &u);

    ElementVect<double> Err;

    double getAverage() const { return average_; }

    const Mesh &getMesh() const { return *theMesh_; }

private :

   Mesh *theMesh_;
   double average_;
   ErrorInEstimator e_;
   void elementT3_(const NodeVect<double> &u, Vect<double> &M, NodeVect<Point<double> > &b);
};


ostream& operator<<(ostream& s, const Estimator &r);

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani