------- Comment #14 from dir at lanl dot gov 2007-06-21 12:57 ------- >What is actually the expected result? Depending on the compiler and compiler >setting, I get completely different results for the second triangular matrix. >(The first matrix remains always the same.)
What the program does in this section is multiply the diagonal of the matrix by 2 with the line - tt=tt+tt and so (0.8000E+01 -> 0.1600E+02), (0.1000E+02 -> 0.2000E+02) and (0.1300E+02 -> 0.2600E+02 ) and the other three elememts should stay the same. In debug mode, most compilers will give the correct answer, but the addition is sometimes being optmized away. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=32393