получить размер gsl_vector_view - PullRequest
Новая
       22

получить размер gsl_vector_view

0 голосов
/ 05 января 2019

Как мне получить размер gsl_vector_view? У меня есть следующая функция для выборки многомерного гауссовского вектора 

void mvnrnd(gsl_vector *x, gsl_matrix *Sigma,gsl_vector *Mu, int K, const gsl_rng *seed){

    gsl_matrix *A = gsl_matrix_alloc(K, K);
    gsl_matrix_memcpy(A, Sigma);
    gsl_linalg_cholesky_decomp (A);

    for (int k=0; k<K; k++){
        gsl_vector_set (x, k, gsl_ran_ugaussian (seed));
        }

    gsl_blas_dtrmv (CblasLower, CblasNoTrans, CblasNonUnit, A, x);
    gsl_vector_add (x, Mu);

    }

Чтобы передать векторы и матрицы, я делаю следующее: 

gsl_matrix *MuH    = gsl_matrix_calloc(Kest*Kest,1);

gsl_matrix *vecH   = gsl_matrix_calloc(Kest*Kest,1);
gsl_matrix *SigmaH = gsl_matrix_calloc(Kest*Kest,Kest*Kest);
gsl_matrix_memcpy (SigmaH, &Q_view.matrix);
gsl_vector_view vecH_view = gsl_matrix_subcolumn (vecH, 0, 0, Kest*Kest);
gsl_vector_view MuH_view =  gsl_matrix_column (MuB, 0);
mvnrnd(&vecH_view.vector, SigmaH , &MuH_view.vector, Kest*Kest, seed);

но я получаю эту ошибку: 

gsl: ./oper_source.c:27: ERROR: vectors must have same length
Default GSL error handler invoked.

Итак, я хотел бы напечатать gsl_vector_view. Как это сделать? Я использовал & vecH_view.vector->size, и он возвращает ошибку.

1 Ответ

0 голосов
/ 07 января 2019

«Метод проб и ошибок» всегда работает для меня. Ниже приведен минимальный пример: 

#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_matrix.h>
int main()
{
    gsl_matrix * y = gsl_matrix_calloc(6,6);
    gsl_vector_view z = gsl_matrix_column(y,0);
    printf("size = %ld\n",(&z.vector)->size);
    gsl_matrix_free(y);
}

$ gcc -Wall -pedantic-errors -O0 se.c -lgsl -lgslcblas -lm && ./a.out
size = 6
...