W dniu 2012-04-01 19:37, M.M. pisze:

> W oryginalnym kodzie mam coś innego. Mam dane zero-jedynkowe i buduję z nich
> układ równań normalnych. Dane mam upakowane. Upakowanie polega
> na tym że zamiast zer i jedynek trzymam numery pozycji na których
> są jedynki. Dane są trudne, mają znamiona funkcji losowej, dlatego

Jedna ze standardowych metod zapisu macierzy rzadkich,
dość rozsądnie.

> w przykładzie dałem rand. Mniej/więcej coś takiego:

Ale te indeksy jedynek trzymasz chyba posortowane
po długiej wspołrzednej dla każej 'krótkiej' osobno?

> for( int i=0 ; i<N ; i++ )
> for( int j=i ; j<N ; j++ )
> macierz[ jedynki[i] * rozmiar_wiersza + jedynki[j] ] ++;


'Macierz' to już ta X^t X?
Oznaczenia
http://en.wikipedia.org/wiki/Numerical_methods_for_l
inear_least_squares#The_general_problem


A skompresowane za pomocą indeksów jest samo X?

Coś tu jest nie tak;)

pzdr
bartekltg