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
|
#include "mmvec.hpp"
#include <iostream>
#include <complex>
int main(int argc, char *argv[]) {
// N dimensional vectors
std::cout << "N dimensional (int) vectors" << std::endl;
mm::vec<int, 5> u {3, 2, 1, 0, 1};
mm::vec<int, 5> v {1, 2, 3, 4, 5};
std::cout << "u = " << u << std::endl;
std::cout << "v = " << v << std::endl;
std::cout << std::endl;
// basic operations
std::cout << "u + v = " << u + v << std::endl;
std::cout << "u - v = " << u - v << std::endl;
std::cout << "2 * v = " << 2 * v << std::endl;
std::cout << "v * 2 = " << v * 2 << std::endl;
std::cout << "u * v = " << u * v << std::endl;
std::cout << std::endl;
u += v;
// three dimensional vectors
std::cout << "three dimensional (double) vectors" << std::endl;
mm::vec3<double> a {1, 2, 3};
mm::vec3<double> b {3, 2, 1};
std::cout << "a = " << a << std::endl;
std::cout << "b = " << b << std::endl;
std::cout << std::endl;
std::cout << "a x b = " << mm::vec3<double>::cross(a, b) << std::endl;
std::cout << "zenith(a) = " << a.zenith() << std::endl;
std::cout << "azimuth(a) = " << a.azimuth() << std::endl;
std::cout << "spherical(a) = " << a.spherical() << std::endl;
std::cout << std::endl;
a -= b;
// two dimensional vector
std::cout << "two dimensional (complex) vectors" << std::endl;
mm::vec2<std::complex<float>> j {{1, 2}, {3, -1}};
mm::vec2<std::complex<float>> k { 5, {-2, 1}};
std::cout << "j = " << j << std::endl;
std::cout << "k = " << k << std::endl;
std::cout << std::endl;
std::cout << "j x k = " << mm::vec2<std::complex<float>>::cross(j, k) << std::endl;
std::cout << "angle(j) = " << j.angle() << std::endl;
std::cout << "polar(j) = " << j.polar() << std::endl;
std::cout << std::endl;
j *= 10;
return 0;
}
|
