File: BoundingSpherePy.cpp

package info (click to toggle)
esys-particle 2.1-4
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 7,284 kB
  • sloc: cpp: 77,304; python: 5,647; makefile: 1,176; sh: 10
file content (161 lines) | stat: -rw-r--r-- 4,731 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
/////////////////////////////////////////////////////////////
//                                                         //
// 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          //
//                                                         //
/////////////////////////////////////////////////////////////


#include <boost/version.hpp>
#include <boost/python.hpp>
#include <iostream>
#include <sstream>
#include "Python/esys/lsm/util/BoundingSpherePy.h"

using namespace boost::python;
using namespace esys::lsm;

namespace esys
{
  namespace lsm
  {
    BoundingSpherePy::BoundingSpherePy() : BoundingSphere()
    {
    }

    BoundingSpherePy::BoundingSpherePy(const Vec3Py &centrePt, double radius)
      : BoundingSphere(centrePt, radius)
    {
    }

    BoundingSpherePy::BoundingSpherePy(const Vec3 &centrePt, double radius)
      : BoundingSphere(centrePt, radius)
    {
    }

    BoundingSpherePy::BoundingSpherePy(const BoundingSpherePy &v)
      : BoundingSphere(v)
    {
    }

    BoundingSpherePy::BoundingSpherePy(const BoundingSphere &v)
      : BoundingSphere(v)
    {
    }

    BoundingSpherePy::BoundingSpherePy(const boost::python::object &pyCentre, double radius)
    {
      if (esys::lsm::bpu::len(pyCentre) == 3)
      {
        *this = BoundingSpherePy(Vec3Py(pyCentre), radius);
      }
      else
      {
        const std::string objectString =
          boost::python::extract<std::string>(boost::python::str(pyCentre))();
        std::stringstream msg;
        msg << "Could not extract (x,y,z) elements from: " << objectString;
        throw runtime_error(msg.str());
      }
    }

    bool BoundingSpherePy::operator==(const BoundingSpherePy &bSphere) const
    {
      return BoundingSphere::operator ==(bSphere);
    }

    Vec3Py BoundingSpherePy::getCentrePy() const
    {
      return getCentre();
    }

    class BoundingSpherePyPickleSuite : public boost::python::pickle_suite
    {
    public:
      static
      boost::python::tuple
      getinitargs(BoundingSpherePy const& bSphere)
      {
          return boost::python::make_tuple(bSphere.getCentrePy(), bSphere.getRadius());
      }
    };

    using boost::python::arg;
    void exportBoundingSphere()
    {
      // Check that Boost 1.34.0 or higher is being used.
      // If so, disable auto-generation of C++ signatures for Epydoc
      // (which stumbles over indentation in the auto-generated strings).
      #if ((BOOST_VERSION / 100000 >= 1) \
          && (BOOST_VERSION / 100 % 1000 >= 34)) \
          || (BOOST_VERSION / 100000 >= 2)
        boost::python::docstring_options no_autogen(true,false);
      #endif

      class_<esys::lsm::BoundingSpherePy>(
        "BoundingSphere",
        "A sphere - centre point and radius.",
        init<>()
      )
        .def(
          init<const BoundingSpherePy &>(
            ( arg("bSphere") )
          )
        )
        .def(
          init<const object &, double>(
            ( arg("seq"), arg("radius") )
          )
        )
        .def(
          init<const Vec3Py &, double>(
            ( arg("centrePt")=Vec3Py(0,0,0), arg("radius")=0.0 ),
            "Creates sphere with specified centre and radius.\n"
            "@type centrePt: L{Vec3}\n"
            "@kwarg centrePt: Centre of sphere.\n"
            "@type radius: float\n"
            "@kwarg radius: Radius of sphere."
          )
        )
        .def(
          "getCentre",
          &BoundingSpherePy::getCentrePy,
          "Returns centre point of this sphere.\n"
          "@rtype: L{Vec3}"
        )
        .def(
          "getCenter",
          &BoundingSpherePy::getCentrePy,
          "Returns centre point of this sphere.\n"
          "@rtype: L{Vec3}"
        )
        .def(
          "getRadius",
          &BoundingSpherePy::getRadius,
          "Returns radius of this sphere.\n"
          "@rtype: float"
        )
        .def(
          self == self//, "Tests equality of spheres - equality of centre point and radius."
        )
        .def(self_ns::str(self))
        .def_pickle(BoundingSpherePyPickleSuite())
      ;
    }
  }
}

std::ostream &operator<<(std::ostream &oStream, const esys::lsm::BoundingSpherePy &bSphere)
{
  oStream
    << bSphere.getCentre()[0]
    << " " << bSphere.getCentre()[1]
    << " " << bSphere.getCentre()[2]
    << " " << bSphere.getRadius();
  return oStream;
}