structured_block.h

Go to the documentation of this file.

//
// Copyright (c) 2014 Limit Point Systems, Inc. 
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//   http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.
//


#ifndef STRUCTURED_BLOCK_H
#define STRUCTURED_BLOCK_H

#ifndef SHEAF_DLL_SPEC_H
#include "SheafSystem/sheaf_dll_spec.h"
#endif

#ifndef HOMOGENEOUS_BLOCK_H
#include "SheafSystem/homogeneous_block.h"
#endif

#ifndef CHART_POINT_H
#include "SheafSystem/chart_point.h"
#endif

#ifndef STD_CMATH_H
#include "SheafSystem/std_cmath.h"
#endif

namespace fiber_bundle
{
  
using namespace sheaf;

class chart_point_3d;

class SHEAF_DLL_SPEC structured_block : public homogeneous_block
{
  // ===========================================================
  // ===========================================================

public:

  static host_type& new_host(namespace_type& xns, 
                 const poset_path& xhost_path, 
                 const poset_path& xschema_path, 
                 int xmax_db,
                 bool xauto_access);

  static host_type& new_host(namespace_type& xns,
                 const poset_path& xhost_path,
                 int xmax_db,
                 bool xauto_access);
  

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  structured_block();

  virtual ~structured_block();

protected:

  virtual void init_handle_data_members();

  index_space_handle* _elements_id_space;

  void refine_coordinate(chart_point_coord_type& xcoord,
                         scoped_index::pod_type& xchart_id,
                         size_type xchart_id_ub) const;

  static inline double TOLERANCE()
  {
    static const double result = 10.0*std::numeric_limits<double>::epsilon();
    return result;
  };

  static inline bool is_close_enough(double x1, double x2)
  {
    return (std::fabs(x2 - x1) <= TOLERANCE());
  };

private:



  // ===========================================================
  // ===========================================================

public:


  chart_point* refine_point(const chart_point& xpt) const;

  virtual void refine_point_pa(const chart_point& xpt, chart_point& result) const = 0;

  virtual void refine_point_pa(const chart_point& xpt,
                               block<chart_point_3d>& result) const = 0;

  chart_point* unrefine_point(const chart_point& xpt) const;

  virtual void unrefine_point_pa(const chart_point& xpt, chart_point& result) const = 0;

  virtual void refine() = 0;

  virtual void unrefine() = 0;

  virtual void simplify() = 0;

  virtual void unsimplify() = 0;

protected:

private:



  // ===========================================================
  // ===========================================================

public:

protected:

private:



  // ===========================================================
  // ===========================================================

public:

protected:

private:



  // ===========================================================
  // ===========================================================

public:

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  virtual structured_block* clone() const = 0;

  inline structured_block* clone(bool xnew_state, bool xauto_access = true) const
  {
    return static_cast<structured_block*>(homogeneous_block::clone(xnew_state, xauto_access));
  }

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  inline structured_block* greatest_jem() const
  {
    return static_cast<structured_block*>(base_space_member::greatest_jem());
  };

  inline structured_block* least_jem() const
  {
    return static_cast<structured_block*>(base_space_member::least_jem());
  };

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  inline structured_block* p_join(abstract_poset_member* other) const
  {
    return static_cast<structured_block*>(base_space_member::p_join(other));
  };

  inline structured_block* p_meet(abstract_poset_member* other)
  {
    return static_cast<structured_block*>(base_space_member::p_meet(other));
  };

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  inline structured_block* l_join(abstract_poset_member* other, bool xnew_jem = true)
  {
    return static_cast<structured_block*>(base_space_member::l_join(other, xnew_jem));
  };

  inline structured_block* l_meet(abstract_poset_member* other, bool xnew_jem = true)
  {
    return static_cast<structured_block*>(base_space_member::l_meet(other, xnew_jem));
  };

  inline structured_block* l_not(bool xnew_jem = true) const
  {
    return static_cast<structured_block*>(base_space_member::l_not(xnew_jem));
  };

protected:

private:



  // ===========================================================
  // ===========================================================

public:

  virtual bool invariant() const;

  virtual bool is_ancestor_of(const any* other) const;

  virtual structured_block& operator=(const abstract_poset_member& xother);

  structured_block& operator=(const structured_block& xother);

protected:

private:

};

//==============================================================================
// NON-MEMBER FUNCTIONS
//==============================================================================

SHEAF_DLL_SPEC size_t deep_size(const structured_block& x0, bool xinclude_shallow = true);

 
} // namespace fiber_bundle


#endif // STRUCTURED_BLOCK_H