/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/geometry/geometry_plane_box.h

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_GEOMETRY_GEOMETRY_PLANE_BOX_H
#define OPENTISSUE_CORE_GEOMETRY_GEOMETRY_PLANE_BOX_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <OpenTissue/core/geometry/geometry_aabb.h>
#include <OpenTissue/core/geometry/geometry_plane.h>

namespace OpenTissue
{
  namespace geometry
  {

    template<typename math_types>
    class PlaneBox
    {
    public:

      typedef typename math_types::value_traits      value_traits;
      typedef typename math_types::vector3_type      vector3_type;
      typedef typename math_types::real_type         real_type;
      typedef Plane<math_types>                      plane_type;

    protected:

      typedef AABB<math_types>                       aabb_type;

      aabb_type   m_box;
      plane_type  m_plane;
      int         m_max_tick;       
      int         m_cur_tick;       
      int         m_axis;           


      vector3_type m_p[4][3];      
      vector3_type m_n[3];         

    public:

      plane_type   const & plane() const { return m_plane;       }
      aabb_type    const & box()   const { return m_box;         }
      vector3_type const & p0()    const { return m_p[0][m_axis];}
      vector3_type const & p1()    const { return m_p[1][m_axis];}
      vector3_type const & p2()    const { return m_p[2][m_axis];}
      vector3_type const & p3()    const { return m_p[3][m_axis];}
      vector3_type const &  n()    const { return m_n[m_axis];   }

    public:

      PlaneBox()
        : m_max_tick(100)
        , m_cur_tick(50)
        , m_axis(0)
      {
        m_n[0] = vector3_type(1,0,0);
        m_n[1] = vector3_type(0,1,0);
        m_n[2] = vector3_type(0,0,1);
      }

      void init( vector3_type const & min_coord, vector3_type const & max_coord )
      {
        m_box.set(min_coord, max_coord);
        m_cur_tick = m_max_tick/2;
        update();
      }

      void increment()
      {
        ++m_cur_tick;
        if(m_cur_tick>m_max_tick)
          m_cur_tick = m_max_tick;
        update();
      }

      void decrement()
      {
        --m_cur_tick;
        if(m_cur_tick<0)
          m_cur_tick = 0;
        update();
      }

      void set_x_axis(){ m_axis = 0; update(); }
      void set_y_axis(){ m_axis = 1; update(); }
      void set_z_axis(){ m_axis = 2; update(); }

    protected:

      void update()
      {
        real_type fraction = (value_traits::one()*m_cur_tick / m_max_tick);

        // x-base
        m_p[0][0] = m_box.min();
        m_p[1][0] = m_box.min() + vector3_type( value_traits::zero(), m_box.h(),            value_traits::zero() );
        m_p[2][0] = m_box.min() + vector3_type( value_traits::zero(), m_box.h(),            m_box.d()            );
        m_p[3][0] = m_box.min() + vector3_type( value_traits::zero(), value_traits::zero(), m_box.d()            );

        // y-base
        m_p[0][1] = m_box.min();
        m_p[1][1] = m_box.min() + vector3_type( value_traits::zero(), value_traits::zero(), m_box.d()            );
        m_p[2][1] = m_box.min() + vector3_type( m_box.w(),            value_traits::zero(), m_box.d()            );
        m_p[3][1] = m_box.min() + vector3_type( m_box.w(),            value_traits::zero(), value_traits::zero() );

        // z-base
        m_p[0][2] = m_box.min();
        m_p[1][2] = m_box.min() + vector3_type( m_box.w(),            value_traits::zero(), value_traits::zero() );
        m_p[2][2] = m_box.min() + vector3_type( m_box.w(),            m_box.h(),            value_traits::zero() );
        m_p[3][2] = m_box.min() + vector3_type( value_traits::zero(), m_box.h(),            value_traits::zero() );

        m_p[0][0](0) += m_box.w()*fraction;
        m_p[1][0](0) += m_box.w()*fraction;
        m_p[2][0](0) += m_box.w()*fraction;
        m_p[3][0](0) += m_box.w()*fraction;

        m_p[0][1](1) += m_box.h()*fraction;
        m_p[1][1](1) += m_box.h()*fraction;
        m_p[2][1](1) += m_box.h()*fraction;
        m_p[3][1](1) += m_box.h()*fraction;

        m_p[0][2](2) += m_box.d()*fraction;
        m_p[1][2](2) += m_box.d()*fraction;
        m_p[2][2](2) += m_box.d()*fraction;
        m_p[3][2](2) += m_box.d()*fraction;

        m_plane.set( m_n[m_axis], m_p[0][m_axis]);

      }

    };

  } // namespace geometry
} // namespace OpenTissue

// OPENTISSUE_CORE_GEOMETRY_GEOMETRY_PLANE_BOX_H
#endif