IntersectionVolCalculator.h

//                                                         //
// Copyright (c) 2003-2011 by The University of Queensland //
// Earth Systems Science Computational Centre (ESSCC)      //
// http://www.uq.edu.au/esscc                              //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //


#ifndef ESYS_LSM_IMPLINTERSECTIONVOLCALCULATOR_H
#define ESYS_LSM_IMPLINTERSECTIONVOLCALCULATOR_H

#include <math.h>

namespace esys
{
  namespace lsm
  {
    namespace impl
    {
      double square(double val);

      template <int tmplDim, typename TmplVec>
      class DimBasicBox
      {
      public:
        typedef TmplVec Vec;
        DimBasicBox(const Vec &minPt, const Vec &maxPt);
  
        const Vec &getMinPt() const;
  
        const Vec &getMaxPt() const;
  
        double getVolume() const;
        
        template <typename TmplSphere>
        bool intersectsWith(const TmplSphere &sphere) const;
  
        bool intersectsWith(const Vec &pt) const;
  
        template <typename TmplSphere>
        bool contains(const TmplSphere &sphere) const;
  
      private:
        Vec m_minPt;
        Vec m_maxPt;
      };
  
      template <int tmplDim, typename TmplVec>
      class DimPlane
      {
      public:
        typedef TmplVec Vec;
  
        static double norm(const Vec &pt);
  
        static double dot(const Vec &p1, const Vec &p2);
        
        DimPlane();
  
        DimPlane(const Vec &normal, const Vec &pt);
  
        DimPlane(const DimPlane &plane);
  
        DimPlane &operator=(const DimPlane &plane);
  
        double getSignedDistanceTo(const Vec &pt) const;
  
        double getDistanceTo(const Vec &pt) const;
  
        const Vec &getNormal() const;
  
      private:
        Vec    m_normal;
        Vec    m_pt;
        double m_invNormalNorm;
      };
  
      template <int tmplDim, typename TmplVec>
      class DimBasicSphere
      {
      public:
        typedef TmplVec                 Vec;
        typedef DimPlane<tmplDim, Vec> Plane;
  
        static const double FOUR_THIRDS_PI;
        static const double ONE_THIRD_PI;

        DimBasicSphere();
        
        DimBasicSphere(const Vec &centrePt, double radius);

        DimBasicSphere(const DimBasicSphere &sphere);

        DimBasicSphere &operator=(const DimBasicSphere &sphere);

        double getRadius() const;

        const Vec &getCentre() const;

        double getVolume() const;

        double getVolume(const Vec &minPt, const Vec &maxPt,DimBasicSphere();