GCC f. AVR: SQRT von _Fract

[Thomas (Gast)]

Zum Ziehen der Wurzel aus einer _Fract-Variablen, kann man sqrt() 
nehmen. Dies bläht jedoch den Code um über 500 Bytes auf, die ich nicht 
habe (ATTiny13).

Gibt es eine platzsparende Möglichkeit eine Wurzel zu ziehen (max. 200 
Bytes)?

[Jim M. (turboj)]

Entweder Festkomma + Lookup Table oder größeren µC mit mehr Flash 
verwenden.

[Thomas (Gast)]

Hier ein Algorithmus der aus einer Integer die Wurzel zieht. Er ist 
schnell (kleiner 200 Takte) und kurz (ca. 100 Bytes).


Quelle: 
https://stackoverflow.com/questions/1100090/looking-for-an-efficient-integer-square-root-algorithm-for-arm-thumb2

uint16_t SQRT16(uint16_t op)
{
  uint16_t res = 0;
  uint16_t one = 1 << 14; // The second-to-top bit is set: 1L<<30 for long

  // "one" starts at the highest power of four <= the argument.
  while (one > op)
  {
    one >>= 2;
  }
  
  while (one != 0)
  {
    if (op >= res + one)
    {
      op -= res + one;
      res += 2 * one;
    }
    res >>= 1;
    one >>= 2;
  }
  return (res);
}

[Jim M. (turboj)]

Thomas schrieb:
> Hier ein Algorithmus der aus einer Integer die Wurzel zieht

Dann hat man aber immer noch das Problem aus dem _Fract (16-bittig 
Fließkomma) einen vernüntigen Integer zu machen. Wäre das trivial, hätte 
der OP sicher nicht gefragt.

[Thomas (Gast)]

Jim M. schrieb:
> Thomas schrieb:
>> Hier ein Algorithmus der aus einer Integer die Wurzel zieht
>
> Dann hat man aber immer noch das Problem aus dem _Fract (16-bittig
> Fließkomma) einen vernüntigen Integer zu machen. Wäre das trivial, hätte
> der OP sicher nicht gefragt.

Das ist trivial. Nur Shiften. _Fract-Vaiable z. B. mit  2^14 
multiplizieren, obigen Algorithnus ausführen, Ergebnis durch 2^7 
dividieren.

Beispiel:

0,25 * 2^14 -> 4096
sqrt(4096)  -> 64
64 / 2^7    -> 0,5