NumericalPhysics/vec.h at master · cznull/NumericalPhysics · GitHub

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
template <int n>
struct doublen
{
	double x[n];
};

template <int n>
double dot(doublen<n> x1, doublen<n> x2);

template <int n>
double norm(doublen<n> x1);

template <int n>
doublen<n> operator-(doublen<n> x1, doublen<n> x2);

template <int n>
doublen<n> operator+(doublen<n> x1, doublen<n> x2);

struct double2 {
	double x, y;
};

double dot(double2 x1, double2 x2);

double norm(double2 x1);

double vol(double2 x1, double2 x2);

double2 operator-(double2 x1, double2 x2);

double2 operator+(double2 x1, double2 x2);

double2 operator*(double s, double2 x);

double2 operator*(double2 x, double s);

struct double3 {
	double x, y, z;
};

double dot(double3 x1, double3 x2);

double norm(double3 x1);

double vol(double3 x1, double3 x2, double3 x3);

double3 cross(double3 x1, double3 x2);

double3 operator-(double3 x1, double3 x2);

double3 operator+(double3 x1, double3 x2);

double3 operator*(double3 x, double s);

struct float3 {
	float x, y, z;
	float3(float a, float b, float c) :x(a),y(b),z(c){
	}
	float3(double a, double b, double c) :x(a), y(b), z(c) {
	}

	float3(double3 vertex) {
		x = vertex.x;
		y = vertex.y;
		z = vertex.z;
	}
	float3() {
	}
};

float dot(float3 x1, float3 x2);

float norm(float3 x1);

float vol(float3 x1, float3 x2, float3 x3);

float3 cross(float3 x1, float3 x2);

float3 operator-(float3 x1, float3 x2);

float3 operator+(float3 x1, float3 x2);

float3 operator*(float3 x, float s);