Mercurial > hg > segmenter-vamp-plugin

annotate armadillo-2.4.4/include/armadillo_bits/eop_core_meat.hpp @ 0:8b6102e2a9b0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Armadillo Library
author maxzanoni76 <max.zanoni@eecs.qmul.ac.uk>
date Wed, 11 Apr 2012 09:27:06 +0100
parents
children
rev   line source
max@0 1 // Copyright (C) 2010-2011 NICTA (www.nicta.com.au)
max@0 2 // Copyright (C) 2010-2011 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 eop_core
max@0 15 //! @{
max@0 16
max@0 17
max@0 18 #undef arma_applier_1
max@0 19 #undef arma_applier_2
max@0 20 #undef arma_applier_3
max@0 21 #undef operatorA
max@0 22
max@0 23 #define arma_applier_1(operatorA) \
max@0 24 {\
max@0 25 uword i,j;\
max@0 26 \
max@0 27 for(i=0, j=1; j<n_elem; i+=2, j+=2)\
max@0 28 {\
max@0 29 eT tmp_i = P[i];\
max@0 30 eT tmp_j = P[j];\
max@0 31 \
max@0 32 tmp_i = eop_core<eop_type>::process(tmp_i, k);\
max@0 33 tmp_j = eop_core<eop_type>::process(tmp_j, k);\
max@0 34 \
max@0 35 out_mem[i] operatorA tmp_i;\
max@0 36 out_mem[j] operatorA tmp_j;\
max@0 37 }\
max@0 38 \
max@0 39 if(i < n_elem)\
max@0 40 {\
max@0 41 out_mem[i] operatorA eop_core<eop_type>::process(P[i], k);\
max@0 42 }\
max@0 43 }
max@0 44
max@0 45
max@0 46 #define arma_applier_2(operatorA) \
max@0 47 {\
max@0 48 uword count = 0;\
max@0 49 \
max@0 50 for(uword col=0; col<n_cols; ++col)\
max@0 51 {\
max@0 52 uword i,j;\
max@0 53 \
max@0 54 for(i=0, j=1; j<n_rows; i+=2, j+=2, count+=2)\
max@0 55 {\
max@0 56 eT tmp_i = P.at(i,col);\
max@0 57 eT tmp_j = P.at(j,col);\
max@0 58 \
max@0 59 tmp_i = eop_core<eop_type>::process(tmp_i, k);\
max@0 60 tmp_j = eop_core<eop_type>::process(tmp_j, k);\
max@0 61 \
max@0 62 out_mem[count ] operatorA tmp_i;\
max@0 63 out_mem[count+1] operatorA tmp_j;\
max@0 64 }\
max@0 65 \
max@0 66 if(i < n_rows)\
max@0 67 {\
max@0 68 out_mem[count] operatorA eop_core<eop_type>::process(P.at(i,col), k);\
max@0 69 ++count;\
max@0 70 }\
max@0 71 }\
max@0 72 }
max@0 73
max@0 74
max@0 75
max@0 76 #define arma_applier_3(operatorA) \
max@0 77 {\
max@0 78 uword count = 0;\
max@0 79 \
max@0 80 for(uword slice=0; slice<n_slices; ++slice)\
max@0 81 {\
max@0 82 for(uword col=0; col<n_cols; ++col)\
max@0 83 {\
max@0 84 uword i,j;\
max@0 85 \
max@0 86 for(i=0, j=1; j<n_rows; i+=2, j+=2, count+=2)\
max@0 87 {\
max@0 88 eT tmp_i = P.at(i,col,slice);\
max@0 89 eT tmp_j = P.at(j,col,slice);\
max@0 90 \
max@0 91 tmp_i = eop_core<eop_type>::process(tmp_i, k);\
max@0 92 tmp_j = eop_core<eop_type>::process(tmp_j, k);\
max@0 93 \
max@0 94 out_mem[count ] operatorA tmp_i;\
max@0 95 out_mem[count+1] operatorA tmp_j;\
max@0 96 }\
max@0 97 \
max@0 98 if(i < n_rows)\
max@0 99 {\
max@0 100 out_mem[count] operatorA eop_core<eop_type>::process(P.at(i,col,slice), k);\
max@0 101 ++count;\
max@0 102 }\
max@0 103 }\
max@0 104 }\
max@0 105 }
max@0 106
max@0 107
max@0 108
max@0 109 //
max@0 110 // matrices
max@0 111
max@0 112
max@0 113
max@0 114 template<typename eop_type>
max@0 115 template<typename T1>
max@0 116 arma_hot
max@0 117 inline
max@0 118 void
max@0 119 eop_core<eop_type>::apply(Mat<typename T1::elem_type>& out, const eOp<T1, eop_type>& x)
max@0 120 {
max@0 121 arma_extra_debug_sigprint();
max@0 122
max@0 123 typedef typename T1::elem_type eT;
max@0 124
max@0 125 const uword n_rows = out.n_rows;
max@0 126 const uword n_cols = out.n_cols;
max@0 127 const uword n_elem = out.n_elem;
max@0 128
max@0 129 // NOTE: we're assuming that the matrix has already been set to the correct size and there is no aliasing;
max@0 130 // size setting and alias checking is done by either the Mat contructor or operator=()
max@0 131
max@0 132 const eT k = x.aux;
max@0 133 eT* out_mem = out.memptr();
max@0 134
max@0 135 if(Proxy<T1>::prefer_at_accessor == false)
max@0 136 {
max@0 137 typename Proxy<T1>::ea_type P = x.P.get_ea();
max@0 138
max@0 139 arma_applier_1(=);
max@0 140 }
max@0 141 else
max@0 142 {
max@0 143 const Proxy<T1>& P = x.P;
max@0 144
max@0 145 arma_applier_2(=);
max@0 146 }
max@0 147 }
max@0 148
max@0 149
max@0 150
max@0 151 template<typename eop_type>
max@0 152 template<typename T1>
max@0 153 arma_hot
max@0 154 inline
max@0 155 void
max@0 156 eop_core<eop_type>::apply_inplace_plus(Mat<typename T1::elem_type>& out, const eOp<T1, eop_type>& x)
max@0 157 {
max@0 158 arma_extra_debug_sigprint();
max@0 159
max@0 160 typedef typename T1::elem_type eT;
max@0 161
max@0 162 const uword n_rows = x.get_n_rows();
max@0 163 const uword n_cols = x.get_n_cols();
max@0 164
max@0 165 arma_debug_assert_same_size(out.n_rows, out.n_cols, n_rows, n_cols, "addition");
max@0 166
max@0 167 eT* out_mem = out.memptr();
max@0 168 const uword n_elem = out.n_elem;
max@0 169
max@0 170 const eT k = x.aux;
max@0 171
max@0 172 if(Proxy<T1>::prefer_at_accessor == false)
max@0 173 {
max@0 174 typename Proxy<T1>::ea_type P = x.P.get_ea();
max@0 175
max@0 176 arma_applier_1(+=);
max@0 177 }
max@0 178 else
max@0 179 {
max@0 180 const Proxy<T1>& P = x.P;
max@0 181
max@0 182 arma_applier_2(+=);
max@0 183 }
max@0 184 }
max@0 185
max@0 186
max@0 187
max@0 188 template<typename eop_type>
max@0 189 template<typename T1>
max@0 190 arma_hot
max@0 191 inline
max@0 192 void
max@0 193 eop_core<eop_type>::apply_inplace_minus(Mat<typename T1::elem_type>& out, const eOp<T1, eop_type>& x)
max@0 194 {
max@0 195 arma_extra_debug_sigprint();
max@0 196
max@0 197 typedef typename T1::elem_type eT;
max@0 198
max@0 199 const uword n_rows = x.get_n_rows();
max@0 200 const uword n_cols = x.get_n_cols();
max@0 201
max@0 202 arma_debug_assert_same_size(out.n_rows, out.n_cols, n_rows, n_cols, "subtraction");
max@0 203
max@0 204 eT* out_mem = out.memptr();
max@0 205 const uword n_elem = out.n_elem;
max@0 206
max@0 207 const eT k = x.aux;
max@0 208
max@0 209 if(Proxy<T1>::prefer_at_accessor == false)
max@0 210 {
max@0 211 typename Proxy<T1>::ea_type P = x.P.get_ea();
max@0 212
max@0 213 arma_applier_1(-=);
max@0 214 }
max@0 215 else
max@0 216 {
max@0 217 const Proxy<T1>& P = x.P;
max@0 218
max@0 219 arma_applier_2(-=);
max@0 220 }
max@0 221 }
max@0 222
max@0 223
max@0 224
max@0 225 template<typename eop_type>
max@0 226 template<typename T1>
max@0 227 arma_hot
max@0 228 inline
max@0 229 void
max@0 230 eop_core<eop_type>::apply_inplace_schur(Mat<typename T1::elem_type>& out, const eOp<T1, eop_type>& x)
max@0 231 {
max@0 232 arma_extra_debug_sigprint();
max@0 233
max@0 234 typedef typename T1::elem_type eT;
max@0 235
max@0 236 const uword n_rows = x.get_n_rows();
max@0 237 const uword n_cols = x.get_n_cols();
max@0 238
max@0 239 arma_debug_assert_same_size(out.n_rows, out.n_cols, n_rows, n_cols, "element-wise multiplication");
max@0 240
max@0 241 eT* out_mem = out.memptr();
max@0 242 const uword n_elem = out.n_elem;
max@0 243
max@0 244 const eT k = x.aux;
max@0 245
max@0 246 if(Proxy<T1>::prefer_at_accessor == false)
max@0 247 {
max@0 248 typename Proxy<T1>::ea_type P = x.P.get_ea();
max@0 249
max@0 250 arma_applier_1(*=);
max@0 251 }
max@0 252 else
max@0 253 {
max@0 254 const Proxy<T1>& P = x.P;
max@0 255
max@0 256 arma_applier_2(*=);
max@0 257 }
max@0 258 }
max@0 259
max@0 260
max@0 261
max@0 262 template<typename eop_type>
max@0 263 template<typename T1>
max@0 264 arma_hot
max@0 265 inline
max@0 266 void
max@0 267 eop_core<eop_type>::apply_inplace_div(Mat<typename T1::elem_type>& out, const eOp<T1, eop_type>& x)
max@0 268 {
max@0 269 arma_extra_debug_sigprint();
max@0 270
max@0 271 typedef typename T1::elem_type eT;
max@0 272
max@0 273 const uword n_rows = x.get_n_rows();
max@0 274 const uword n_cols = x.get_n_cols();
max@0 275
max@0 276 arma_debug_assert_same_size(out.n_rows, out.n_cols, n_rows, n_cols, "element-wise division");
max@0 277
max@0 278 eT* out_mem = out.memptr();
max@0 279 const uword n_elem = out.n_elem;
max@0 280
max@0 281 const eT k = x.aux;
max@0 282
max@0 283 if(Proxy<T1>::prefer_at_accessor == false)
max@0 284 {
max@0 285 typename Proxy<T1>::ea_type P = x.P.get_ea();
max@0 286
max@0 287 arma_applier_1(/=);
max@0 288 }
max@0 289 else
max@0 290 {
max@0 291 const Proxy<T1>& P = x.P;
max@0 292
max@0 293 arma_applier_2(/=);
max@0 294 }
max@0 295 }
max@0 296
max@0 297
max@0 298
max@0 299 //
max@0 300 // cubes
max@0 301
max@0 302
max@0 303
max@0 304 template<typename eop_type>
max@0 305 template<typename T1>
max@0 306 arma_hot
max@0 307 inline
max@0 308 void
max@0 309 eop_core<eop_type>::apply(Cube<typename T1::elem_type>& out, const eOpCube<T1, eop_type>& x)
max@0 310 {
max@0 311 arma_extra_debug_sigprint();
max@0 312
max@0 313 typedef typename T1::elem_type eT;
max@0 314
max@0 315 const uword n_rows = out.n_rows;
max@0 316 const uword n_cols = out.n_cols;
max@0 317 const uword n_slices = out.n_slices;
max@0 318 const uword n_elem = out.n_elem;
max@0 319
max@0 320 // NOTE: we're assuming that the matrix has already been set to the correct size and there is no aliasing;
max@0 321 // size setting and alias checking is done by either the Mat contructor or operator=()
max@0 322
max@0 323 const eT k = x.aux;
max@0 324 eT* out_mem = out.memptr();
max@0 325
max@0 326 if(ProxyCube<T1>::prefer_at_accessor == false)
max@0 327 {
max@0 328 typename ProxyCube<T1>::ea_type P = x.P.get_ea();
max@0 329
max@0 330 arma_applier_1(=);
max@0 331 }
max@0 332 else
max@0 333 {
max@0 334 const ProxyCube<T1>& P = x.P;
max@0 335
max@0 336 arma_applier_3(=);
max@0 337 }
max@0 338 }
max@0 339
max@0 340
max@0 341
max@0 342 template<typename eop_type>
max@0 343 template<typename T1>
max@0 344 arma_hot
max@0 345 inline
max@0 346 void
max@0 347 eop_core<eop_type>::apply_inplace_plus(Cube<typename T1::elem_type>& out, const eOpCube<T1, eop_type>& x)
max@0 348 {
max@0 349 arma_extra_debug_sigprint();
max@0 350
max@0 351 typedef typename T1::elem_type eT;
max@0 352
max@0 353 const uword n_rows = x.get_n_rows();
max@0 354 const uword n_cols = x.get_n_cols();
max@0 355 const uword n_slices = x.get_n_slices();
max@0 356
max@0 357 arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "addition");
max@0 358
max@0 359 eT* out_mem = out.memptr();
max@0 360 const uword n_elem = out.n_elem;
max@0 361
max@0 362 const eT k = x.aux;
max@0 363
max@0 364 if(ProxyCube<T1>::prefer_at_accessor == false)
max@0 365 {
max@0 366 typename ProxyCube<T1>::ea_type P = x.P.get_ea();
max@0 367
max@0 368 arma_applier_1(+=);
max@0 369 }
max@0 370 else
max@0 371 {
max@0 372 const ProxyCube<T1>& P = x.P;
max@0 373
max@0 374 arma_applier_3(+=);
max@0 375 }
max@0 376 }
max@0 377
max@0 378
max@0 379
max@0 380 template<typename eop_type>
max@0 381 template<typename T1>
max@0 382 arma_hot
max@0 383 inline
max@0 384 void
max@0 385 eop_core<eop_type>::apply_inplace_minus(Cube<typename T1::elem_type>& out, const eOpCube<T1, eop_type>& x)
max@0 386 {
max@0 387 arma_extra_debug_sigprint();
max@0 388
max@0 389 typedef typename T1::elem_type eT;
max@0 390
max@0 391 const uword n_rows = x.get_n_rows();
max@0 392 const uword n_cols = x.get_n_cols();
max@0 393 const uword n_slices = x.get_n_slices();
max@0 394
max@0 395 arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "subtraction");
max@0 396
max@0 397 eT* out_mem = out.memptr();
max@0 398 const uword n_elem = out.n_elem;
max@0 399
max@0 400 const eT k = x.aux;
max@0 401
max@0 402 if(ProxyCube<T1>::prefer_at_accessor == false)
max@0 403 {
max@0 404 typename ProxyCube<T1>::ea_type P = x.P.get_ea();
max@0 405
max@0 406 arma_applier_1(-=);
max@0 407 }
max@0 408 else
max@0 409 {
max@0 410 const ProxyCube<T1>& P = x.P;
max@0 411
max@0 412 arma_applier_3(-=);
max@0 413 }
max@0 414 }
max@0 415
max@0 416
max@0 417
max@0 418 template<typename eop_type>
max@0 419 template<typename T1>
max@0 420 arma_hot
max@0 421 inline
max@0 422 void
max@0 423 eop_core<eop_type>::apply_inplace_schur(Cube<typename T1::elem_type>& out, const eOpCube<T1, eop_type>& x)
max@0 424 {
max@0 425 arma_extra_debug_sigprint();
max@0 426
max@0 427 typedef typename T1::elem_type eT;
max@0 428
max@0 429 const uword n_rows = x.get_n_rows();
max@0 430 const uword n_cols = x.get_n_cols();
max@0 431 const uword n_slices = x.get_n_slices();
max@0 432
max@0 433 arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "element-wise multiplication");
max@0 434
max@0 435 eT* out_mem = out.memptr();
max@0 436 const uword n_elem = out.n_elem;
max@0 437
max@0 438 const eT k = x.aux;
max@0 439
max@0 440 if(ProxyCube<T1>::prefer_at_accessor == false)
max@0 441 {
max@0 442 typename ProxyCube<T1>::ea_type P = x.P.get_ea();
max@0 443
max@0 444 arma_applier_1(*=);
max@0 445 }
max@0 446 else
max@0 447 {
max@0 448 const ProxyCube<T1>& P = x.P;
max@0 449
max@0 450 arma_applier_3(*=);
max@0 451 }
max@0 452 }
max@0 453
max@0 454
max@0 455
max@0 456 template<typename eop_type>
max@0 457 template<typename T1>
max@0 458 arma_hot
max@0 459 inline
max@0 460 void
max@0 461 eop_core<eop_type>::apply_inplace_div(Cube<typename T1::elem_type>& out, const eOpCube<T1, eop_type>& x)
max@0 462 {
max@0 463 arma_extra_debug_sigprint();
max@0 464
max@0 465 typedef typename T1::elem_type eT;
max@0 466
max@0 467 const uword n_rows = x.get_n_rows();
max@0 468 const uword n_cols = x.get_n_cols();
max@0 469 const uword n_slices = x.get_n_slices();
max@0 470
max@0 471 arma_debug_assert_same_size(out.n_rows, out.n_cols, out.n_slices, n_rows, n_cols, n_slices, "element-wise division");
max@0 472
max@0 473 eT* out_mem = out.memptr();
max@0 474 const uword n_elem = out.n_elem;
max@0 475
max@0 476 const eT k = x.aux;
max@0 477
max@0 478 if(ProxyCube<T1>::prefer_at_accessor == false)
max@0 479 {
max@0 480 typename ProxyCube<T1>::ea_type P = x.P.get_ea();
max@0 481
max@0 482 arma_applier_1(/=);
max@0 483 }
max@0 484 else
max@0 485 {
max@0 486 const ProxyCube<T1>& P = x.P;
max@0 487
max@0 488 arma_applier_3(/=);
max@0 489 }
max@0 490 }
max@0 491
max@0 492
max@0 493
max@0 494 //
max@0 495 // common
max@0 496
max@0 497
max@0 498
max@0 499 template<typename eop_type>
max@0 500 template<typename eT>
max@0 501 arma_hot
max@0 502 arma_pure
max@0 503 arma_inline
max@0 504 eT
max@0 505 eop_core<eop_type>::process(const eT val, const eT k)
max@0 506 {
max@0 507 arma_ignore(val);
max@0 508 arma_ignore(k);
max@0 509
max@0 510 arma_stop("eop_core::process(): unhandled eop_type");
max@0 511 return eT(0);
max@0 512 }
max@0 513
max@0 514
max@0 515
max@0 516 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 517 eop_core<eop_scalar_plus >::process(const eT val, const eT k) { return val + k; }
max@0 518
max@0 519 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 520 eop_core<eop_scalar_minus_pre >::process(const eT val, const eT k) { return k - val; }
max@0 521
max@0 522 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 523 eop_core<eop_scalar_minus_post>::process(const eT val, const eT k) { return val - k; }
max@0 524
max@0 525 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 526 eop_core<eop_scalar_times >::process(const eT val, const eT k) { return val * k; }
max@0 527
max@0 528 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 529 eop_core<eop_scalar_div_pre >::process(const eT val, const eT k) { return k / val; }
max@0 530
max@0 531 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 532 eop_core<eop_scalar_div_post >::process(const eT val, const eT k) { return val / k; }
max@0 533
max@0 534 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 535 eop_core<eop_square >::process(const eT val, const eT ) { return val*val; }
max@0 536
max@0 537 template<> template<typename eT> arma_hot arma_const arma_inline eT
max@0 538 eop_core<eop_neg >::process(const eT val, const eT ) { return eop_aux::neg(val); }
max@0 539
max@0 540 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 541 eop_core<eop_sqrt >::process(const eT val, const eT ) { return eop_aux::sqrt(val); }
max@0 542
max@0 543 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 544 eop_core<eop_log >::process(const eT val, const eT ) { return eop_aux::log(val); }
max@0 545
max@0 546 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 547 eop_core<eop_log2 >::process(const eT val, const eT ) { return eop_aux::log2(val); }
max@0 548
max@0 549 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 550 eop_core<eop_log10 >::process(const eT val, const eT ) { return eop_aux::log10(val); }
max@0 551
max@0 552 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 553 eop_core<eop_trunc_log >::process(const eT val, const eT ) { return arma::trunc_log(val); }
max@0 554
max@0 555 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 556 eop_core<eop_exp >::process(const eT val, const eT ) { return eop_aux::exp(val); }
max@0 557
max@0 558 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 559 eop_core<eop_exp2 >::process(const eT val, const eT ) { return eop_aux::exp2(val); }
max@0 560
max@0 561 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 562 eop_core<eop_exp10 >::process(const eT val, const eT ) { return eop_aux::exp10(val); }
max@0 563
max@0 564 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 565 eop_core<eop_trunc_exp >::process(const eT val, const eT ) { return arma::trunc_exp(val); }
max@0 566
max@0 567 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 568 eop_core<eop_cos >::process(const eT val, const eT ) { return eop_aux::cos(val); }
max@0 569
max@0 570 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 571 eop_core<eop_sin >::process(const eT val, const eT ) { return eop_aux::sin(val); }
max@0 572
max@0 573 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 574 eop_core<eop_tan >::process(const eT val, const eT ) { return eop_aux::tan(val); }
max@0 575
max@0 576 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 577 eop_core<eop_acos >::process(const eT val, const eT ) { return eop_aux::acos(val); }
max@0 578
max@0 579 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 580 eop_core<eop_asin >::process(const eT val, const eT ) { return eop_aux::asin(val); }
max@0 581
max@0 582 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 583 eop_core<eop_atan >::process(const eT val, const eT ) { return eop_aux::atan(val); }
max@0 584
max@0 585 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 586 eop_core<eop_cosh >::process(const eT val, const eT ) { return eop_aux::cosh(val); }
max@0 587
max@0 588 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 589 eop_core<eop_sinh >::process(const eT val, const eT ) { return eop_aux::sinh(val); }
max@0 590
max@0 591 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 592 eop_core<eop_tanh >::process(const eT val, const eT ) { return eop_aux::tanh(val); }
max@0 593
max@0 594 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 595 eop_core<eop_acosh >::process(const eT val, const eT ) { return eop_aux::acosh(val); }
max@0 596
max@0 597 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 598 eop_core<eop_asinh >::process(const eT val, const eT ) { return eop_aux::asinh(val); }
max@0 599
max@0 600 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 601 eop_core<eop_atanh >::process(const eT val, const eT ) { return eop_aux::atanh(val); }
max@0 602
max@0 603 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 604 eop_core<eop_eps >::process(const eT val, const eT ) { return eop_aux::direct_eps(val); }
max@0 605
max@0 606 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 607 eop_core<eop_abs >::process(const eT val, const eT ) { return eop_aux::arma_abs(val); }
max@0 608
max@0 609 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 610 eop_core<eop_conj >::process(const eT val, const eT ) { return eop_aux::conj(val); }
max@0 611
max@0 612 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 613 eop_core<eop_pow >::process(const eT val, const eT k) { return eop_aux::pow(val, k); }
max@0 614
max@0 615 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 616 eop_core<eop_floor >::process(const eT val, const eT ) { return eop_aux::floor(val); }
max@0 617
max@0 618 template<> template<typename eT> arma_hot arma_pure arma_inline eT
max@0 619 eop_core<eop_ceil >::process(const eT val, const eT ) { return eop_aux::ceil(val); }
max@0 620
max@0 621
max@0 622
max@0 623 #undef arma_applier_1
max@0 624 #undef arma_applier_2
max@0 625 #undef arma_applier_3
max@0 626
max@0 627
max@0 628
max@0 629 //! @}