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
|
#include "point.H"
#include "triangle.H"
#include "tetrahedron.H"
#include "IOstreams.H"
using namespace Foam;
int main()
{
triangle<point, point> tri
(
vector(0, 0, 0),
vector(1, 0, 0),
vector(1, 1, 0)
);
Info<< "tri circumCentre = " << tri.circumCentre() << endl;
Info<< "tri circumRadius = " << tri.circumRadius() << endl;
tetrahedron<point, point> tet
(
vector(1, 0, 0),
vector(0, 1, 0),
vector(0, 0, 1),
vector(0.5773502, 0.5773502, 0.5773502)
);
Info<< "tet circumCentre = " << tet.circumCentre() << endl;
Info<< "tet circumRadius = " << tet.circumRadius() << endl;
vector a(Sin);
vector b(Sin);
vector c(Sin);
vector d(Sin);
Info<< "tet circumRadius = "
<< tetrahedron<point, point>(a, b, c, d).circumRadius() << endl;
return 0;
}
|
