geometric_shapes.cpp 4.57 KB
Newer Older
sachinc's avatar
sachinc committed
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
/*
 * Software License Agreement (BSD License)
 *
 *  Copyright (c) 2011, Willow Garage, Inc.
 *  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.
 *   * Neither the name of Willow Garage, Inc. nor the names of its
 *     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 */


jpan's avatar
jpan committed
38
39
#include "fcl/shape/geometric_shapes.h"
#include "fcl/shape/geometric_shapes_utility.h"
sachinc's avatar
sachinc committed
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

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;
  }
}

void Plane::unitNormalTest()
{
102
  FCL_REAL l = n.length();
sachinc's avatar
sachinc committed
103
104
  if(l > 0)
  {
105
    FCL_REAL inv_l = 1.0 / l;
106
    n *= inv_l;
sachinc's avatar
sachinc committed
107
108
109
110
    d *= inv_l;
  }
  else
  {
111
    n.setValue(1, 0, 0);
sachinc's avatar
sachinc committed
112
113
114
115
    d = 0;
  }
}

116
117
118

void Box::computeLocalAABB()
{
119
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
120
121
  aabb_center = aabb_local.center();
  aabb_radius = (aabb_local.min_ - aabb_center).length();
122
123
124
125
}

void Sphere::computeLocalAABB()
{
126
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
127
  aabb_center = aabb_local.center();
128
129
130
131
132
  aabb_radius = radius;
}

void Capsule::computeLocalAABB()
{
133
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
134
135
  aabb_center = aabb_local.center();
  aabb_radius = (aabb_local.min_ - aabb_center).length();
136
137
138
139
}

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

void Cylinder::computeLocalAABB()
{
147
  computeBV<AABB>(*this, Transform3f(), aabb_local);
jpan's avatar
jpan committed
148
149
  aabb_center = aabb_local.center();
  aabb_radius = (aabb_local.min_ - aabb_center).length();
150
151
152
153
}

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

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

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


sachinc's avatar
sachinc committed
174
}