2011-10-13 2 views
0

gsl_matrix_set을 16 번 작성하는 대신 gsl_matrix 데이터 포인터에 값을 할당하기위한 바로 가기를 갖고 싶지만 작동시키지 못합니다. 작은 자체에 포함 된 샘플 :gsl_matrix 구조체 포인터에 -> 사용하기

#include <stdio.h> 
#include <stdlib.h> 
#include <gsl/gsl_linalg.h> 
#include <gsl/gsl_matrix.h> 

gsl_matrix* get_rates(double k) { 

    int nChar = 4; 
    gsl_matrix *rates = gsl_matrix_calloc(nChar, nChar); 

//TODO: this apparently fails 
    rates->data = (double[]) {-1, k/(k+2), 1/(k+2), 1/(k+2), 
     k/(k+2), -1, 1/(k+2), 1/(k+2), 
     1/(k+2), 1/(k+2), -1, k/(k+2), 
     1/(k+2), 1/(k+2), k/(k+2), -1}; 

    return (rates); 

} 

int main(void) { 

    gsl_matrix *rates; 

    rates = get_rates(2.0); 

    gsl_matrix_fprintf (stdout, rates, " %f"); 

     return(0); 
} 

반환 :

-1.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 

답변