[LUGOS] Control-Alt-delete

[Metod Kozelj]

Howdy!

igor wrote:

>>/Kaj pa tisti choleskyjev razcep, ki ga ponuja BLAS? Streljam na slepo, ker to ni moje področje, ampak se mi zdi, da če ti algoritmi že obstajajo, verjetno obstajajo tudi optimizirane verzije. Vem, da obstajata knjižnici LAPACK in BLAS za linearno algebro tudi na Linuxu, vendar nisem še nikdar delal z njima. /
>>
> Kolikor sem spegal, sta BLAS in LAPACK napisana v navadnem C-ju.  
> Nekje sem nasel celo
> funkcije napisane v fortranu.  Te pa ne morejo biti hitrejse od 
> tistih, ki jih znam sam napisati v
> C-ju  (Choleskyijev razcep se napise v dveh zankah oziroma desetih 
> vrsticah).  Mislim, da so
> oni morali napisati algoritem prav v asemblerju.  Bi pa bil prav 
> vesel, ce bi me kdo opozoril,
> ce spet streljam mimo.


Eh no, stvar lahko napišeš tako ali drugače. Verjamem, da jaz npr. 
Fourierovo transformacijo napišem drugače (beri: bistveno manj 
optimalno) kot pa je izvedena v FFT proceduri od Numerical Recipes (in 
C). Morda bi pa Cholesky-jev razcep, napisan v 25 vrsticah C-ja, v 
resnici deloval 300% hitreje?

-- 
Peace!
  Mkx

---- perl -e 'print $i=pack(c5,(41*2),sqrt(7056),(unpack(c,H)-2),oct(115),10);'