Как выполнить матрицу обратной операции, используя ускоренную структуру?

Я хотел бы найти обратную матрицу.

Я знаю, что это включает в себя сначала факторизацию LU, а затем этап инверсии, но я не могу найти требуемую функцию путем поиска в документах Apple по 10,7!

Это похоже на полезный постИнверсия симметричной матрицы в C с использованием CBLAS / LAPACK, указывая на то, чтоsgetrf_ а такжеsgetri_ функции должны быть использованы. Однако в поиске этих терминов я ничего не нахожу в документах Xcode.

У кого-нибудь есть код котельной плиты для этой матричной операции?

Ответы на вопрос(1)

Решение Вопроса

Я полагаю, что Apple не документирует код LAPACK, потому что они просто реализуют стандартный интерфейсnetlib.org, Обидно, что вы не можете искать имена этих функций из встроенных документов Xcode, однако решение довольно простое: просто укажите имя функции в URL, например. заdgetrf_() идти к,http://www.netlib.org/clapack/what/double/dgetrf.c.

Чтобы инвертировать матрицу, нужны две функции LAPACK:dgetrf_(), который выполняет факторизацию LU, иdgetri_() который принимает выходные данные предыдущей функции и выполняет фактическую инверсию.

Я создал стандартный проект приложения с использованием XCode, добавил Accelerate Framework, создал два файла C: matinv.h, matinv.c и отредактировал файл main.m, чтобы удалить объекты Cocoa:

// main.m

#import "matinv.h"

int main(int argc, char *argv[])
{
    int N = 3;
    double A[N*N];
    A[0] = 1; A[1] = 1; A[2] = 7;
    A[3] = 1; A[4] = 2; A[5] = 1;
    A[6] = 1; A[7] = 1; A[8] = 3;
    matrix_invert(N, A);
    //        [ -1.25  -1.0  3.25 ]
    // A^-1 = [  0.5       1.0  -1.5  ]
    //        [  0.25   0.0 -0.25 ] 
    return 0;
}

Теперь заголовочный файл,

//  matinv.h

int matrix_invert(int N, double *matrix);

а затем исходный файл,

int matrix_invert(int N, double *matrix) {

    int error=0;
    int *pivot = malloc(N*sizeof(int)); // LAPACK requires MIN(M,N), here M==N, so N will do fine.
    double *workspace = malloc(N*sizeof(double));

    /*  LU factorisation */
    dgetrf_(&N, &N, matrix, &N, pivot, &error);

    if (error != 0) {
        NSLog(@"Error 1");
        free(pivot);
        free(workspace);
        return error;
    }

    /*  matrix inversion */
    dgetri_(&N, matrix, &N, pivot, workspace, &N, &error);

    if (error != 0) {
        NSLog(@"Error 2");
        free(pivot);
        free(workspace);
        return error;
    }

    free(pivot);
    free(workspace);
    return error;
}
 Daniel Farrell22 июн. 2013 г., 04:33
Благодарю. Исправлено выше.
 12 июн. 2013 г., 10:32
О, я нашел это: sgetrf_ и sgetri_ (S для "одинарной точности"?)
 04 июл. 2012 г., 04:02
Каноническая ссылка на LAPACK - это руководство пользователя LAPACK. (netlib.org/lapack/lug)
 12 июн. 2013 г., 08:47
У меня возникают проблемы при сканировании этой загадочной (по меньшей мере) библиотеки, которая является LAPACK. Как я могу адаптировать этот код к плавающей запятой одинарной точности?
 21 июн. 2013 г., 21:50
Спасибо за ваш код, это помогло мне. В вашем matrix_invert есть ошибка, если есть ошибка, вы пропустили выпуск pivot и рабочей области и вернулись раньше времени.

Ваш ответ на вопрос