ScalarDiscretizationTypeBase.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2005-2008 Stéphane Del Pino

//  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; either version 2, or (at your option)
//  any later version.

//  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.  

//  $Id: ScalarDiscretizationTypeBase.hpp,v 1.4 2008/01/31 21:38:15 delpinux Exp $

#ifndef SCALAR_DISCRETIZATION_TYPE_BASE_HPP
#define SCALAR_DISCRETIZATION_TYPE_BASE_HPP

#include <ErrorHandler.hpp>
#include <string>

class ScalarDiscretizationTypeBase
{
public:
  enum Type {
    undefined       =-1,        
    functionLike    = 0,        
    // finite element-like discretizations
    lagrangianFEM0  =10,        
    lagrangianFEM1  =11,        
    lagrangianFEM2  =12,        
    // finite element-like discretizations
    DGFEM0          =20,        
    DGFEM1          =21,        
    DGFEM2          =22,        
    // Spectral-like discretizations
    spectralLegendre=50         
  };

private:
  Type __type;          
  virtual void __instanciable() const = 0; 
public:
  static Type getDefault(const Type& type)
  {
    return (type==undefined)?lagrangianFEM1:type;
  }

  const ScalarDiscretizationTypeBase&
  operator=(const ScalarDiscretizationTypeBase::Type& type)
  {
    __type=type;
    return *this;
  }

  static std::string name(const ScalarDiscretizationTypeBase::Type& type)
  {
    switch (type) {
    case lagrangianFEM0:  return "FEM-P0";
    case lagrangianFEM1:  return "FEM-P1";
    case lagrangianFEM2:  return "FEM-P2";
    case DGFEM0:          return "DG-P0";
    case DGFEM1:          return "DG-P1";
    case DGFEM2:          return "DG-P2";
    case spectralLegendre:return "Legendre";
    case undefined:       return "undefined";
    case functionLike:    return "function-like";
    }
    throw ErrorHandler(__FILE__,__LINE__,
                       "unknown Discretization type",
                       ErrorHandler::unexpected);
  }

  static std::string name(const ScalarDiscretizationTypeBase& discretization)
  {
    return ScalarDiscretizationTypeBase::name(discretization.type());
  }

  const Type& type() const
  {
    return __type;
  }

  explicit ScalarDiscretizationTypeBase(const ScalarDiscretizationTypeBase::Type& type)
    : __type(type)
  {
    ;
  }

  explicit ScalarDiscretizationTypeBase(const ScalarDiscretizationTypeBase& d)
    : __type(d.__type)
  {
    ;
  }

  virtual ~ScalarDiscretizationTypeBase()
  {
    ;
  }
};

#endif // SCALAR_DISCRETIZATION_TYPE_BASE_HPP

---