limit_order.cpp

#include "global.h"
#include "quad_all.h"
#include "limit_order.h"

namespace Hermes
{
  namespace Hermes2D
  {
    static unsigned short default_order_table_tri[] =
    {
      0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,
      17, 18, 19, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20,
      20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20,
      20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20,
      20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20
    };

#ifdef EXTREME_QUAD
    static unsigned short default_order_table_quad[] =
    {
      1, 1, 3, 3, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15,
      17, 17, 19, 19, 21, 21, 23, 23, 25, 25, 27, 27, 29, 29, 31, 31,
      33, 33, 35, 35, 37, 37, 39, 39, 41, 41, 43, 43, 45, 45, 47, 47,
      49, 49, 51, 51, 53, 53, 55, 55, 57, 57, 59, 59, 61, 61, 63, 63,
      65, 65, 67, 67, 69, 69, 71, 71, 73, 73, 75, 75, 77, 77, 79, 79,
      81, 81, 83, 83, 85, 85, 87, 87, 89, 89, 91, 91, 93, 93, 95, 95,
      97, 97, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99
    };
#else
    static unsigned short default_order_table_quad[] =
    {
      1, 1, 3, 3, 5, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15, 17,
      17, 19, 19, 21, 21, 23, 23, 24, 24, 24, 24, 24, 24, 24,
      24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24,
      24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24,
      24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
    };
#endif

    static unsigned short* g_order_table_quad = default_order_table_quad;
    static unsigned short* g_order_table_tri = default_order_table_tri;
    static bool warned_order = false;

    HERMES_API unsigned short  g_max_order;
    HERMES_API unsigned short  g_safe_max_order;

    HERMES_API void set_order_limit_table(unsigned short* tri_table, unsigned short* quad_table, unsigned short n)
    {
      if (n < 24)
        throw Hermes::Exceptions::Exception("Hermes::Order limit tables must have at least 24 entries.");
      g_order_table_tri = tri_table;
      g_order_table_quad = quad_table;
    }

    HERMES_API void update_limit_table(ElementMode2D mode)
    {
      g_max_order = g_quad_2d_std.get_max_order(mode);
      g_safe_max_order = g_quad_2d_std.get_safe_max_order(mode);
    }

    HERMES_API void reset_warn_order()
    {
      warned_order = false;
    }

    HERMES_API void warn_order()
    {
      if (HermesCommonApi.get_integral_param_value(Hermes::showInternalWarnings))
        if (!warned_order)
        {
#pragma omp critical (warn_oder)
        if (!warned_order)
        {
          Hermes::Mixins::Loggable::Static::warn("Warning: Not enough integration rules for exact integration.");
          warned_order = true;
        }
        }
    }

    HERMES_API void limit_order(int& o, ElementMode2D mode)
    {
      if (o >= g_quad_2d_std.get_safe_max_order(mode))
      {
        o = g_quad_2d_std.get_safe_max_order(mode);
        warn_order();
      }
      if (mode == HERMES_MODE_TRIANGLE)
        o = g_order_table_tri[o];
      else
        o = g_order_table_quad[o];
    }

    HERMES_API void limit_order_nowarn(int& o, ElementMode2D mode)
    {
      if (o > g_quad_2d_std.get_safe_max_order(mode))
        o = g_quad_2d_std.get_safe_max_order(mode);
      if (mode == HERMES_MODE_TRIANGLE)
        o = g_order_table_tri[o];
      else
        o = g_order_table_quad[o];
    }
  }
}