global.h

Go to the documentation of this file.

// This file is part of Hermes2D
//
// Copyright (c) 2009 hp-FEM group at the University of Nevada, Reno (UNR).
// Email: hpfem-group@unr.edu, home page: http://hpfem.org/.
//
// Hermes2D 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 of the License,
// or (at your option) any later version.
//
// Hermes2D 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 Hermes2D; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
#ifndef __H2D_COMMON_H_
#define __H2D_COMMON_H_

#include "hermes_common.h"
#include "config.h"

#define H2D_MAX_ELEMENT_SONS 4 
#define H2D_MAX_NUMBER_EDGES 4 
#define H2D_MAX_NUMBER_VERTICES 4 


#define H2D_NUM_MODES 2 
#define H2D_SOLUTION_ELEMENT_CACHE_SIZE 2 
#define H2D_MAX_NODE_ID 10000000
#define H2D_MAX_SOLUTION_COMPONENTS 2

#define HERMES_ONE NULL
#define HERMES_DEFAULT_FUNCTION NULL
#define HERMES_DEFAULT_SPLINE NULL

#define HERMES_TOTAL_ERROR_REL  0x00  


#define HERMES_TOTAL_ERROR_ABS  0x01  


#define HERMES_ELEMENT_ERROR_REL 0x00 




#define HERMES_ELEMENT_ERROR_ABS 0x10 




    static const int H2D_TRF_TRI_NUM = 4;
    static const int H2D_TRF_QUAD_NUM = 8;
    static const int H2D_TRF_NUM = (H2D_TRF_QUAD_NUM + 1);
    static const int H2D_TRF_IDENTITY = H2D_TRF_QUAD_NUM;

#define H2DRS_ASSUMED_MAX_CANDS 512 

//TODO: find out why 20 used used, should'n be there 2*(H2DRS_MAX_ORDER+1)
#define H2DRS_INTR_GIP_ORDER 20 
#define H2DRS_MAX_ORDER_INC 1 

#define H2DRS_SCORE_DIFF_ZERO 1E-13 

#define H2DRS_ORDER_ANY -1 

# define H2DRS_DEFAULT_ERR_WEIGHT_H 2.0 
# define H2DRS_DEFAULT_ERR_WEIGHT_P 1.0 
# define H2DRS_DEFAULT_ERR_WEIGHT_ANISO 1.414214 

namespace Hermes
{
  namespace Hermes2D
  {
    class RefMap;
    template<typename Scalar> class DiscreteProblem;
    template<typename Scalar> class Space;
    template<typename Scalar> class WeakForm;
    template<typename Scalar> class MeshFunction;
    template<typename Scalar> class Solution;
    class Quad2D;
    class Quad1DStd;
    class Quad2DStd;

    const int H2D_ORDER_BITS = 5;
    const int H2D_ORDER_MASK = (1 << H2D_ORDER_BITS) - 1;

    #define H2D_GET_H_ORDER(encoded_order) ((encoded_order) & H2D_ORDER_MASK)
    #define H2D_GET_V_ORDER(encoded_order) ((encoded_order) >> H2D_ORDER_BITS)
    #define H2D_MAKE_QUAD_ORDER(h_encoded_order, v_encoded_order) (((v_encoded_order) << H2D_ORDER_BITS) + (h_encoded_order))
    #define H2D_MAKE_EDGE_ORDER(mode, edge, order) ((mode == HERMES_MODE_TRIANGLE || edge == 0 || edge == 2) ? H2D_GET_H_ORDER(order) : H2D_GET_V_ORDER(order))

    template<typename Scalar>
    class HERMES_API Global : public Hermes::Mixins::Loggable
    {
    public:
      Global() : Hermes::Mixins::Loggable() {};
      friend void warn_order();
    public:
      // Note: coarse mesh sln has to be first, then
      // ref_sln (because the abs. error is divided
      // by the norm of the latter).
      static double calc_rel_error(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, int norm_type);

      static double calc_abs_error(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, int norm_type);

      static double calc_norm(MeshFunction<Scalar>* sln, int norm_type);

      static double calc_norms(Hermes::vector<Solution<Scalar>*> slns);
      static double calc_abs_errors(Hermes::vector<Solution<Scalar>*> slns1, Hermes::vector<Solution<Scalar>*> slns2);
      static double calc_rel_errors(Hermes::vector<Solution<Scalar>*> slns1, Hermes::vector<Solution<Scalar>*> slns2);

      static double error_fn_l2(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, RefMap* ru, RefMap* rv);
      static double norm_fn_l2(MeshFunction<Scalar>* sln, RefMap* ru);

      static double error_fn_h1(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, RefMap* ru, RefMap* rv);
      static double norm_fn_h1(MeshFunction<Scalar>* sln, RefMap* ru);

      static double error_fn_hc(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, RefMap* ru, RefMap* rv);
      static double norm_fn_hc(MeshFunction<Scalar>* sln, RefMap* ru);

      static double error_fn_hcl2(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, RefMap* ru, RefMap* rv);
      static double norm_fn_hcl2(MeshFunction<Scalar>* sln, RefMap* ru);

      static double error_fn_hdiv(MeshFunction<Scalar>* sln1, MeshFunction<Scalar>* sln2, RefMap* ru, RefMap* rv);
      static double norm_fn_hdiv(MeshFunction<Scalar>* sln, RefMap* ru);

      static double get_l2_norm(Vector<Scalar>* vec);
    };

    enum ProjNormType
    {
      HERMES_L2_NORM,
      HERMES_H1_NORM,
      HERMES_H1_SEMINORM,
      HERMES_HCURL_NORM,
      HERMES_HDIV_NORM,
      HERMES_UNSET_NORM
    };

    enum ElementMode2D {
      HERMES_MODE_TRIANGLE = 0,
      HERMES_MODE_QUAD = 1
    };
  }
}
#endif