-
Path: news-archive.icm.edu.pl!news.gazeta.pl!not-for-mail
From: " M.M." <m...@N...gazeta.pl>
Newsgroups: pl.comp.programming
Subject: Re: dalsza optymalizacja
Date: Sun, 1 Apr 2012 20:29:53 +0000 (UTC)
Organization: "Portal Gazeta.pl -> http://www.gazeta.pl"
Lines: 53
Message-ID: <jlads1$h1i$1@inews.gazeta.pl>
References: <jl73ie$b6f$1@inews.gazeta.pl> <jl9aaf$326$1@inews.gazeta.pl>
<jl9al3$ee7$1@inews.gazeta.pl> <jl9dfl$a80$1@inews.gazeta.pl>
<jl9e7f$epj$1@node2.news.atman.pl> <jl9ftm$j00$1@inews.gazeta.pl>
<jl9h5r$i0o$1@node2.news.atman.pl> <jl9ju9$o28$1@inews.gazeta.pl>
<jl9nc4$oha$1@node2.news.atman.pl> <jl9pbd$s10$1@inews.gazeta.pl>
<jl9pjl$qts$1@node2.news.atman.pl> <jl9udd$jka$1@inews.gazeta.pl>
<jla0s1$37m$1@node2.news.atman.pl> <jla3o5$hv7$1@inews.gazeta.pl>
<jla5pu$8kp$1@node2.news.atman.pl> <jla87q$2on$1@inews.gazeta.pl>
<jlabvl$fau$1@node2.news.atman.pl>
NNTP-Posting-Host: localhost
Content-Type: text/plain; charset=ISO-8859-2
Content-Transfer-Encoding: 8bit
X-Trace: inews.gazeta.pl 1333312193 17458 172.20.26.236 (1 Apr 2012 20:29:53 GMT)
X-Complaints-To: u...@a...pl
NNTP-Posting-Date: Sun, 1 Apr 2012 20:29:53 +0000 (UTC)
X-User: mariotti
X-Forwarded-For: 89.229.34.123
X-Remote-IP: localhost
Xref: news-archive.icm.edu.pl pl.comp.programming:196521
[ ukryj nagłówki ]bartekltg <b...@g...com> napisał(a):
> TrochÄ skromnie to opisaĹeĹ i nie do koca widzÄ, jak to robisz.
>
> Jak dokĹadnie zapisujesz X i ma m*n. jedynki[] to tablica
> z ktĂłrymi pozycjami?
Pisząc z pamięci:
Dla jednej pary x i y:
x to wektor zer i jedynek.
y to wartość skalarna którą chcemy aproksymować.
m to macierz układu równań normalnych
x ma rozmiar x.size // notacja za Cormenem
m ma rozmiar x.size wierszy i (x.size+1) kolumn, ostatnia kolumna to wyrazy
wolne równania.
Na początku m jest wyzerowane.
Rozszerzam x przez przylaczenie na koniec y:
x[x.size] = y;
x.size = x.size+1
potem w dwóch pętlach:
for( i=0 ; i<x.size-1 ; i++ )
for( j=0 ; j<x.size ; j++ )
m[i][j] += x[i] * x[j];
I tak dla każdego wektora. Potem oddzielna sprawa rozwiązać ten układ
równań.
Pierwsza optymalizacja:
Ze względu na to że macierz jest symetryczna, to można wyliczyć tylko
jeden trójkąt:
for( i=0 ; i<x.size-1 ; i++ )
for( j=i ; j<x.size ; j++ )
m[i][j] += x[i] * x[j];
I druga optymalizacja:
ze względu na to że dane to zera i jedynki, zapamiętuję w
jedynki[] pozycje jedynek (posortowane) i wychodzi:
for( i=0 ; i<jedynki.size ; i++ ) {
for( j=i ; j<jedynki.size ; j++ )
m[ jedynki[i] ][ jedynki[j] ] ++ ; // x[ jedynki[i] ] * x[ jedynki[j] ];
m[i][x.size-1] += y;
}
Można coś ulepszyć?
Pozdrawiam
--
Wysłano z serwisu Usenet w portalu Gazeta.pl -> http://www.gazeta.pl/usenet/
Następne wpisy z tego wątku
- 01.04.12 22:44 bartekltg
- 01.04.12 22:49 bartekltg
- 01.04.12 22:50
- 01.04.12 23:56 M.M.
- 01.04.12 23:59 bartekltg
- 02.04.12 00:11 M.M.
- 02.04.12 00:13 M.M.
- 02.04.12 01:20 bartekltg
- 02.04.12 02:58 M.M.
- 02.04.12 08:25
- 02.04.12 10:40 zażółcony
- 02.04.12 10:46
- 02.04.12 11:46
- 02.04.12 14:58 bartekltg
- 02.04.12 15:00 bartekltg
Najnowsze wątki z tej grupy
- Alg. kompresji LZW
- Popr. 14. Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- Arch. Prog. Nieuprzywilejowanych w pełnej wer. na nowej s. WWW energokod.pl
- 7. Raport Totaliztyczny: Sprawa Qt Group wer. 424
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
Najnowsze wątki
- 2025-02-17 Kraków => MS Dynamics 365BC/NAV Developer <=
- 2025-02-17 Chrzanów => Programista NodeJS <=
- 2025-02-17 Warszawa => Node.js / Fullstack Developer <=
- 2025-02-17 Białystok => System Architect (Java background) <=
- 2025-02-17 Białystok => Solution Architect (Java background) <=
- 2025-02-17 Gliwice => Team Lead / Tribe Lead FrontEnd <=
- 2025-02-17 Gdańsk => PHP Developer <=
- 2025-02-17 Warszawa => Senior ASP.NET Developer <=
- 2025-02-17 Gliwice => Business Development Manager - Network and Network Security
- 2025-02-17 Mińsk Mazowiecki => Area Sales Manager OZE <=
- 2025-02-17 Odśnieżanie samochodu
- 2025-02-17 Katowice => Regionalny Kierownik Sprzedaży (OZE) <=
- 2025-02-17 Dęblin => JavaScript / Node / Fullstack Developer <=
- 2025-02-17 Pompiarze...
- 2025-02-16 PV teraz