2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 найти пример на clapack
Сообщение27.12.2008, 14:04 
поставил библиотеку 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

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group