Svenska ElektronikForumet

superx skrev:Min erfarenhet är tyvärr att kompilatorerna är rätt korkade med sånt här. Läs här t.ex.
http://mekonik.wordpress.com/2009/03/18 ... plication/

#include <stdint.h>

uint32_t foo(uint16_t a, uint16_t b)
{
  return (uint32_t)a * (uint32_t)b;
}

$ avr-gcc -mmcu=atmega48 -O2 -S multest.c -o -
        .file   "multest.c"
__SREG__ = 0x3f
__SP_H__ = 0x3e
__SP_L__ = 0x3d
__CCP__ = 0x34
__tmp_reg__ = 0
__zero_reg__ = 1
        .text
.global foo
        .type   foo, @function
foo:
/* prologue: function */
/* frame size = 0 */
/* stack size = 0 */
.L__stack_usage = 0
        movw r20,r24
        movw r18,r22
        rcall __umulhisi3
/* epilogue start */
        ret
        .size   foo, .-foo

$ avr-objdump -d /usr/lib/gcc/avr/4.5.3/avr4/libgcc.a
.
.
00000000 <__umulhisi3>:
   0:   53 9f           mul     r21, r19
   2:   c0 01           movw    r24, r0
   4:   42 9f           mul     r20, r18
   6:   b0 01           movw    r22, r0
   8:   52 9f           mul     r21, r18
   a:   70 0d           add     r23, r0
   c:   81 1d           adc     r24, r1
   e:   11 24           eor     r1, r1
  10:   91 1d           adc     r25, r1
  12:   34 9f           mul     r19, r20
  14:   70 0d           add     r23, r0
  16:   81 1d           adc     r24, r1
  18:   11 24           eor     r1, r1
  1a:   91 1d           adc     r25, r1
  1c:   08 95           ret

superx skrev:Testade lite med avr-gcc här.

Kod: Markera allt

uint16_t interp1(const uint8_t x, const uint8_t A, const uint8_t B, const uint8_t M)
{
  uint8_t AB = A-B;
  uint16_t V = B + (AB)*x/M;

  return V;
}

Ger 1 multiplikation

Kod: Markera allt

uint16_t interp2(const uint8_t x, const uint8_t A, const uint8_t B, const uint8_t M)
{
  uint16_t V = B + (A-B)*x/M;

  return V;
}

Ger 3 multiplikationer. Jag antar att det beror på att man i det övre fallet själv tar ansvar för att A-B inte slår runt. Går detta att utläsa av standarden?

EDIT: Ändrade AB från signed till unsigned

superx skrev:Spännande exempel! Gäller att se upp när man blandar signed och unsigned verkar det som.

Testade lite med avr-gcc här.

Kod: Markera allt

uint16_t interp1(const uint8_t x, const uint8_t A, const uint8_t B, const uint8_t M)
{
  uint8_t AB = A-B;
  uint16_t V = B + (AB)*x/M;

  return V;
}

Ger 1 multiplikation

Kod: Markera allt

uint16_t interp2(const uint8_t x, const uint8_t A, const uint8_t B, const uint8_t M)
{
  uint16_t V = B + (A-B)*x/M;

  return V;
}

Ger 3 multiplikationer. Jag antar att det beror på att man i det övre fallet själv tar ansvar för att A-B inte slår runt. Går detta att utläsa av standarden?

EDIT: Ändrade AB från signed till unsigned

superx skrev:Jag har egentligen mest haft problem av den här typen på plattformar med extra stor ackumulator och MAC-instruktioner (multiply-accumulate), där man gärna vill behålla alla delresultaten i ackumulatorn för att undvika rundningar.