| 68 |
|
RectMatrix() { |
| 69 |
|
for (unsigned int i = 0; i < Row; i++) |
| 70 |
|
for (unsigned int j = 0; j < Col; j++) |
| 71 |
< |
data_[i][j] = 0.0; |
| 71 |
> |
this->data_[i][j] = 0.0; |
| 72 |
|
} |
| 73 |
|
|
| 74 |
|
/** Constructs and initializes every element of this matrix to a scalar */ |
| 75 |
|
RectMatrix(Real s) { |
| 76 |
|
for (unsigned int i = 0; i < Row; i++) |
| 77 |
|
for (unsigned int j = 0; j < Col; j++) |
| 78 |
< |
data_[i][j] = s; |
| 78 |
> |
this->data_[i][j] = s; |
| 79 |
|
} |
| 80 |
|
|
| 81 |
|
RectMatrix(Real* array) { |
| 82 |
|
for (unsigned int i = 0; i < Row; i++) |
| 83 |
|
for (unsigned int j = 0; j < Col; j++) |
| 84 |
< |
data_[i][j] = array[i * Row + j]; |
| 84 |
> |
this->data_[i][j] = array[i * Row + j]; |
| 85 |
|
} |
| 86 |
|
|
| 87 |
|
/** copy constructor */ |
| 99 |
|
|
| 100 |
|
for (unsigned int i = 0; i < Row; i++) |
| 101 |
|
for (unsigned int j = 0; j < Col; j++) |
| 102 |
< |
data_[i][j] = m.data_[i][j]; |
| 102 |
> |
this->data_[i][j] = m.data_[i][j]; |
| 103 |
|
return *this; |
| 104 |
|
} |
| 105 |
|
|
| 111 |
|
*/ |
| 112 |
|
Real& operator()(unsigned int i, unsigned int j) { |
| 113 |
|
//assert( i < Row && j < Col); |
| 114 |
< |
return data_[i][j]; |
| 114 |
> |
return this->data_[i][j]; |
| 115 |
|
} |
| 116 |
|
|
| 117 |
|
/** |
| 122 |
|
*/ |
| 123 |
|
Real operator()(unsigned int i, unsigned int j) const { |
| 124 |
|
|
| 125 |
< |
return data_[i][j]; |
| 125 |
> |
return this->data_[i][j]; |
| 126 |
|
} |
| 127 |
|
|
| 128 |
|
/** |
| 132 |
|
void getArray(Real* array) { |
| 133 |
|
for (unsigned int i = 0; i < Row; i++) { |
| 134 |
|
for (unsigned int j = 0; j < Col; j++) { |
| 135 |
< |
array[i * Row + j] = data_[i][j]; |
| 135 |
> |
array[i * Row + j] = this->data_[i][j]; |
| 136 |
|
} |
| 137 |
|
} |
| 138 |
|
} |
| 140 |
|
|
| 141 |
|
/** Returns the pointer of internal array */ |
| 142 |
|
Real* getArrayPointer() { |
| 143 |
< |
return &data_[0][0]; |
| 143 |
> |
return &this->data_[0][0]; |
| 144 |
|
} |
| 145 |
|
|
| 146 |
|
/** |
| 152 |
|
Vector<Real, Row> v; |
| 153 |
|
|
| 154 |
|
for (unsigned int i = 0; i < Row; i++) |
| 155 |
< |
v[i] = data_[row][i]; |
| 155 |
> |
v[i] = this->data_[row][i]; |
| 156 |
|
|
| 157 |
|
return v; |
| 158 |
|
} |
| 165 |
|
void setRow(unsigned int row, const Vector<Real, Row>& v) { |
| 166 |
|
|
| 167 |
|
for (unsigned int i = 0; i < Row; i++) |
| 168 |
< |
data_[row][i] = v[i]; |
| 168 |
> |
this->data_[row][i] = v[i]; |
| 169 |
|
} |
| 170 |
|
|
| 171 |
|
/** |
| 177 |
|
Vector<Real, Col> v; |
| 178 |
|
|
| 179 |
|
for (unsigned int j = 0; j < Col; j++) |
| 180 |
< |
v[j] = data_[j][col]; |
| 180 |
> |
v[j] = this->data_[j][col]; |
| 181 |
|
|
| 182 |
|
return v; |
| 183 |
|
} |
| 190 |
|
void setColumn(unsigned int col, const Vector<Real, Col>& v){ |
| 191 |
|
|
| 192 |
|
for (unsigned int j = 0; j < Col; j++) |
| 193 |
< |
data_[j][col] = v[j]; |
| 193 |
> |
this->data_[j][col] = v[j]; |
| 194 |
|
} |
| 195 |
|
|
| 196 |
|
/** |
| 202 |
|
assert(i < Row && j < Row); |
| 203 |
|
|
| 204 |
|
for (unsigned int k = 0; k < Col; k++) |
| 205 |
< |
std::swap(data_[i][k], data_[j][k]); |
| 205 |
> |
std::swap(this->data_[i][k], this->data_[j][k]); |
| 206 |
|
} |
| 207 |
|
|
| 208 |
|
/** |
| 214 |
|
assert(i < Col && j < Col); |
| 215 |
|
|
| 216 |
|
for (unsigned int k = 0; k < Row; k++) |
| 217 |
< |
std::swap(data_[k][i], data_[k][j]); |
| 217 |
> |
std::swap(this->data_[k][i], this->data_[k][j]); |
| 218 |
|
} |
| 219 |
|
|
| 220 |
|
/** |
| 227 |
|
bool operator ==(const RectMatrix<Real, Row, Col>& m) { |
| 228 |
|
for (unsigned int i = 0; i < Row; i++) |
| 229 |
|
for (unsigned int j = 0; j < Col; j++) |
| 230 |
< |
if (!equal(data_[i][j], m.data_[i][j])) |
| 230 |
> |
if (!equal(this->data_[i][j], m.data_[i][j])) |
| 231 |
|
return false; |
| 232 |
|
|
| 233 |
|
return true; |
| 246 |
|
inline void negate() { |
| 247 |
|
for (unsigned int i = 0; i < Row; i++) |
| 248 |
|
for (unsigned int j = 0; j < Col; j++) |
| 249 |
< |
data_[i][j] = -data_[i][j]; |
| 249 |
> |
this->data_[i][j] = -this->data_[i][j]; |
| 250 |
|
} |
| 251 |
|
|
| 252 |
|
/** |
| 256 |
|
inline void negate(const RectMatrix<Real, Row, Col>& m) { |
| 257 |
|
for (unsigned int i = 0; i < Row; i++) |
| 258 |
|
for (unsigned int j = 0; j < Col; j++) |
| 259 |
< |
data_[i][j] = -m.data_[i][j]; |
| 259 |
> |
this->data_[i][j] = -m.data_[i][j]; |
| 260 |
|
} |
| 261 |
|
|
| 262 |
|
/** |
| 266 |
|
inline void add( const RectMatrix<Real, Row, Col>& m ) { |
| 267 |
|
for (unsigned int i = 0; i < Row; i++) |
| 268 |
|
for (unsigned int j = 0; j < Col; j++) |
| 269 |
< |
data_[i][j] += m.data_[i][j]; |
| 269 |
> |
this->data_[i][j] += m.data_[i][j]; |
| 270 |
|
} |
| 271 |
|
|
| 272 |
|
/** |
| 277 |
|
inline void add( const RectMatrix<Real, Row, Col>& m1, const RectMatrix<Real, Row, Col>& m2 ) { |
| 278 |
|
for (unsigned int i = 0; i < Row; i++) |
| 279 |
|
for (unsigned int j = 0; j < Col; j++) |
| 280 |
< |
data_[i][j] = m1.data_[i][j] + m2.data_[i][j]; |
| 280 |
> |
this->data_[i][j] = m1.data_[i][j] + m2.data_[i][j]; |
| 281 |
|
} |
| 282 |
|
|
| 283 |
|
/** |
| 287 |
|
inline void sub( const RectMatrix<Real, Row, Col>& m ) { |
| 288 |
|
for (unsigned int i = 0; i < Row; i++) |
| 289 |
|
for (unsigned int j = 0; j < Col; j++) |
| 290 |
< |
data_[i][j] -= m.data_[i][j]; |
| 290 |
> |
this->data_[i][j] -= m.data_[i][j]; |
| 291 |
|
} |
| 292 |
|
|
| 293 |
|
/** |
| 298 |
|
inline void sub( const RectMatrix<Real, Row, Col>& m1, const RectMatrix<Real, Row, Col>& m2){ |
| 299 |
|
for (unsigned int i = 0; i < Row; i++) |
| 300 |
|
for (unsigned int j = 0; j < Col; j++) |
| 301 |
< |
data_[i][j] = m1.data_[i][j] - m2.data_[i][j]; |
| 301 |
> |
this->data_[i][j] = m1.data_[i][j] - m2.data_[i][j]; |
| 302 |
|
} |
| 303 |
|
|
| 304 |
|
/** |
| 308 |
|
inline void mul( Real s ) { |
| 309 |
|
for (unsigned int i = 0; i < Row; i++) |
| 310 |
|
for (unsigned int j = 0; j < Col; j++) |
| 311 |
< |
data_[i][j] *= s; |
| 311 |
> |
this->data_[i][j] *= s; |
| 312 |
|
} |
| 313 |
|
|
| 314 |
|
/** |
| 319 |
|
inline void mul( Real s, const RectMatrix<Real, Row, Col>& m ) { |
| 320 |
|
for (unsigned int i = 0; i < Row; i++) |
| 321 |
|
for (unsigned int j = 0; j < Col; j++) |
| 322 |
< |
data_[i][j] = s * m.data_[i][j]; |
| 322 |
> |
this->data_[i][j] = s * m.data_[i][j]; |
| 323 |
|
} |
| 324 |
|
|
| 325 |
|
/** |
| 329 |
|
inline void div( Real s) { |
| 330 |
|
for (unsigned int i = 0; i < Row; i++) |
| 331 |
|
for (unsigned int j = 0; j < Col; j++) |
| 332 |
< |
data_[i][j] /= s; |
| 332 |
> |
this->data_[i][j] /= s; |
| 333 |
|
} |
| 334 |
|
|
| 335 |
|
/** |
| 340 |
|
inline void div( Real s, const RectMatrix<Real, Row, Col>& m ) { |
| 341 |
|
for (unsigned int i = 0; i < Row; i++) |
| 342 |
|
for (unsigned int j = 0; j < Col; j++) |
| 343 |
< |
data_[i][j] = m.data_[i][j] / s; |
| 343 |
> |
this->data_[i][j] = m.data_[i][j] / s; |
| 344 |
|
} |
| 345 |
|
|
| 346 |
|
/** |
| 385 |
|
|
| 386 |
|
for (unsigned int i = 0; i < Row; i++) |
| 387 |
|
for (unsigned int j = 0; j < Col; j++) |
| 388 |
< |
result(j, i) = data_[i][j]; |
| 388 |
> |
result(j, i) = this->data_[i][j]; |
| 389 |
|
|
| 390 |
|
return result; |
| 391 |
|
} |