Witam serdecznie,
Oczywiscie temat jest znany, ale nie znalazlem na forum lub w Internecie dokladnie tego co by mi bylo potrzebne. Otoz mam dosc spore macierze do odwracania za pomoca programu napisanego w Fortranie. Oczywiscie "spore" jest wzgledne, w moim przypadku sa do rzedy od kilku do kilku tysiecy. Sa to macierze, ktore pojawiaja sie w analizie metoda najmniejszych kwadratow. W tym celu zaimplementowalem algorytm QR odwracajacy macierz (znalazlem go na strone http://beige.ucs.indiana.edu/B673/node32.html). Jest to w moim przypadku istotne bo jednoczenie chce znajdowac wartosci i wektory wlasne i ew. je filtrowac. Jednakze, co mnie troche zmartwilo to algorytm ten daje gorsze wyniki niz testy zrobione za pomoca najprostrzego algorytmu Gaussa. Gorsze w sensie, ze sa zauwazalne roznice miedzy obiema metodami, przy czym z uwagi na ostateczne wyniki wiem, ze metoda Gaussa dziala dobrze. Podejrzewam, ze metoda oparta na faktoryzacji QR nie jest stabila numerycznie... Z gory bede wdzieczny za sugestie dotyczace jaka inna metode ew. uzyc albo na co zwracac uwage przy implementacji. Oczywiscie moge w dowolnym momencie napisac wiecej i ew. przeslac dane liczbowe.
Pozdrawiam,
Radek