/usr/include/libwildmagic/Wm5Torus3.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 100 101 102 103 104 105 106 107 108 109 110 | // 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>
Torus3<Real>::Torus3 ()
{
}
//----------------------------------------------------------------------------
template <typename Real>
Torus3<Real>::~Torus3 ()
{
}
//----------------------------------------------------------------------------
template <typename Real>
Torus3<Real>::Torus3 (Real outerRadius, Real innerRadius)
{
// assert: outerRadius > 0 && innerRadius > 0
OuterRadius = outerRadius;
InnerRadius = innerRadius;
}
//----------------------------------------------------------------------------
template <typename Real>
Vector3<Real> Torus3<Real>::GetPosition (Real s, Real t) const
{
Real twoPiS = Math<Real>::TWO_PI*s;
Real twoPiT = Math<Real>::TWO_PI*t;
Real cosTwoPiS = Math<Real>::Cos(twoPiS);
Real sinTwoPiS = Math<Real>::Sin(twoPiS);
Real cosTwoPiT = Math<Real>::Cos(twoPiT);
Real sinTwoPiT = Math<Real>::Sin(twoPiT);
Real maxRadius = OuterRadius + InnerRadius*cosTwoPiT;
Vector3<Real> position(
maxRadius*cosTwoPiS,
maxRadius*sinTwoPiS,
InnerRadius*sinTwoPiT);
return position;
}
//----------------------------------------------------------------------------
template <typename Real>
Vector3<Real> Torus3<Real>::GetNormal (Real s, Real t) const
{
Real twoPiS = Math<Real>::TWO_PI*s;
Real cosTwoPiS = Math<Real>::Cos(twoPiS);
Real sinTwoPiS = Math<Real>::Sin(twoPiS);
Vector3<Real> position = GetPosition(s, t);
Vector3<Real> normal(
position.X() - OuterRadius*cosTwoPiS,
position.Y() - OuterRadius*sinTwoPiS,
position.Z());
normal.Normalize();
return normal;
}
//----------------------------------------------------------------------------
template <typename Real>
void Torus3<Real>::GetParameters (const Vector3<Real>& position, Real& s,
Real& t) const
{
Real radius = Math<Real>::Sqrt(position.X()*position.X() +
position.Y()*position.Y());
Real angle;
if (radius < Math<Real>::ZERO_TOLERANCE)
{
s = (Real)0;
}
else
{
angle = Math<Real>::ATan2(position.Y(), position.X());
if (angle >= (Real)0)
{
s = angle*Math<Real>::INV_TWO_PI;
}
else
{
s = (Real)1 + angle*Math<Real>::INV_TWO_PI;
}
}
Real diff = radius - OuterRadius;
if (Math<Real>::FAbs(diff) < Math<Real>::ZERO_TOLERANCE
&& Math<Real>::FAbs(position.Z()) < Math<Real>::ZERO_TOLERANCE)
{
t = (Real)0;
}
else
{
angle = Math<Real>::ATan2(position.Z(), diff);
if (angle >= (Real)0)
{
t = angle*Math<Real>::INV_TWO_PI;
}
else
{
t = (Real)1 + angle*Math<Real>::INV_TWO_PI;
}
}
}
//----------------------------------------------------------------------------
|