Mercurial > hg > segmenter-vamp-plugin

annotate armadillo-2.4.4/include/armadillo_bits/op_relational_meat.hpp @ 5:79b343f3e4b8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
In thi version the problem of letters assigned to each segment has been solved.
author maxzanoni76 <max.zanoni@eecs.qmul.ac.uk>
date Wed, 11 Apr 2012 13:48:13 +0100
parents 8b6102e2a9b0
children
rev   line source
max@0 1 // Copyright (C) 2009-2012 NICTA (www.nicta.com.au)
max@0 2 // Copyright (C) 2009-2012 Conrad Sanderson
max@0 3 //
max@0 4 // This file is part of the Armadillo C++ library.
max@0 5 // It is provided without any warranty of fitness
max@0 6 // for any purpose. You can redistribute this file
max@0 7 // and/or modify it under the terms of the GNU
max@0 8 // Lesser General Public License (LGPL) as published
max@0 9 // by the Free Software Foundation, either version 3
max@0 10 // of the License or (at your option) any later version.
max@0 11 // (see http://www.opensource.org/licenses for more info)
max@0 12
max@0 13
max@0 14 //! \addtogroup op_relational
max@0 15 //! @{
max@0 16
max@0 17
max@0 18 #undef operator_rel
max@0 19
max@0 20 #undef arma_applier_mat_pre
max@0 21 #undef arma_applier_mat_post
max@0 22
max@0 23 #undef arma_applier_cube_pre
max@0 24 #undef arma_applier_cube_post
max@0 25
max@0 26
max@0 27 #define arma_applier_mat_pre(operator_rel) \
max@0 28 {\
max@0 29 typedef typename T1::elem_type eT;\
max@0 30 typedef typename Proxy<T1>::ea_type ea_type;\
max@0 31 \
max@0 32 const eT val = X.aux;\
max@0 33 \
max@0 34 const Proxy<T1> P(X.m);\
max@0 35 \
max@0 36 const uword n_rows = P.get_n_rows();\
max@0 37 const uword n_cols = P.get_n_cols();\
max@0 38 \
max@0 39 const bool bad_alias = ( Proxy<T1>::has_subview && P.is_alias(out) );\
max@0 40 \
max@0 41 if(bad_alias == false)\
max@0 42 {\
max@0 43 out.set_size(n_rows, n_cols);\
max@0 44 \
max@0 45 uword* out_mem = out.memptr();\
max@0 46 \
max@0 47 if(Proxy<T1>::prefer_at_accessor == false)\
max@0 48 {\
max@0 49 ea_type PA = P.get_ea();\
max@0 50 const uword n_elem = out.n_elem;\
max@0 51 \
max@0 52 for(uword i=0; i<n_elem; ++i)\
max@0 53 {\
max@0 54 out_mem[i] = (val operator_rel PA[i]) ? uword(1) : uword(0);\
max@0 55 }\
max@0 56 }\
max@0 57 else\
max@0 58 {\
max@0 59 uword count = 0;\
max@0 60 \
max@0 61 for(uword col=0; col < n_cols; ++col)\
max@0 62 for(uword row=0; row < n_rows; ++row, ++count)\
max@0 63 {\
max@0 64 out_mem[count] = (val operator_rel P.at(row,col)) ? uword(1) : uword(0);\
max@0 65 }\
max@0 66 }\
max@0 67 }\
max@0 68 else\
max@0 69 {\
max@0 70 const unwrap<typename Proxy<T1>::stored_type> tmp(P.Q);\
max@0 71 \
max@0 72 out = (val) operator_rel (tmp.M);\
max@0 73 }\
max@0 74 }
max@0 75
max@0 76
max@0 77
max@0 78 #define arma_applier_mat_post(operator_rel) \
max@0 79 {\
max@0 80 typedef typename T1::elem_type eT;\
max@0 81 typedef typename Proxy<T1>::ea_type ea_type;\
max@0 82 \
max@0 83 const eT val = X.aux;\
max@0 84 \
max@0 85 const Proxy<T1> P(X.m);\
max@0 86 \
max@0 87 const uword n_rows = P.get_n_rows();\
max@0 88 const uword n_cols = P.get_n_cols();\
max@0 89 \
max@0 90 const bool bad_alias = ( Proxy<T1>::has_subview && P.is_alias(out) );\
max@0 91 \
max@0 92 if(bad_alias == false)\
max@0 93 {\
max@0 94 out.set_size(n_rows, n_cols);\
max@0 95 \
max@0 96 uword* out_mem = out.memptr();\
max@0 97 \
max@0 98 if(Proxy<T1>::prefer_at_accessor == false)\
max@0 99 {\
max@0 100 ea_type PA = P.get_ea();\
max@0 101 const uword n_elem = out.n_elem;\
max@0 102 \
max@0 103 for(uword i=0; i<n_elem; ++i)\
max@0 104 {\
max@0 105 out_mem[i] = (PA[i] operator_rel val) ? uword(1) : uword(0);\
max@0 106 }\
max@0 107 }\
max@0 108 else\
max@0 109 {\
max@0 110 uword count = 0;\
max@0 111 \
max@0 112 for(uword col=0; col < n_cols; ++col)\
max@0 113 for(uword row=0; row < n_rows; ++row, ++count)\
max@0 114 {\
max@0 115 out_mem[count] = (P.at(row,col) operator_rel val) ? uword(1) : uword(0);\
max@0 116 }\
max@0 117 }\
max@0 118 }\
max@0 119 else\
max@0 120 {\
max@0 121 const unwrap<typename Proxy<T1>::stored_type> tmp(P.Q);\
max@0 122 \
max@0 123 out = (tmp.M) operator_rel (val);\
max@0 124 }\
max@0 125 }
max@0 126
max@0 127
max@0 128
max@0 129 #define arma_applier_cube_pre(operator_rel) \
max@0 130 {\
max@0 131 typedef typename T1::elem_type eT;\
max@0 132 typedef typename ProxyCube<T1>::ea_type ea_type;\
max@0 133 \
max@0 134 const eT val = X.aux;\
max@0 135 \
max@0 136 const ProxyCube<T1> P(X.m);\
max@0 137 \
max@0 138 const uword n_rows = P.get_n_rows();\
max@0 139 const uword n_cols = P.get_n_cols();\
max@0 140 const uword n_slices = P.get_n_slices();\
max@0 141 \
max@0 142 const bool bad_alias = ( ProxyCube<T1>::has_subview && P.is_alias(out) );\
max@0 143 \
max@0 144 if(bad_alias == false)\
max@0 145 {\
max@0 146 out.set_size(n_rows, n_cols, n_slices);\
max@0 147 \
max@0 148 uword* out_mem = out.memptr();\
max@0 149 \
max@0 150 if(ProxyCube<T1>::prefer_at_accessor == false)\
max@0 151 {\
max@0 152 ea_type PA = P.get_ea();\
max@0 153 const uword n_elem = out.n_elem;\
max@0 154 \
max@0 155 for(uword i=0; i<n_elem; ++i)\
max@0 156 {\
max@0 157 out_mem[i] = (val operator_rel PA[i]) ? uword(1) : uword(0);\
max@0 158 }\
max@0 159 }\
max@0 160 else\
max@0 161 {\
max@0 162 uword count = 0;\
max@0 163 \
max@0 164 for(uword slice=0; slice < n_slices; ++slice)\
max@0 165 for(uword col=0; col < n_cols; ++col)\
max@0 166 for(uword row=0; row < n_rows; ++row, ++count)\
max@0 167 {\
max@0 168 out_mem[count] = (val operator_rel P.at(row,col,slice)) ? uword(1) : uword(0);\
max@0 169 }\
max@0 170 }\
max@0 171 }\
max@0 172 else\
max@0 173 {\
max@0 174 const unwrap_cube<typename ProxyCube<T1>::stored_type> tmp(P.Q);\
max@0 175 \
max@0 176 out = (val) operator_rel (tmp.M);\
max@0 177 }\
max@0 178 }
max@0 179
max@0 180
max@0 181
max@0 182 #define arma_applier_cube_post(operator_rel) \
max@0 183 {\
max@0 184 typedef typename T1::elem_type eT;\
max@0 185 typedef typename ProxyCube<T1>::ea_type ea_type;\
max@0 186 \
max@0 187 const eT val = X.aux;\
max@0 188 \
max@0 189 const ProxyCube<T1> P(X.m);\
max@0 190 \
max@0 191 const uword n_rows = P.get_n_rows();\
max@0 192 const uword n_cols = P.get_n_cols();\
max@0 193 const uword n_slices = P.get_n_slices();\
max@0 194 \
max@0 195 const bool bad_alias = ( ProxyCube<T1>::has_subview && P.is_alias(out) );\
max@0 196 \
max@0 197 if(bad_alias == false)\
max@0 198 {\
max@0 199 out.set_size(n_rows, n_cols, n_slices);\
max@0 200 \
max@0 201 uword* out_mem = out.memptr();\
max@0 202 \
max@0 203 if(ProxyCube<T1>::prefer_at_accessor == false)\
max@0 204 {\
max@0 205 ea_type PA = P.get_ea();\
max@0 206 const uword n_elem = out.n_elem;\
max@0 207 \
max@0 208 for(uword i=0; i<n_elem; ++i)\
max@0 209 {\
max@0 210 out_mem[i] = (PA[i] operator_rel val) ? uword(1) : uword(0);\
max@0 211 }\
max@0 212 }\
max@0 213 else\
max@0 214 {\
max@0 215 uword count = 0;\
max@0 216 \
max@0 217 for(uword slice=0; slice < n_slices; ++slice)\
max@0 218 for(uword col=0; col < n_cols; ++col)\
max@0 219 for(uword row=0; row < n_rows; ++row, ++count)\
max@0 220 {\
max@0 221 out_mem[count] = (P.at(row,col,slice) operator_rel val) ? uword(1) : uword(0);\
max@0 222 }\
max@0 223 }\
max@0 224 }\
max@0 225 else\
max@0 226 {\
max@0 227 const unwrap_cube<typename ProxyCube<T1>::stored_type> tmp(P.Q);\
max@0 228 \
max@0 229 out = (tmp.M) operator_rel (val);\
max@0 230 }\
max@0 231 }
max@0 232
max@0 233
max@0 234
max@0 235 template<typename T1>
max@0 236 inline
max@0 237 void
max@0 238 op_rel_lt_pre::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_lt_pre>& X)
max@0 239 {
max@0 240 arma_extra_debug_sigprint();
max@0 241
max@0 242 arma_applier_mat_pre( < );
max@0 243 }
max@0 244
max@0 245
max@0 246
max@0 247 template<typename T1>
max@0 248 inline
max@0 249 void
max@0 250 op_rel_gt_pre::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_gt_pre>& X)
max@0 251 {
max@0 252 arma_extra_debug_sigprint();
max@0 253
max@0 254 arma_applier_mat_pre( > );
max@0 255 }
max@0 256
max@0 257
max@0 258
max@0 259 template<typename T1>
max@0 260 inline
max@0 261 void
max@0 262 op_rel_lteq_pre::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_lteq_pre>& X)
max@0 263 {
max@0 264 arma_extra_debug_sigprint();
max@0 265
max@0 266 arma_applier_mat_pre( <= );
max@0 267 }
max@0 268
max@0 269
max@0 270
max@0 271 template<typename T1>
max@0 272 inline
max@0 273 void
max@0 274 op_rel_gteq_pre::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_gteq_pre>& X)
max@0 275 {
max@0 276 arma_extra_debug_sigprint();
max@0 277
max@0 278 arma_applier_mat_pre( >= );
max@0 279 }
max@0 280
max@0 281
max@0 282
max@0 283 template<typename T1>
max@0 284 inline
max@0 285 void
max@0 286 op_rel_lt_post::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_lt_post>& X)
max@0 287 {
max@0 288 arma_extra_debug_sigprint();
max@0 289
max@0 290 arma_applier_mat_post( < );
max@0 291 }
max@0 292
max@0 293
max@0 294
max@0 295 template<typename T1>
max@0 296 inline
max@0 297 void
max@0 298 op_rel_gt_post::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_gt_post>& X)
max@0 299 {
max@0 300 arma_extra_debug_sigprint();
max@0 301
max@0 302 arma_applier_mat_post( > );
max@0 303 }
max@0 304
max@0 305
max@0 306
max@0 307 template<typename T1>
max@0 308 inline
max@0 309 void
max@0 310 op_rel_lteq_post::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_lteq_post>& X)
max@0 311 {
max@0 312 arma_extra_debug_sigprint();
max@0 313
max@0 314 arma_applier_mat_post( <= );
max@0 315 }
max@0 316
max@0 317
max@0 318
max@0 319 template<typename T1>
max@0 320 inline
max@0 321 void
max@0 322 op_rel_gteq_post::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_gteq_post>& X)
max@0 323 {
max@0 324 arma_extra_debug_sigprint();
max@0 325
max@0 326 arma_applier_mat_post( >= );
max@0 327 }
max@0 328
max@0 329
max@0 330
max@0 331 template<typename T1>
max@0 332 inline
max@0 333 void
max@0 334 op_rel_eq::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_eq>& X)
max@0 335 {
max@0 336 arma_extra_debug_sigprint();
max@0 337
max@0 338 arma_applier_mat_post( == );
max@0 339 }
max@0 340
max@0 341
max@0 342
max@0 343 template<typename T1>
max@0 344 inline
max@0 345 void
max@0 346 op_rel_noteq::apply(Mat<uword>& out, const mtOp<uword, T1, op_rel_noteq>& X)
max@0 347 {
max@0 348 arma_extra_debug_sigprint();
max@0 349
max@0 350 arma_applier_mat_post( != );
max@0 351 }
max@0 352
max@0 353
max@0 354
max@0 355 //
max@0 356 //
max@0 357 //
max@0 358
max@0 359
max@0 360
max@0 361 template<typename T1>
max@0 362 inline
max@0 363 void
max@0 364 op_rel_lt_pre::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_lt_pre>& X)
max@0 365 {
max@0 366 arma_extra_debug_sigprint();
max@0 367
max@0 368 arma_applier_cube_pre( < );
max@0 369 }
max@0 370
max@0 371
max@0 372
max@0 373 template<typename T1>
max@0 374 inline
max@0 375 void
max@0 376 op_rel_gt_pre::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_gt_pre>& X)
max@0 377 {
max@0 378 arma_extra_debug_sigprint();
max@0 379
max@0 380 arma_applier_cube_pre( > );
max@0 381 }
max@0 382
max@0 383
max@0 384
max@0 385 template<typename T1>
max@0 386 inline
max@0 387 void
max@0 388 op_rel_lteq_pre::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_lteq_pre>& X)
max@0 389 {
max@0 390 arma_extra_debug_sigprint();
max@0 391
max@0 392 arma_applier_cube_pre( <= );
max@0 393 }
max@0 394
max@0 395
max@0 396
max@0 397 template<typename T1>
max@0 398 inline
max@0 399 void
max@0 400 op_rel_gteq_pre::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_gteq_pre>& X)
max@0 401 {
max@0 402 arma_extra_debug_sigprint();
max@0 403
max@0 404 arma_applier_cube_pre( >= );
max@0 405 }
max@0 406
max@0 407
max@0 408
max@0 409 template<typename T1>
max@0 410 inline
max@0 411 void
max@0 412 op_rel_lt_post::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_lt_post>& X)
max@0 413 {
max@0 414 arma_extra_debug_sigprint();
max@0 415
max@0 416 arma_applier_cube_post( < );
max@0 417 }
max@0 418
max@0 419
max@0 420
max@0 421 template<typename T1>
max@0 422 inline
max@0 423 void
max@0 424 op_rel_gt_post::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_gt_post>& X)
max@0 425 {
max@0 426 arma_extra_debug_sigprint();
max@0 427
max@0 428 arma_applier_cube_post( > );
max@0 429 }
max@0 430
max@0 431
max@0 432
max@0 433 template<typename T1>
max@0 434 inline
max@0 435 void
max@0 436 op_rel_lteq_post::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_lteq_post>& X)
max@0 437 {
max@0 438 arma_extra_debug_sigprint();
max@0 439
max@0 440 arma_applier_cube_post( <= );
max@0 441 }
max@0 442
max@0 443
max@0 444
max@0 445 template<typename T1>
max@0 446 inline
max@0 447 void
max@0 448 op_rel_gteq_post::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_gteq_post>& X)
max@0 449 {
max@0 450 arma_extra_debug_sigprint();
max@0 451
max@0 452 arma_applier_cube_post( >= );
max@0 453 }
max@0 454
max@0 455
max@0 456
max@0 457 template<typename T1>
max@0 458 inline
max@0 459 void
max@0 460 op_rel_eq::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_eq>& X)
max@0 461 {
max@0 462 arma_extra_debug_sigprint();
max@0 463
max@0 464 arma_applier_cube_post( == );
max@0 465 }
max@0 466
max@0 467
max@0 468
max@0 469 template<typename T1>
max@0 470 inline
max@0 471 void
max@0 472 op_rel_noteq::apply(Cube<uword>& out, const mtOpCube<uword, T1, op_rel_noteq>& X)
max@0 473 {
max@0 474 arma_extra_debug_sigprint();
max@0 475
max@0 476 arma_applier_cube_post( != );
max@0 477 }
max@0 478
max@0 479
max@0 480
max@0 481 #undef arma_applier_mat_pre
max@0 482 #undef arma_applier_mat_post
max@0 483
max@0 484 #undef arma_applier_cube_pre
max@0 485 #undef arma_applier_cube_post
max@0 486
max@0 487
max@0 488
max@0 489 //! @}