-
Notifications
You must be signed in to change notification settings - Fork 0
/
vector.h
111 lines (94 loc) · 2.08 KB
/
vector.h
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
#ifndef VECTOR_H_
#define VECTOR_H_
#include <iostream>
#include <memory>
#include <cstring>
template<class R>
class Vector {
public:
Vector() {
#ifndef NDEBUG
for (int i=0; i<R::Dim(); ++i) {
data_[i] = 0xDEADBEEF;
}
#endif // NDEBUG
}
inline int& operator[](int i) {
return data_[i];
}
inline int operator[](int i) const {
return data_[i];
}
inline bool operator==(const Vector<R>& rhs) const {
for (int i=0; i<R::Dim(); ++i) {
if ((*this)[i] != rhs[i]) return false;
}
return true;
}
inline bool operator!=(const Vector<R>& rhs) const {
for (int i=0; i<R::Dim(); ++i) {
if ((*this)[i] != rhs[i]) return true;
}
return false;
}
inline bool IsWithinDeviation() const {
for (int i=0; i<R::Dim(); ++i) {
if (0 > (*this)[i] || (*this)[i] >= R::Deviation(i)) return false;
}
return true;
}
inline int L1Norm() const {
int norm = 0;
for (int i=0; i<R::Dim(); ++i) {
norm += std::abs((*this)[i]);
}
return norm;
}
inline Vector<R>& operator+=(const Vector<R>& rhs) {
for (int i=0; i<R::Dim(); ++i) {
(*this)[i] += rhs[i];
}
return *this;
}
inline Vector<R>& operator-=(const Vector<R>& rhs) {
for (int i=0; i<R::Dim(); ++i) {
(*this)[i] -= rhs[i];
}
return *this;
}
inline double Volume() const {
double v = 0.0;
for (int i=0; i<R::Dim(); ++i) {
v += R::Size(i) * (*this)[i];
}
return v;
}
inline void Reset() {
std::memset(data_, 0, sizeof(data_));
}
inline bool IsPositive() const {
for (int i = 0; i < R::Dim(); ++i) {
if ((*this)[i] < 0) return false;
}
return true;
}
inline bool IsZero() const {
for (int i = 0; i < R::Dim(); ++i) {
if ((*this)[i] != 0) return false;
}
return true;
}
private:
int data_[R::Dim()];
};
template<class R>
std::ostream& operator<<(std::ostream& os, const Vector<R>& obj) {
os << "(";
for (int i=0; i<R::Dim(); ++i) {
os << obj[i];
if (i < R::Dim() - 1) os << ", ";
}
os << ")";
return os;
}
#endif // VECTOR_H_