------- Comment #1 from yuri at tsoft dot com 2005-12-20 05:34 ------- actually it's the defect in this case: result is not used. But runtimes are very different in any case. 44.9s on 4.x.x vs. 0m2.371s on 3.4.4
--- begin corrected testcase ----------------------------------- #include <iostream> using namespace std; typedef float __v2df __attribute__ ((__vector_size__ (16))); typedef __v2df __m128; static __inline __m128 _mm_sub_pd (__m128 __A, __m128 __B) { return (__m128)__builtin_ia32_subps ((__v2df)__A, (__v2df)__B); } static __inline __m128 _mm_add_pd (__m128 __A, __m128 __B) { return (__m128)__builtin_ia32_addps ((__v2df)__A, (__v2df)__B); } static __inline __m128 _mm_setr_ps (float __Z, float __Y, float __X, float __W) { return __extension__ (__m128)(__v2df){ __Z, __Y, __X, __W }; } struct FF { __m128 d; __inline FF() { } __inline FF(__m128 new_d) : d(new_d) { } __inline FF(float f) : d(_mm_setr_ps(f, f, f, f)) { } __inline FF operator+(FF other) { return (FF(_mm_add_pd(d,other.d))); } __inline FF operator-(FF other) { return (FF(_mm_sub_pd(d,other.d))); } }; float f[1024*1024]; union U { __m128 m; float f[4]; }; int main() { int i; FF gtotal(0.f); for (i = 0; i < 1024*1024; i++) { f[i] = 1.f/(1024*1024 + 10 - i); } FF total(0.f); for (int rpt = 0; rpt < 1000; rpt++) { FF p1(0.f), p2(0.), c; __m128 *pf = (__m128*)f; for (i = 0; i < 1024*1024/4; i++) { FF c(*pf++); total = total + c - p2 + p1; p1 = p2; p2 = c; } gtotal = gtotal + total; } U u; u.m = gtotal.d; cout << (u.f[0]) << endl; } --- end corrected testcase ------------------------------------- -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=25500