Вызов DPOTRS из LAPACK в C на GNU / Linux - PullRequest
Новая
       8

Вызов DPOTRS из LAPACK в C на GNU / Linux

0 голосов
/ 26 февраля 2012

Я бы хотел обратиться за помощью к функции DPOTRS из LAPACK при вызове изнутри C. Да, я знаю, что матрица, с которой я пытаюсь работать, является положительно определенной (с собственными значениями 3, 1 на самом деле!)

Сейчас моя функция работает неправильно.Возвращает неверный результат. 

double A[] = {2.0, 1.0, 1.0, 2.0};
double b[] = {1.5, 0.0};
printf("%5.3f %5.3f\n", b[0], b[1]);
info = dpotrs('U',2,1,A,2,b,2);
printf("%d\n", info);
printf("%5.3f %5.3f\n", b[0], b[1]);

Это код.dpotrs импортируется вручную с использованием этого кода: 

static long dpotrs(char UPLO, long N, long NRHS, double* A, long LDA, double* B, long LDB)
{
    extern void dpotrs_(char* UPLOp, long* Np, long* NRHSp, double* A, long* LDAp, double* B, long* LDBp, long* infop);
    long info;
    dpotrs_(&UPLO, &N, &NRHS, A, &LDA, B, &LDB, &info);
    return info;
}

Возвращает неверный результат!Правильный результат будет 1.000, -0.500

Тем не менее, я получаю 0.469, -0.188

Но вдобавок ко всему, информация возвращает 0, как будто все прошло хорошо!* И это, дамы и господа, уму непостижимо.

Заранее спасибо!

1 Ответ

1 голос
/ 26 февраля 2012

Я думаю, вам нужно прочитать документацию: 

*
*  Purpose
*  =======
*
*  DPOTRS solves a system of linear equations A*X = B with a symmetric
*  positive definite matrix A using the Cholesky factorization
*  A = U**T*U or A = L*L**T computed by DPOTRF.
*

Матрица ввода для DPOTRS должна быть уже факторизованной матрицей Холецкого в верхнем или нижнем треугольном формате LAPACK. Вы пропустили шаг. Факторизовать ваш A с помощью DPOTRF сначала , а затем использовать его для решения одной или нескольких правых сторон с помощью DPOTRS.

...