поставил библиотеку clapack для облегчения работы с программированием мат. моделей. Помогите найти (или составить) пример решения СЛАУ или вычисления какого-либо параметра матрицы, например, границ спектра. (Бьюсь уже не одни сутки.)
Код:
# pwd
/opt/CLAPACK-3.1.1.1
# ls
BLAS COPYING INSTALL lapackUbuntu.a README.install tmglibUbuntu.a
blasUbuntu.a F2CLIBS LAPACK-revisions3.1.0.info Makefile SRC
CLAPACK-revisions3.1.1.1.info INCLUDE LAPACK-revisions3.1.1.info make.inc TESTING
Добавлено спустя 33 минуты 6 секунд:Re: найти пример на clapackВообще устроит пример для любой библиотеки на языке C++, с вызовом высокоуровневых функций обработки матриц: решение СЛАУ, вычисление числа обусловленности или собственных чисел, и т.д. Нашел их много lapackcpp, lpp, atlas-dev-cpp, jama ... но запустить нифига не получается ... . Исключение boost::uBLAS но в нем не такие общие операции .
В общем, если есть опыт поделитесь.
Добавлено спустя 1 час 26 минут 58 секунд:появился свой вариант
Код:
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <math.h>
#include <iostream>
main(){
boost::numeric::ublas::matrix<double> A(2,2);
A(0,0)=1;A(0,1)=2;
A(1,0)=3;A(1,1)=2;
// for octave
std::cout << "A=[";
for(size_t i=0; i<2; ++i){
for(size_t j=0; j<2; ++j){
std::cout << A(i,j) << ' ';
}
if(i<1)std::cout << ";";
}
std::cout << "]; disp(eig(A))";
}
Код:
$ c++ test.cc -o test && ./test | octave -q
-1
4