geometric_shapes.cpp 5.03 KB
Newer Older
sachinc's avatar
sachinc committed
1
2
3
/*
 * Software License Agreement (BSD License)
 *
4
5
 *  Copyright (c) 2011-2014, Willow Garage, Inc.
 *  Copyright (c) 2014-2015, Open Source Robotics Foundation
sachinc's avatar
sachinc committed
6
7
8
9
10
11
12
13
14
15
16
17
 *  All rights reserved.
 *
 *  Redistribution and use in source and binary forms, with or without
 *  modification, are permitted provided that the following conditions
 *  are met:
 *
 *   * Redistributions of source code must retain the above copyright
 *     notice, this list of conditions and the following disclaimer.
 *   * Redistributions in binary form must reproduce the above
 *     copyright notice, this list of conditions and the following
 *     disclaimer in the documentation and/or other materials provided
 *     with the distribution.
18
 *   * Neither the name of Open Source Robotics Foundation nor the names of its
sachinc's avatar
sachinc committed
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
 *     contributors may be used to endorse or promote products derived
 *     from this software without specific prior written permission.
 *
 *  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 *  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 *  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
 *  FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
 *  COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
 *  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
 *  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
 *  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
 *  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
 *  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
 *  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 *  POSSIBILITY OF SUCH DAMAGE.
 */

/** \author Jia Pan */


39
40
#include <hpp/fcl/shape/geometric_shapes.h>
#include <hpp/fcl/shape/geometric_shapes_utility.h>
sachinc's avatar
sachinc committed
41

42
43
namespace hpp
{
sachinc's avatar
sachinc committed
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
namespace fcl
{

void Convex::fillEdges()
{
  int* points_in_poly = polygons;
  if(edges) delete [] edges;

  int num_edges_alloc = 0;
  for(int i = 0; i < num_planes; ++i)
  {
    num_edges_alloc += *points_in_poly;
    points_in_poly += (*points_in_poly + 1);
  }

  edges = new Edge[num_edges_alloc];

  points_in_poly = polygons;
  int* index = polygons + 1;
  num_edges = 0;
  Edge e;
  bool isinset;
  for(int i = 0; i < num_planes; ++i)
  {
    for(int j = 0; j < *points_in_poly; ++j)
    {
      e.first = std::min(index[j], index[(j+1)%*points_in_poly]);
      e.second = std::max(index[j], index[(j+1)%*points_in_poly]);
      isinset = false;
      for(int k = 0; k < num_edges; ++k)
      {
        if((edges[k].first == e.first) && (edges[k].second == e.second))
        {
          isinset = true;
          break;
        }
      }

      if(!isinset)
      {
        edges[num_edges].first = e.first;
        edges[num_edges].second = e.second;
        ++num_edges;
      }
    }

    points_in_poly += (*points_in_poly + 1);
    index = points_in_poly + 1;
  }

  if(num_edges < num_edges_alloc)
  {
    Edge* tmp = new Edge[num_edges];
    memcpy(tmp, edges, num_edges * sizeof(Edge));
    delete [] edges;
    edges = tmp;
  }
}

103
104
void Halfspace::unitNormalTest()
{
105
  FCL_REAL l = n.norm();
106
107
108
109
110
111
112
113
  if(l > 0)
  {
    FCL_REAL inv_l = 1.0 / l;
    n *= inv_l;
    d *= inv_l;
  }
  else
  {
114
    n << 1, 0, 0;
115
116
117
118
    d = 0;
  }  
}

sachinc's avatar
sachinc committed
119
120
void Plane::unitNormalTest()
{
121
  FCL_REAL l = n.norm();
sachinc's avatar
sachinc committed
122
123
  if(l > 0)
  {
124
    FCL_REAL inv_l = 1.0 / l;
125
    n *= inv_l;
sachinc's avatar
sachinc committed
126
127
128
129
    d *= inv_l;
  }
  else
  {
130
    n << 1, 0, 0;
sachinc's avatar
sachinc committed
131
132
133
134
    d = 0;
  }
}

135
136
137

void Box::computeLocalAABB()
{
138
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
139
  aabb_center = aabb_local.center();
140
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
141
142
143
144
}

void Sphere::computeLocalAABB()
{
145
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
146
  aabb_center = aabb_local.center();
147
148
149
150
151
  aabb_radius = radius;
}

void Capsule::computeLocalAABB()
{
152
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
153
  aabb_center = aabb_local.center();
154
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
155
156
157
158
}

void Cone::computeLocalAABB()
{
159
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
160
  aabb_center = aabb_local.center();
161
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
162
163
164
165
}

void Cylinder::computeLocalAABB()
{
166
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
167
  aabb_center = aabb_local.center();
168
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
169
170
171
172
}

void Convex::computeLocalAABB()
{
173
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
174
  aabb_center = aabb_local.center();
175
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
176
177
}

178
179
180
181
void Halfspace::computeLocalAABB()
{
  computeBV<AABB>(*this, Transform3f(), aabb_local);
  aabb_center = aabb_local.center();
182
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
183
184
}

185
186
void Plane::computeLocalAABB()
{
187
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
188
  aabb_center = aabb_local.center();
189
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
190
191
}

192
void TriangleP::computeLocalAABB()
193
{
194
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
195
  aabb_center = aabb_local.center();
196
  aabb_radius = (aabb_local.min_ - aabb_center).norm();
197
198
199
}


sachinc's avatar
sachinc committed
200
}
201
202

} // namespace hpp