[Megoldva: -ffloat-store] Egyenletrendszer megoldás, i686 vs x86-64

A lenti program megoldana egy lineáris egyenletrendszert. A megoldás: 1 2 3 4. A kód olyan, amilyen. A probléma az, hogy i686-on rossz lesz az eredmény. Nevezetesen: 1.07407 2 3 3.48148 . x86_64-en és sparc-on jól futott/fordult le. Érdekesség: ha kihagyom az oszlopcserés részt(HIBA innen-HIBA eddig között), akkor jó lesz i686-on is. Jo akkor is, ha 7 helyett 6 es 16 helyett 15 van. Mellékelem a kódot, ha valakinek van egy kis ideje, akkor futtassa le, és írja meg, és hogy működött-e. Írja le hozzászólásban, hogy milyen környezetben használta. Bármilyen ötletet szívesen veszek.

Köszönöm a segítséget:
István

main.cpp

#include<iostream>
#include "gauss.h"

int main()
{
    mat A={{1.0,1.0,1.0,1.0},
           {1.0,3.0,1.0,4.0},
           {7.0,1.0,1.0,1.0},
           {3.0,3.0,1.0,3.0}};
    vect b={10.0,26.0,16.0,24.0};
    vect x;
    gauss(&A,&b,&x);
    for (int i=0;i<N;++i)  std::cout << x[i] <<" ";
    std::cout << std::endl;
    return 0;
}

gauss.h:

#ifndef GAUSS_H
#define GAUSS_H

#include<iostream>
#include<cmath>

#define N 4

typedef double vect[N];
typedef double mat[N][N];


int gauss(mat *A, vect *b, vect *x)
{
    
    int i,j,k;
    int xn[N];
    int mx,my;
    double d;
        
    for(i=0;i<N;i++) xn[i]=i;
    for(i=0;i<N;i++)
    {
        
        // foelem kivalasztas
        mx=i;
        my=i;
        for(j=i;j<N;j++)
        {
            for(k=i;k<N;k++)
            {
                if(fabs((*A)[j][k])>fabs((*A)[my][mx]))
                {
                    my=j;
                    mx=k;                    
                }
            }
        }

        if( my != i) // sorcsere
        {
            for(j=i;j<N;j++)
            {
                d=(*A)[i][j];
                (*A)[i][j]=(*A)[my][j];
                (*A)[my][j]=d; 
            }
            d=(*b)[i];
            (*b)[i]=(*b)[my];
            (*b)[my]=d;        
        } 

/* HIBA innen */
        
       if( mx != i ) // oszlopcsere
        {
            for(j=i;j<N;j++)
            {
                d=(*A)[j][i];
                (*A)[j][i]=(*A)[j][mx];
                (*A)[j][mx]=d;    
            }
            k=xn[i];
            xn[i]=xn[mx];
            xn[mx]=k;
        }
/* HIBA Eddig  */
        if( (*A)[i][i] == 0.0)
        {
            std::cerr << "Gauss Error" << std::endl;
            return 10;
        }
        d=(*A)[i][i];
        
        for(j=i;j<N;j++) (*A)[i][j]=(*A)[i][j]/d;
        (*b)[i]=(*b)[i]/d;
    
        if( i < N-1)
        {
            for(k=i+1;k<N;k++)
            {
                d=(*A)[k][i];
                for(j=i;j<N;j++)
                {
                    (*A)[k][j]=(*A)[k][j]-d*(*A)[i][j]/(*A)[i][i];
                }   
                (*b)[k]=(*b)[k]-d*((*b)[i]);        
            }
        }
    }
    
    for(i=N-2;i>0;i--)
    {
        for(j=N-1;j>i;j--)
        {
            (*b)[i]=(*b)[i]-(*A)[j][i]*(*b)[j];
        }    
    }
    
    for(i=N-2;i>=0;i--)
    {
        d=0.0;
        for(j=N-1;j>i;j--)
        {
            d+=(*A)[i][j]*(*b)[j];
        }    
        (*b)[i]-=d;
    }
    
    for(i=0;i<N;i++)
    {
        (*x)[xn[i]]=(*b)[i];
    }
    return 0;

}

#endif