Lapacke в C: использование dsymm, ошибка компоновщика - PullRequest
Новая
       74

Lapacke в C: использование dsymm, ошибка компоновщика

0 голосов
/ 05 июня 2018

** РЕДАКТИРОВАТЬ, включая обновление в конце вопроса.

Когда я пытаюсь использовать dsymm из LAPACKE, я получаю сообщение об ошибке «Ошибка команды компоновщика».(Итак, возможно, я неправильно компилирую код?) Вот код, о котором идет речь: 

#include "matrix_multiplication_attempt.h"
#include <stdio.h>
#include "lapacke.h"
int main ( )
{

/* 3x3 A matrix (symmetric) */
double a[] = {1,2,6,
              2,3,1,
              6,1,4};
/* 3x3 B Matrix */
double b[] = {2,3,4,
              3,6,7,
              4,7,4};
/* 3x3 C Matrix */
double c[] = {0,0,0,
              0,0,0,
              0,0,0};


char side, uplo;
int M,N, lda, ldb, ldc;
double alpha, beta,info;

side= 'L';
uplo='L';
M=3;
N=3;
alpha=1.0;
beta=0.0;
lda=3;
ldb=3;
info=8.0;

info=cblas_dsymm(side, uplo,
              M,N, alpha, a,lda,
              b,ldb,beta,c,ldc);
return info;

И вот ошибка, которую я получаю:

matrix_multiplication_attempt.c:51:10: warning: implicit declaration of function 'lapacke_dsymm' is invalid in C99 [-Wimplicit-function-declaration] info=lapacke_dsymm(side, uplo, ^ 1 warning generated. Undefined symbols for architecture x86_64: "_lapacke_dsymm", referenced from: _main in matrix_multiplication_attempt-e2c0b9.o ld: symbol(s) not found for architecture x86_64 clang: error: linker command failed with exit code 1 (use -v to see invocation)

У меня нет опыта работы с C, поэтому, возможно, есть простая ошибка.Любое понимание очень ценится.

РЕДАКТИРОВАТЬ: Спасибо, Фрэнсис, проблема была в том, что dsymm не является частью LAPACK, но в BLAS и ваше решение работает.Теперь мой обновленный вопрос: как я могу связать LAPACKE, чтобы я мог использовать как LAPACKE, так и BLAS подпрограммы в одном файле?Одна из моих попыток компиляции:

gcc matrix_multiplication_attempt.c -o matrix_multiplication_attempt -lblas -Wall -I/usr/local/opt/lapack/include -L/usr/local/opt/lapack/lib/ -llapacke

, но это дает ошибку Undefined symbols for architecture x86_64: "_cblas_dsymm", referenced from: _main in matrix_multiplication_attempt-76b8f6.o ld: symbol(s) not found for architecture x86_64 clang: error: linker command failed with exit code 1 (use -v to see invocation)

1 Ответ

0 голосов
/ 05 июня 2018

В отличие от dsyev, подпрограмма dsymm отсутствует в библиотеке Лапака.Действительно, это реализовано в библиотеке Blas, в BLAS / SRC / dsymm.f .Как следствие, в LAPACKE нет интерфейса с этой подпрограммой.

Лучшая альтернатива для программ на C - это ссылка на CBLAS и использование cblas_dsymm().

Здесьпример кода: 

#include<stdio.h>
#include<math.h>
#include<complex.h>
#include "cblas.h"

int main ( int argc, char **argv ){
    int i,j;

    int m=2;
    int n=3;
    double alpha=1.0;
    double beta=1.0;
    double AA[2][2]={{2,0},{0,1}};
    double BB[2][3]={{1,0,0},{1,0,0}};
    double CC[2][3]={{1,0,0},{1,0,0}};
    int lda=2;
    int ldb=n;
    int ldc=n;
    cblas_dsymm(CblasRowMajor,CblasLeft ,CblasUpper,m,n,alpha, &AA[0][0], lda, &BB[0][0], ldb, beta, &CC[0][0], ldc);

    for(i=0;i<m;++i)
    {
        for(j=0;j<n;++j)
            printf("%lf " ,CC[i][j]);
        printf("\n");
    }
    return 0;
}

Я удивлен, увидев, что я успешно скомпилировал его: gcc main.c -o main -lblas -Wall.Действительно, я думал, что необходимо связать cblas с помощью gcc main.c -o main -lcblas -lblas -Wall, но не в моем случае.

...