gsl_matrix 的初始化

Question

我正在使用 GNU 科学库，我想用值初始化一个矩阵，但我不明白如果没有循环该怎么做:

这个有效：

  gsl_matrix * m = gsl_matrix_alloc (3, 3);

  for (int i = 0; i < 3; i++)
    for (int j = 0; j < 3; j++)
      gsl_matrix_set (m, i, j, 100*i + j);

但我想做这样的事情:

    double data[] = { i  , 1/2.0, 1/3.0,
                      1/2.0, 1/3.0, 1/4.0,
                      1/3.0, 1/4.0, 1/5.0};

    gsl_matrix mat = gsl_matrix_from_array(data); // does not exists

有办法吗？

Answer 1

您可以使用 std::copy() 来执行此操作。

#include <iostream>
#include <gsl/gsl_matrix.h>

using namespace std;

int main(void){

    double data[] = {1,2,3,4};
    gsl_matrix *m = gsl_matrix_alloc(2,2);
    copy(data, data+4, m->data);

    printf("Matrix:\n%f %f\n%f %f\n", gsl_matrix_get(m, 0, 0), gsl_matrix_get(m, 0, 1),
        gsl_matrix_get(m, 1, 0), gsl_matrix_get(m, 1, 1));
    return 0;
}

gsl_matrix 的初始化

initialization of a gsl_matrix

c++

gsl