Re: FORTRAN && C

["Nefedov U." <nefedov@nusun.jinr.ru>]

On Wed, 30 Oct 2002, Alexey Ozeritsky wrote:

> 
> проблема в том, что после этого
> ничего не работает.
> например функция для решения ур-я AX=B,
> с помощью LU разложения, если я не ошибаюсь:
> dgesv(N, NRHS, A, LDA, IPIV, B, LDB, INFO )
> 
> в прграмме пишу
> double *A,*X;
> int *C;
> int n;
> int INFO;
> .......
> n=5;
> ......
> A=malloc(n*n*sizeof(double));
> X=malloc(n*n*sizeof(double));
> C=malloc(n*n*sizeof(double));
> .......
> dgesv_(n,n,A,n,C,X,n,INFO);
> .......
> (в этом случае должно происходить обращение матрицы)
> 
> и получаю segmentation fault
> #1  0x080493c6 in main (argc=0, argv=0x0) at main.c:77
> 77              dgesv_(n,1,A,n,C,X,n,INFO);
> 
> 

  В фортране все аргументы передаются по адресу.
  Правильно писать

 dgesv_(&n,&n,A,&n,C,X,&n,&INFO);
 
  А вообще то я пользуюсь "обертками": либо lapack++, 
  либо mtl. 

  lapack++ - маленькая библиотека и будет полезна
  чтобы смотреть как там осуществляется вызов lapack-овских
  программок. Ответ я подсмотрел именно в ней.

  Успехов. 
  Юра.