Компиляция кода C ++ с помощью библиотеки числовых привязок Boost для решения линейной системы Ax = b
-
11-09-2019 - |
Вопрос
Я использую Привязки числовой библиотеки для Boost uBLAS чтобы решить простую линейную систему:
#include<boost/numeric/ublas/matrix.hpp>
#include<boost/numeric/ublas/io.hpp>
#include<boost/numeric/bindings/traits/ublas_matrix.hpp>
#include<boost/numeric/bindings/lapack/gesv.hpp>
#include <boost/numeric/bindings/traits/ublas_vector2.hpp>
namespace ublas = boost::numeric::ublas;
namespace lapack= boost::numeric::bindings::lapack;
int main()
{
ublas::matrix<float,ublas::column_major> A(3,3);
ublas::vector<float> b(3);
for(unsigned i=0;i < A.size1();i++)
for(unsigned j =0;j < A.size2();j++)
{
std::cout << "enter element "<<i << j << std::endl;
std::cin >> A(i,j);
}
std::cout << A << std::endl;
b(0) = 21; b(1) = 1; b(2) = 17;
lapack::gesv(A,b);
std::cout << b << std::endl;
return 0;
}
Я попытался скомпилировать его с помощью следующей команды:
g++ -I/home/foolb/.boost/include/boost-1_38 -I/home/foolb/.boostnumbind/include/boost-numeric-bindings solve_Axb_byhand.cc -o solve_Axb_byhand
но сбой со следующей ошибкой:
/media/disk/tmp/ccbd973l.o: In function `boost::numeric::bindings::lapack::detail::gesv(int, int, float*, int, int*, float*, int, int*)':
solve_Axb_byhand2.cc:(.text._ZN5boost7numeric8bindings6lapack6detail4gesvEiiPfiPiS4_iS5_[boost::numeric::bindings::lapack::detail::gesv(int, int, float*, int, int*, float*, int, int*)]+0x59): undefined reference to `sgesv_'
collect2: ld returned 1 exit status
Что не так с моим подходом в коде?
Решение
sgesv_ - это символ из библиотеки LAPACK, вам нужно будет сослаться на него.уБЛАС просто привязывается к этому, я думаю.
Хотя я тоже не знаю названия библиотеки :)
Другие советы
Извините, если это не соответствует действительности, но я не вижу, чтобы вы ссылались на библиотеки boost в вашей команде g ++.Я вижу, что вы включаете пути поиска, но нет явного включения самих скомпилированных библиотек Boost;что-то вроде -lboost (боюсь, я не знаю точного формата, который вам понадобится, и это вполне может зависеть от местоположения).
При связывании с библиотекой числовой привязки boost вы можете связать с параметром
-Lpath/to/lapack -llapack -Lpath/to/blas -lblas -lgfortran
в gcc4
-Lpath/to/lapack -llapack -Lpath/to/blas -lblas -lg2c
в gcc3