şamandıra boyutları ve çift

iki vektörleri (şamandıra değeri) nokta ürünün farklı makinelerde farklı sonuçlar veren hesaplama:şamandıra boyutları ve çift

6.102435302 (Win7 x64, compiler VS12 version 17.00.50727.1) 
6.140244007 (Win7 x64, MinGW, gcc version 5.3.0)

kodudur:

#include <iostream> 
#include <iterator> 
#include <fstream> 
#include <vector> 
#include <iomanip> 
#include <algorithm> 

int main(int argc, char** argv){ 

std::ifstream is("test.txt"); 
std::istream_iterator<float> start(is), end; 
std::vector<float> numbers(start, end); 
std::cout << "Read " << numbers.size() << " numbers" << std::endl; 

float product = 0; 
for (int i = 0; i <= numbers.size() - 1; i++) 
    product += (numbers[i])*(numbers[i]); // += means add to product 
std::cout << std::setprecision(10) << product << std::endl; 

std::cin.get(); 
}

metin.txt : Ben herhalde

, bu şekilde şamandıranın büyüklüğüne ve çift bağlıdır ve bu makineler için aynı değildir. Her iki bilgisayar için de aynı çıkışı yapmak mümkün mü?

ben (ilk sonuç 6,102435302 ile) birinci makineye erişimi yok ama (Numpy ile) aşağıdaki python kodu ile aynı sonucu üretebilir:

test = np.loadtxt(test_file, dtype=np.float32) 
result = test.dot(test)

kaynak

2016-04-04 KTBFFH

Kodunuzu, giriş ve derleyici ayarlarınızı gönderebilir misiniz? – Micka

Ve senin [mcve] nerede? –

Normalde, son birkaç bitde kayan nokta doğruluğu varsayılan olarak çeşitli derleyici optimizasyon ayarları ile garanti edilemez, ancak değerleriniz _so far_ apart Bu sorunun sorunu olduğunu düşünmüyorum. Kodunuzu gönderin. –

fark kadar büyük double yerine float kullanılarak açıklanmalıdır. Kodunuzdaki gerçek hataları arayın. Ya da hesabınız derece dengesizdir, bu durumda neler olup bittiğini anlayana kadar nelerin olup bittiğini incelemenize gerek yoktur.

Her iki derleyici için de aynı çıktıyı almak kolaydır - sonucu sıfır olarak ayarlamanız yeterlidir. Ama ne istiyorsan sonucu sonucunu elde etmektir. Yanlış bir sonuç, güvenilemeyecek bir sonuç var ve hangisinin hangisi olduğunu bilmiyorsunuz. Sonuçları aynı yapmak sadece bu durumu kapsayacak, ancak herhangi bir problemi çözmeyecektir.

kaynak

2016-04-04 14:03:49 gnasher729

Büyük olasılıkla bu, büyüklükte çok yakın ama kutupsallıkta tersi olan büyük ürün terimlerinin sonucudur. Birleştirildiğinde, göreli hata artar. Bu genellikle, hastalıklı matrislerle de görülür. – doug

@doug: terimlerin hepsi kare olsa da, hepsi olumlu. –

Ah, kod ve verileri eklediğiniz için teşekkürler. Bu oldukça yaygındır. Bir akümülatöre uzun bir sayı listesi eklediğinizde, başlangıç verileri, sondaki verilerden daha doğru bir şekilde eklenir, çünkü akümülatör büyüklükte daha büyük olur. 1.333333 ve .0000333333 eklemeyi düşünün. Veriler, tam makine hassaslığına yuvarlanmış ve kesilmediğinde, ancak kırpılmışsa, daha hızlı eklenmiş bir yanlılık alabilirsiniz. Şamandıraların vektörleri üzerinde işlemler biriktirirken daima bir çift kullanırım. – doug

Muhtemelen VS12 otomatik döngüyü vectorized (ve sonra sonucu yanlış yazdıysanız) görünüyor.

sen şekilde vektörleştirilmiş döngü çalıştırırsanız:

float product = 0; 
for (int i = 0; i < numbers.size(); i+= 4) 
{ 
    __m128 val = *(__m128*)(&numbers[i]); 
    auto res = _mm_dp_ps (val, val, 255); 
    float result; 
    _mm_store_ss(&result, res); 

    product += result; 
}

Sonra çıkmak sonucudur:

6,14024353

Bu ama görünüşe 6.102435302 senin ilk sonucu olarak aynıdır transkripsiyon yaparken 4'u kaçırmışsınız gibi.

En azından en iyi açıklaması ile ortaya çıkabilirim. Zaten bu soruda çok uzun süre geçirdim :-)

kaynak

2016-04-04 15:33:48

Makul bir hipoteze yol açan çok güzel bir dedektif çalışması! –

şamandıra boyutları ve çift

cevap

İlgili konular