monolish  0.14.0
MONOlithic LIner equation Solvers for Highly-parallel architecture
blas/vector/vecsub.hpp
Go to the documentation of this file.
1 #pragma once
2 
3 namespace monolish {
4 namespace {
5 template <typename F1, typename F2, typename F3>
6 void vecsub_core(const F1 &a, const F2 &b, F3 &y) {
7  Logger &logger = Logger::get_instance();
8  logger.func_in(monolish_func);
9 
10  const size_t aoffset = a.get_offset();
11  const size_t boffset = b.get_offset();
12  const size_t yoffset = y.get_offset();
13 
14  // err
15  assert(util::is_same_size(a, b, y));
16  assert(util::is_same_device_mem_stat(a, b, y));
17 
18  internal::vsub(y.size(), a.data() + aoffset, b.data() + boffset,
19  y.data() + yoffset, y.get_device_mem_stat());
20 
21  logger.func_out();
22 }
23 } // namespace
24 } // namespace monolish
monolish::util::is_same_size
bool is_same_size(const T &x, const U &y)
compare size of vector or 1Dview (same as is_same_structure())
Definition: monolish_common.hpp:358
monolish_func
#define monolish_func
Definition: monolish_logger.hpp:9
monolish::util::is_same_device_mem_stat
bool is_same_device_mem_stat(const T &arg1, const U &arg2)
compare same device memory status
Definition: monolish_common.hpp:431
monolish
Definition: monolish_matrix_blas.hpp:9
monolish::Logger::get_instance
static Logger & get_instance()
Definition: monolish_logger.hpp:42