[Algorytmy] Implementacja potęgowania x^y

[dzindzinier]

Witam serdecznie.

Wczoraj zderzyłem się twardo z ziemią. Biblioteka matematyki stałoprzecinkowej w moim DSP nie obsługuje żadnego potęgowania za wyjątkiem \(\displaystyle{ e^{x}}\). Potęgowania zmiennoprzecinkowego użyć nie mogę ze względu na wydajność.

Czy szanowni Koledzy/Koleżanki Matematycy mogliby naprowadzić na jakiekolwiek rozwiązanie ?

Myślałem o aproksymacji potęgowania \(\displaystyle{ x^{y}}\) szeregiem. Czy mogę liczyć na podanie szeregu o ile taka aproksymacja jest w ogóle możliwa ?

Czy znane są jakieś inne metody ? Pobieżnie zagłębiłem się w źródła znajdujące się w bibliotece glibc lecz nie tłumaczą one w żaden sposób algorytmu postępowania...


Z góry dziękuję.
Bartek

[tomkoder]

Może to Ci wystarczy?

[dzindzinier]

Cześć.

Niestety to jest jeden z pierwszych algorytmów które sprawdzałem. Cała idea polega na tym aby podnosić dowolną liczbę rzeczywistą do dowolnej rzeczywistej potęgi. Prawie jak na zasadzie funkcji pow() tylko z liczbami stałoprzecinkowymi z ułamkiem. Liczby ułamkowe stałoprzecinkowe stosuje się tam gdzie nie ma sprzętowej jednostki zmiennoprzecinkowej. Zasadę działania można złapać wyszukując "IQmath stellaris" albo "iqmath c28xx".

Pozdrawiam!

[Zordon]

\(\displaystyle{ x^y=e^{y\ln x}}\)