Armadillo не выполняет оптимизацию подкачки - PullRequest
Новая
       12

Armadillo не выполняет оптимизацию подкачки

0 голосов
/ 11 ноября 2018

Как поменять объекты броненосца, например, arma::vec с, без обмена содержимым? 

void f5()
{
  arma::vec x(10);
  arma::vec y(10);
  std::cout << &x[2] << ", " << &y[2] << "\n";
  x.swap(y);
  std::cout << &x[2] << ", " << &y[2];
}

Вышеприведенный код выводит 

0x24fbe50, 0x24fbef0
0x24fbe50, 0x24fbef0

Спасибо!

Ответы [ 2 ]

0 голосов
/ 11 ноября 2018

Небольшое расширение ответа @ OZ17. Armadillo, кажется, хранит данные с размерами <16 локально <code>mem_local и более крупными в области, обозначенной mem

From GDB:
> p x
{
  <arma::Mat<double>> = {
    <arma::Base<double, arma::Mat<double> >> = {
    <arma::Base_inv_yes<arma::Mat<double> >> = {<No data fields>}, 
    <arma::Base_eval_Mat<double, arma::Mat<double> >> = {<No data fields>}, 
    <arma::Base_trans_default<arma::Mat<double> >> = {<No data fields>}, <No data fields>}, 
    members of arma::Mat<double>: 
    n_rows = 1, 
    n_cols = 10, 
    n_elem = 10, 
    vec_state = 2, 
    mem_state = 0, 
    mem = 0x7fffffffb830, 
    mem_local = {0 <repeats 16 times>},   
    static is_col = false, 
    static is_row = false
  }, 
  members of arma::Row<double>: 
  static is_col = false, 
  static is_row = false
}

и небольшой пример для его визуализации: 

arma::rowvec x(10,arma::fill::ones);
arma::rowvec y(10,arma::fill::zeros);
std::cout << "Size=10" << std::endl;
std::cout << "&x=" << x.memptr() << ", x[0..4]=" << x.subvec(1,5);
std::cout << "&y=" << y.memptr() << ", y[0..4]=" << y.subvec(1,5);
x.swap(y);
std::cout << "x.swap(y)" << std::endl;
std::cout << "&x=" << x.memptr() << ", x[0..4]=" << x.subvec(1,5);
std::cout << "&y=" << y.memptr() << ", y[0..4]=" << y.subvec(1,5);

arma::rowvec x2(17,arma::fill::ones);
arma::rowvec y2(17,arma::fill::zeros);
std::cout << "\nSize=17" << std::endl;
std::cout << "&x=" << x2.memptr() << ", x[0..4]=" << x2.subvec(1,5);
std::cout << "&y=" << y2.memptr() << ", y[0..4]=" << y2.subvec(1,5);
x2.swap(y2);
std::cout << "x.swap(y)" << std::endl;
std::cout << "&x=" << x2.memptr() << ", x[0..4]=" << x2.subvec(1,5);
std::cout << "&y=" << y2.memptr() << ", y[0..4]=" << y2.subvec(1,5);

Вывод из примера показывает, что в обоих случаях содержимое поменяно местами, но для небольших массивов оно поменял местную область памяти, а для более крупного случая поменял местами указатель mem. 

Size=10
&x=0x7fffffffb830, x[0..4]=   1.0000   1.0000   1.0000   1.0000   1.0000
&y=0x7fffffffb8e0, y[0..4]=        0        0        0        0        0
x.swap(y)
&x=0x7fffffffb830, x[0..4]=        0        0        0        0        0
&y=0x7fffffffb8e0, y[0..4]=   1.0000   1.0000   1.0000   1.0000   1.0000

Size=17
&x=0x5555557d7fd0, x[0..4]=   1.0000   1.0000   1.0000   1.0000   1.0000
&y=0x5555557d8060, y[0..4]=        0        0        0        0        0
x.swap(y)
&x=0x5555557d8060, x[0..4]=        0        0        0        0        0
&y=0x5555557d7fd0, y[0..4]=   1.0000   1.0000   1.0000   1.0000   1.0000
0 голосов
/ 11 ноября 2018

Похоже, что своп броненосца внутренне является memcpy ниже определенного размера массива (согласно op <= 16). </p>

...