[Thunderx]

Originariamente inviato da altairz

...
if (a[i+1][z] = 0)
...
if (a[i-1][3-z] = 0)
...

credo che un bug sia in questo punto

alla fine quel pezzo l'ho tolto modificando la condizione sul coefficiente.....ecco il listato della parte modificata

       for (z=0;z<3;z++) {
	for (i=z;i<3;i++) {
   c = -(a[i+1][z] / a[z][z]);
                for (j=z;j<4;j++) {
                 
                a[i + 1][j] = c*a[i][j] + a[i + 1][j];
		
                b[i + 1][j] =c*b[i][j]+ b[i + 1][j];
}
}	
}
  for (z=0;z<3;z++) {
	for (i=(3-z);i>0;i--){
                   d = -(a[i-1][(3-z)]) / a[3-z][3-z];
	               for (j=(3-z);j> -1;j--){
                        a[i-1][j]=d*a[i][j]+a[i-1][j];
                        b[i-1][j]=d*b[i][j]+b[i-1][j];
}	
}
}

questa è l'uscita

1.000000-1.0000002.000000-1.000000
0.0000001.000000-1.0000001.000000
0.0000000.0000001.000000-1.000000
0.0000000.0000000.0000001.000000

con ingresso

1 1 1 1
0 1 1 1
0 0 1 1
0 0 0 1