-
Data: 2013-03-27 08:31:45
Temat: Re: Nowoczesne procesory - jak to z nimi jest?
Od: "M.M." <m...@g...com> szukaj wiadomości tego autora
[ pokaż wszystkie nagłówki ]W dniu środa, 27 marca 2013 01:38:25 UTC+1 użytkownik bartekltg napisał:
[...]
> A, nie mówiąc już o tym, że sama procedura sumowania może
> doprowadzić do błędów!
> Wynik powinien byc 100 000b, a jest 10 000b
Może rozmawiamy o innym przypadku? Ja mówię mniej/więcej o czymś takim:
#include <cstdlib>
#include <cstdio>
#include <ctime>
typedef double ftyp;
typedef const ftyp cftyp;
#define N (100)
#define M (1024)
#define K (1<<24)
#define L (8)
void init( ftyp v[N] ) {
for( int i=0 ; i<N ; i++ )
v[i] = 0;
for( int i=0 ; i<M ; i++ )
v[rand()%N] += 1.0 / M;
}
void compute( ftyp v[N] ) {
for( int i=0 ; i<K ; i++ ) {
const int s = rand() % N;
const int t = rand() % N;
ftyp d = 0;
for( int j=0 ; j<L ; j++ )
if( rand() % 1 )
d += 1.0 / ( 1<<(j+1) );
cftyp tmp = v[s] * d;
v[s] -= tmp;
v[t] += tmp;
}
}
void test( cftyp v[N] ) {
ftyp sum = 0;
for( int i=0 ; i<N ; i++ )
sum += v[i];
printf( "%0.20lf %s\n" , sum-1.0, sum==1.0?"true":"false" );
}
int main(int argc, char *argv[]) {
srand(time(NULL));
ftyp v[N];
init(v);
compute(v);
test(v);
return 0;
}
Wynik na kompie/kompilatorze przy którym akurat siedzę:
0.00000000000000000000 true
Pozdrawiam
Następne wpisy z tego wątku
- 27.03.13 08:42 Tomasz Kaczanowski
- 27.03.13 08:47 Tomasz Kaczanowski
- 27.03.13 08:54 AK
- 27.03.13 08:55 AK
- 27.03.13 09:42 M.M.
- 27.03.13 09:47 M.M.
- 27.03.13 10:10 Tomasz Kaczanowski
- 27.03.13 10:25 M.M.
- 27.03.13 10:28 Michoo
- 27.03.13 11:16 AK
- 27.03.13 11:20 darekm
- 27.03.13 12:21 M.M.
- 27.03.13 12:39 AK
- 27.03.13 13:05 M.M.
- 27.03.13 14:32 AK
Najnowsze wątki z tej grupy
- Can you activate BMW 48V 10Ah Li-Ion battery, connecting to CAN-USB laptop interface ?
- We Wrocławiu ruszyła Odra 5, pierwszy w Polsce komputer kwantowy z nadprzewodzącymi kubitami
- Ada-Europe - AEiC 2025 early registration deadline imminent
- John Carmack twierdzi, że gdyby gry były optymalizowane, to wystarczyły by stare kompy
- Ada-Europe Int.Conf. Reliable Software Technologies, AEiC 2025
- Linuks od wer. 6.15 przestanie wspierać procesory 486 i będzie wymagać min. Pentium
- ,,Polski przemysł jest w stanie agonalnym" - podkreślił dobitnie, wskazując na brak zamówień.
- Rewolucja w debugowaniu!!! SI analizuje zrzuty pamięci systemu M$ Windows!!!
- Brednie w wiki - hasło Dehomag
- Perfidne ataki krakerów z KRLD na skrypciarzy JS i Pajton
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- Instytut IDEAS może zacząć działać: "Ma to być unikalny w europejskiej skali ośrodek badań nad sztuczną inteligencją."
- U nas propagują modę na SI, a w Chinach naukowcy SI po kolei umierają w wieku 40-50lat
- C++. Podróż Po Języku - komentarz
Najnowsze wątki
- 2025-07-05 Shrek, wzrośnie cena prądu :-)
- 2025-07-05 Warszawa => Senior Frontend Developer (React + React Native) <=
- 2025-07-05 Środa Wielkopolska => SAP FI/CO Internal Consultant <=
- 2025-07-05 Środa Wielkopolska => SAP FI/CO Konsultant wewnętrzny <=
- 2025-07-06 Płacić za wodę ze studni
- 2025-07-06 Kolejny inżynier...
- 2025-07-05 UE ustępuje pod presją Trumpa. Socjaldemokraci: oddaliśmy miliardy
- 2025-07-04 od kiedy można obrażać Plastusia i do kiedy Batyra?
- 2025-07-05 UE ustępuje pod presją Trumpa. Socjaldemokraci: oddaliśmy miliardy
- 2025-07-03 Trybik
- 2025-07-04 Renault Symbioz
- 2025-07-04 Architektura IIIRP: Wyjątkowa, a prymitywniejsza niż stodoła pod zaborami
- 2025-07-04 Warszawa => International Freight Forwarder <=
- 2025-07-04 Wrocław => SAP ABAP Developer <=
- 2025-07-04 Warszawa => Mid/Senior IT Recruiter <=