Übersetzungsfehler in GCC

[Bremer (Gast)]

Hallo,

ich bin gerade dabei eine FFT zu schreiben. Dazu müssen einige Summen 
und Differnzen gebildet werden. Leider funktioniert das beim Compilieren 
nicht. GCC verwirft in der Optimierung die Differenzenbildung. Ich weiß, 
nur leider nicht warum. Der Code in C lautet:

temp.real = fft[0].real;
  fft[0].real = (temp.real + fft[2].real)*8;
  fft[2].real = (temp.real - fft[2].real)*8;
  
  temp.real = fft[1].real;
  fft[1].real = temp.real + fft[3].real;
  fft[1].real *= 8;
  fft[3].imag = temp.real - fft[3].real;
  fft[3].imag *= 8;
  
  fft[4].real *= 8;
  fft[6].imag *= 8;
  // fft[4] = fft[4] + fft[6]
  // fft[6] = fft[4] - fft[6]
  
  temp.real = fft[5].real;
  fft[5].real = fft[5].real - fft[7].imag;  // fft[5] = fft[5] + fft[7]
  fft[5].real *= w4.real;
  fft[7].imag = temp.real + fft[7].imag;  // fft[7] = -fft[5] + fft[7]
  fft[7].imag *= w4.real;

und in Assembler wird er erste Block zu:

  temp.real = fft[0].real;
000003C8  LDD R6,Y+1    Load indirect with displacement 
000003C9  LDD R7,Y+2    Load indirect with displacement 
  fft[0].real = (temp.real + fft[2].real)*8;
000003CA  LDD R4,Y+9    Load indirect with displacement 
000003CB  LDD R5,Y+10    Load indirect with displacement 
000003CC  MOVW R20,R6    Copy register pair 
000003CD  ADD R20,R4    Add without carry 
000003CE  ADC R21,R5    Add with carry 
000003CF  LSL R20    Logical Shift Left 
000003D0  ROL R21    Rotate Left Through Carry 
000003D1  LSL R20    Logical Shift Left 
000003D2  ROL R21    Rotate Left Through Carry 
000003D3  LSL R20    Logical Shift Left 
000003D4  ROL R21    Rotate Left Through Carry

Kann mir einer sagen, welchen Fehler ich mache oder wie ich das Problem 
umgehen kann?

Danke im Voraus für die Antworten.

[Jim M. (turboj)]

Ich sehe keinen Fehler im Assembler Code. Der macht ganz genau das was 
Du in C hingeschrieben hast: Eine Addition und einmal Multiplizieren mit 
8 in Form von 3 Left-Shifts, alles in 16 Bit.

[Bremer (Gast)]

ja die Addition passt, aber die Subtraktion fehlt. Vielleicht war es 
missverständlich. Hier sind noch ein paar Zeilen mehr:

  temp.real = fft[0].real;
000003C8  LDD R6,Y+1    Load indirect with displacement 
000003C9  LDD R7,Y+2    Load indirect with displacement 
  fft[0].real = (temp.real + fft[2].real)*8;
000003CA  LDD R4,Y+9    Load indirect with displacement 
000003CB  LDD R5,Y+10    Load indirect with displacement 
000003CC  MOVW R20,R6    Copy register pair 
000003CD  ADD R20,R4    Add without carry 
000003CE  ADC R21,R5    Add with carry 
000003CF  LSL R20    Logical Shift Left 
000003D0  ROL R21    Rotate Left Through Carry 
000003D1  LSL R20    Logical Shift Left 
000003D2  ROL R21    Rotate Left Through Carry 
000003D3  LSL R20    Logical Shift Left 
000003D4  ROL R21    Rotate Left Through Carry 
  temp.real = fft[1].real;
000003D5  LDD R10,Y+5    Load indirect with displacement 
000003D6  LDD R11,Y+6    Load indirect with displacement 
  fft[1].real = temp.real + fft[3].real;
000003D7  LDD R8,Y+13    Load indirect with displacement 
000003D8  LDD R9,Y+14    Load indirect with displacement 
000003D9  MOVW R22,R10    Copy register pair 
000003DA  ADD R22,R8    Add without carry 
000003DB  ADC R23,R9    Add with carry 
  fft[1].real *= 8;
000003DC  LSL R22    Logical Shift Left 
000003DD  ROL R23    Rotate Left Through Carry 
000003DE  LSL R22    Logical Shift Left 
000003DF  ROL R23    Rotate Left Through Carry 
000003E0  LSL R22    Logical Shift Left 
000003E1  ROL R23    Rotate Left Through Carry 
  fft[4].real *= 8;
000003E2  LDD R26,Y+17    Load indirect with displacement 
000003E3  LDD R27,Y+18    Load indirect with displacement 
000003E4  LSL R26    Logical Shift Left 
000003E5  ROL R27    Rotate Left Through Carry 
000003E6  LSL R26    Logical Shift Left 
000003E7  ROL R27    Rotate Left Through Carry 
000003E8  LSL R26    Logical Shift Left 
000003E9  ROL R27    Rotate Left Through Carry 
  fft[6].imag *= 8;
000003EA  LDD R30,Y+27    Load indirect with displacement 
000003EB  LDD R31,Y+28    Load indirect with displacement 
000003EC  LSL R30    Logical Shift Left 
000003ED  ROL R31    Rotate Left Through Carry 
000003EE  LSL R30    Logical Shift Left 
000003EF  ROL R31    Rotate Left Through Carry 
000003F0  LSL R30    Logical Shift Left 
000003F1  ROL R31    Rotate Left Through Carry 
  temp.real = fft[5].real;

[Karl H. (kbuchegg)]

Bremer schrieb:
> ja die Addition passt, aber die Subtraktion fehlt. Vielleicht war es
> missverständlich. Hier sind noch ein paar Zeilen mehr:

Noch mehr.

Denn da fehlt noch mehr. Auch das Zurückschreiben ins Array hab ich noch 
nicht gesehen.
Irgendwo muss das Zeugs ja sein.

[Gerald M. (gerald_m17)]

Vielleicht benutzt du die berechneten Variablen nicht, weshalb sie 
wegoptimiert werden

[Bremer (Gast)]

Guten Morgen,

hier ist der gesammte Quellcode der Funktion.

[Mw E. (Firma: fritzler-avr.de) (fritzler)]

blind?
Such in deiner .s Datei nach
"fft[2].real = (temp.real - fft[2].real)*8;"
und du wirst es finden.
Er hats einfach nur umsortiert.
Da er in der nachfolgenden Rechnung nicht auf "fft[2].real" zugrift 
berechnet er es eben irgendwo.