------- Comment #4 from dominiq at lps dot ens dot fr 2007-03-11 13:16 ------- Note that the bug affect lapack/blas/{s,d}scal.f. Note also that in the *optimized* dumps of lapack/blas/{c,z}scal.f there are several:
Invalid sum of incoming frequencies xxxx, should be xxxx cscal.f: ... Invalid sum of incoming frequencies 2224, should be 2116 <L19>:; D.1073 = (<unnamed type>) ivtmp.47; CR.63 = MEM[index: D.1073]; CI.64 = MEM[index: D.1073, offset: 4]; MEM[index: D.1073] = CR.61 * CR.63 - CI.62 * CI.64; MEM[index: D.1073, offset: 4] = CR.61 * CI.64 + CI.62 * CR.63; ivtmp.44 = ivtmp.44 + 1; ivtmp.47 = ivtmp.47 + 8; if ((int4) ivtmp.44 == D.968) goto <L14>; else goto <L19>; Invalid sum of incoming frequencies 2322, should be 2430 <L14>:; return; and zscal.f: ... Invalid sum of incoming frequencies 10511, should be 10000 <L24>:; D.1095 = (<unnamed type>) ivtmp.59; CR.21 = MEM[index: D.1095]; CI.22 = MEM[index: D.1095, offset: 8]; MEM[index: D.1095] = CR.64 * CR.21 - CI.65 * CI.22; MEM[index: D.1095, offset: 8] = CR.64 * CI.22 + CI.65 * CR.21; ivtmp.53 = ivtmp.53 + 1; ivtmp.59 = ivtmp.59 + ivtmp.62; if ((int4) ivtmp.53 == D.967) goto <L16>; else goto <L24>; <L29>:; CR.19 = REALPART_EXPR <*za>; CI.20 = IMAGPART_EXPR <*za>; ivtmp.50 = (<unnamed type>) &REALPART_EXPR <(*zx)[0]>; ivtmp.47 = 0; Invalid sum of incoming frequencies 4723, should be 4493 <L23>:; D.1079 = (<unnamed type>) ivtmp.50; CR.25 = MEM[index: D.1079]; CI.26 = MEM[index: D.1079, offset: 8]; MEM[index: D.1079] = CR.19 * CR.25 - CI.20 * CI.26; MEM[index: D.1079, offset: 8] = CR.19 * CI.26 + CI.20 * CR.25; ivtmp.47 = ivtmp.47 + 1; ivtmp.50 = ivtmp.50 + 16; if ((int4) ivtmp.47 == D.967) goto <L16>; else goto <L23>; Invalid sum of incoming frequencies 4420, should be 5160 <L16>:; return; -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31030