boxunions.h

Go to the documentation of this file.

/* 
 * SPDX-FileCopyrightText: 2025 Gesellschaft fuer Anlagen- und Reaktorsicherheit gGmbH
 * SPDX-License-Identifier: EUPL-1.2
 * SPDX-FileContributor: Dmitry Logashenko
 * SPDX-FileContributor: Goethe Universität Frankfurt
 * SPDX-FileType: SOURCE
 * 
 * This file is part of d3f++.
 * d3f++ is an extension for UG4. Licensing information and citation requirements of UG4 are provided in LICENSES/UG4-LGPL_2.1
 */
 
 
 
/*
 * Definition of surface functions by unions of boxes
 */
#ifndef __H__UG__PLUGINS__D3F__BOXUNIONS__
#define __H__UG__PLUGINS__D3F__BOXUNIONS__
 
#include <vector>
 
// ug4 headers
#include "common/common.h"
#include "common/math/ugmath.h"
 
namespace ug {
namespace d3f {
 
class UnionsOfBoxes
{
//  Description of the union of boxes
  class t_box_union_item
  {
  public:
    t_box_union_item * tl; 
    int ident; 
    size_t dir; 
    number dx_half; 
    number dy_half; 
    std::vector<MathVector<2> > beg_corner; 
    std::vector<MathVector<2> > end_corner; 
    
    t_box_union_item
    (
      int the_id,
      size_t the_dir,
      number the_dx, number the_dy,
      number x, number y,
      t_box_union_item * the_tl = NULL)
    : tl (the_tl), ident (the_id), dir (the_dir), dx_half (the_dx/2), dy_half (the_dy/2),
      beg_corner (1), end_corner (1)
    {
      beg_corner[0][0] = x; beg_corner[0][1] = y;
      end_corner[0] = beg_corner[0];
    }
    
    void append (number x, number y);
    
    bool pnt_inside (number x, number y) const;
  };
 
public:
 
  UnionsOfBoxes ()
  : m_dir (2), m_dx (1), m_dy (1), m_box_unions (NULL)
  {};
  
  ~UnionsOfBoxes ()
  {
    t_box_union_item * t = m_box_unions;
    while (t != NULL)
    {
      t_box_union_item * s = t->tl;
      delete t;
      t = s;
    }
  }
  
  void set_grid_sizes
  (
    size_t dir, 
    number dx, 
    number dy 
  )
  {
    m_dir = dir;
    m_dx = dx; m_dy = dy;
  }
  
  void append_point (number x, number y, int ident)
  {
    t_box_union_item * t = union_by_id (ident);
    if (t != NULL)
      t->append (x, y);
    else
      m_box_unions = new t_box_union_item (ident, m_dir, m_dx, m_dy, x, y, m_box_unions);
  }
 
  bool get_id_by_pnt (number x, number y, int & ident) const
  {
    for (t_box_union_item * t = m_box_unions; t != NULL; t = t->tl)
      if (t->pnt_inside (x, y))
      {
        ident = t->ident;
        return true;
      }
    return false;
  }
  
  int get_id_by_pnt_dv (number x, number y, int dv) const
  {
    int ident;
    return (get_id_by_pnt (x, y, ident))? ident : dv;
  }
  
  void load_points_from (const char * file_name);
  
private:
 
  t_box_union_item * union_by_id (int ident) const
  {
    t_box_union_item * t;
    for (t = m_box_unions; t != NULL; t = t->tl)
      if (t->ident == ident)
        break;
    return t;
  }
 
  void load_points_from (std::istream & input);
  
private:
  size_t m_dir; 
  number m_dx; 
  number m_dy; 
  t_box_union_item * m_box_unions; 
};
 
} // namespace d3f
} // end namespace ug
 
#endif // __H__UG__PLUGINS__D3F__BOXUNIONS__
 
/* End of File */