gsl c ++ в проблеме cygwin - PullRequest
Новая
       22

gsl c ++ в проблеме cygwin

0 голосов
/ 25 ноября 2010

Я использую C ++ для решения методом граничных элементов, у меня проблема с моим кодом: 

#include <stdio.h> 
#include <math.h> 
#include <gsl/gsl_linalg.h> 
using namespace std; 
string STRING; 
int i,q,s;  
int const nodes = 16; 
double A[nodes][nodes];  
double b[nodes]; 
 int main(){ 
  for (i=0;i<nodes;i++) 
       { 


           { 
            A[q][i] = 1.;       
            b[q] = 1.; 
           }        
       }

Как только A и b построены, нам нужно решить систему Ax = b с помощью
вычисление обратного значения x = A ^ (- 1) * b 

gsl_matrix_view m = gsl_matrix_view_array (*A, nodes, nodes); for(q=0;q<nodes;q++)
gsl_matrix_view b = gsl_matrix_view_array (b, nodes, nodes); 
gsl_vector *x = gsl_vector_alloc (nodes); 
gsl_permutation * p = gsl_permutation_alloc (nodes); 
gsl_linalg_LU_decomp (&m.matrix, p, &s); 
gsl_linalg_LU_solve (&m.matrix, p, &b.vector, x); 

      return 0;    
     }

Когда я компилирую в Cygwin, используя 

g++ test.cpp -lm -lgsl -o bem.out -L/usr/bin

, я получаю следующую ошибку: 

test.cpp: In function 'int main()':
test.cpp:39:59: error: cannot convert 'gsl_matrix_view' to 'double*' for argument '1' to '_gsl_matrix_view gsl_matrix_view_array(double*, size_t, size_t)'
test.cpp:43:39: error: 'struct gsl_matrix_view' has no member named 'vector'

Я следовал тому же примеру, что и учебник по GSL, но я получаю эти ошибки.Кто-нибудь может помочь?Я был бы очень признателен.

Спасибо!

1 Ответ

0 голосов
/ 17 февраля 2011

Вы повторно объявили b как gsl_matrix_view в той же строке, на которую жалуется компиляция.Вы объявили его как массив double ранее.

...