/usr/include/libwildmagic/Wm5Tetrahedron3.inl is in libwildmagic-dev 5.13-1+b2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | // Geometric Tools, LLC
// Copyright (c) 1998-2014
// Distributed under the Boost Software License, Version 1.0.
// http://www.boost.org/LICENSE_1_0.txt
// http://www.geometrictools.com/License/Boost/LICENSE_1_0.txt
//
// File Version: 5.0.0 (2010/01/01)
//----------------------------------------------------------------------------
template <typename Real>
Tetrahedron3<Real>::Tetrahedron3 ()
{
}
//----------------------------------------------------------------------------
template <typename Real>
Tetrahedron3<Real>::~Tetrahedron3 ()
{
}
//----------------------------------------------------------------------------
template <typename Real>
Tetrahedron3<Real>::Tetrahedron3 (const Vector3<Real>& v0,
const Vector3<Real>& v1, const Vector3<Real>& v2,
const Vector3<Real>& v3)
{
V[0] = v0;
V[1] = v1;
V[2] = v2;
V[3] = v3;
}
//----------------------------------------------------------------------------
template <typename Real>
Tetrahedron3<Real>::Tetrahedron3 (const Vector3<Real> vertex[4])
{
V[0] = vertex[0];
V[1] = vertex[1];
V[2] = vertex[2];
V[3] = vertex[3];
}
//----------------------------------------------------------------------------
template <typename Real>
void Tetrahedron3<Real>::GetFaceIndices (int face, int index[3]) const
{
// assert: 0 <= face <= 3
if (face == 0)
{
index[0] = 0;
index[1] = 2;
index[2] = 1;
}
else if (face == 1)
{
index[0] = 0;
index[1] = 1;
index[2] = 3;
}
else if (face == 2)
{
index[0] = 0;
index[1] = 3;
index[2] = 2;
}
else // face == 3 (no index validation is performed)
{
index[0] = 1;
index[1] = 2;
index[2] = 3;
}
}
//----------------------------------------------------------------------------
template <typename Real>
void Tetrahedron3<Real>::GetPlanes (Plane3<Real> plane[4]) const
{
Vector3<Real> edge10 = V[1] - V[0];
Vector3<Real> edge20 = V[2] - V[0];
Vector3<Real> edge30 = V[3] - V[0];
Vector3<Real> edge21 = V[2] - V[1];
Vector3<Real> edge31 = V[3] - V[1];
plane[0].Normal = edge20.UnitCross(edge10); // <v0,v2,v1>
plane[1].Normal = edge10.UnitCross(edge30); // <v0,v1,v3>
plane[2].Normal = edge30.UnitCross(edge20); // <v0,v3,v2>
plane[3].Normal = edge21.UnitCross(edge31); // <v1,v2,v3>
Real det = edge10.Dot(plane[3].Normal);
if (det < (Real)0)
{
// The normals are inner pointing, reverse their directions.
for (int i = 0; i < 4; ++i)
{
plane[i].Normal = -plane[i].Normal;
}
}
for (int i = 0; i < 4; ++i)
{
plane[i].Constant = V[i].Dot(plane[i].Normal);
}
}
//----------------------------------------------------------------------------
|