On 21 Sty, 09:28, Wit Jakuczun <w...@g...com> wrote:
> On 21 Sty, 09:18, bartekltg <b...@g...com> wrote:
>
> > > Ale mogil sie postarac a wrecz powinni. Wystarczy rozpoznac
> > > ze mnozymy wektor przez macierz i dac goto do wyspecjalizowanej
> > > procedurki. Tak samo jesli mnozymy macierz [NxK] przez [KxM] dla
> > > M i N znacznie wiekszych niz K, to jakas wyspecjalizowana procedura
> > > wymnozylaby duzo szybciej niz ogolna.
>
> > Dla bardzo malych K mozna by sie bawic w rozwijanie petli.
> > Dla wiekszych - nie ma znaczenia przy algorytmie n^3.
>
> Może to rzuci jakieś światło na całą dyskusję (slaj
16ty):http://www-users.mat.uni.torun.pl/~bala/wyklad
_arch/wshop2000.pdf
> Złożoność ciągle jest N^3 ale działa szybciej niż pętelki.
>
> Blokowe mnożenie wykorzystuje ATLAS. Idea jest taka, żeby
> cache maskymalnie efektywnie wykorzystać.
> Ja kiedyś robiłem porównania ATLAS vs ręczny kod i ATLAS był
> dużo szybszy. Z tego co pamiętam to oprócz mnożenia rozwiązywałem
> układy równań.
>
Uzupełnienie:

http://www.mimuw.edu.pl/~przykry/matmult.php

Pozdrawiam,
Wit Jakuczun