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
|
/**************************************************************************/
/* Copyright 2009 Tim Day */
/* */
/* This file is part of Fracplanet */
/* */
/* Fracplanet is free software: you can redistribute it and/or modify */
/* it under the terms of the GNU General Public License as published by */
/* the Free Software Foundation, either version 3 of the License, or */
/* (at your option) any later version. */
/* */
/* Fracplanet is distributed in the hope that it will be useful, */
/* but WITHOUT ANY WARRANTY; without even the implied warranty of */
/* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the */
/* GNU General Public License for more details. */
/* */
/* You should have received a copy of the GNU General Public License */
/* along with Fracplanet. If not, see <http://www.gnu.org/licenses/>. */
/**************************************************************************/
#include "xyz.h"
/*! Table so we can look up by element number.
*/
XYZ::ElementPtr XYZ::element_table[3]={&XYZ::x,&XYZ::y,&XYZ::z};
std::ostream& XYZ::write(std::ostream& out) const
{
return out << x << " " << y << " " << z;
}
const std::string XYZ::format_comma() const
{
std::ostringstream s;
s << x << "," << y << "," << z;
return s.str();
}
const std::string XYZ::format_blender() const
{
std::ostringstream s;
s << x << ", " << y << ", " << z;
return s.str();
}
/*! Also transposes y and z co-ordinates because of POV-Rays idea of up.
*/
const std::string XYZ::format_pov() const
{
std::ostringstream s;
s << "<" << x << "," << z << "," << y << ">";
return s.str();
}
RandomXYZInUnitCube::RandomXYZInUnitCube(Random01& rng)
:XYZ()
{
x=rng();
y=rng();
z=rng();
}
RandomXYZInBox::RandomXYZInBox(Random01& rng,const XYZ& bounds)
:XYZ()
{
x=-bounds.x+2.0*bounds.x*rng();
y=-bounds.y+2.0*bounds.y*rng();
z=-bounds.z+2.0*bounds.z*rng();
}
RandomXYZInSphere::RandomXYZInSphere(Random01& rng,float radius)
:XYZ()
{
do
{
x=2.0*rng()-1.0;
y=2.0*rng()-1.0;
z=2.0*rng()-1.0;
}
while (magnitude2()>1.0);
(*this)*=radius;
}
/*! Can handle case of individual axes being zero.
*/
RandomXYZInEllipsoid::RandomXYZInEllipsoid(Random01& rng,const XYZ& axes)
:XYZ()
{
do
{
assign(RandomXYZInBox(rng,axes));
}
while (
(axes.x==0.0 ? 0.0 : sqr(x/axes.x))
+(axes.y==0.0 ? 0.0 : sqr(y/axes.y))
+(axes.z==0.0 ? 0.0 : sqr(z/axes.z))
>1.0
);
}
RandomXYZSphereNormal::RandomXYZSphereNormal(Random01& rng)
:XYZ(0.0,0.0,0.0)
{
float m2;
do
{
assign(RandomXYZInSphere(rng,1.0));
m2=magnitude2();
}
while (m2==0.0);
(*this)/=sqrtf(m2);
}
|
