Quellcode für sin(x) in Delphi

[Wolfgang Mix]

Mich interessiert, ob die Sinusfunktion in Delphi Pascalcode oder
Assemblercode ist. Wenn Pascalcode, würde mich der Quellcode interessieren.
Wenn das ein Staatsgeheimnis ist, genügt mir die Antwort, ob die Funktion
als Taylorreihe implementiert ist.

Habe leider nur eine PE2005

Gruß

Wolfgang

[sirius]

Ähm, Es wird die Sinusfunktion der FPU genutzt :gruebel:

Edit: Und wie die das macht, frag mal deinen Prozessorhersteller!

[Medium]

Imho wird eine Kombi aus LUT und Interpolation genutzt, die für die maximale Genauigkeit grad ausreichen.

[Wolfgang Mix]

Naja, die werden wohl eher antworten als Embadingsbums.
Danke für die rasche Antwort.
Meine Frage hat sich damit erledigt.

Grüß

Wolfgang

[gammatester]

Zitat:

Zitat von Wolfgang Mix

genügt mir die Antwort, ob die Funktion
als Taylorreihe implementiert ist.

Wie immer es auch implementiert ist, mit ziemlicher Sicherheit nicht als Taylorreihe. Taylorreihen haben sehr schlechte Approximationseigenschaften, besser sind Tschebyschev- oder ratinale Approximationen (Stichwort Remez-Algorithmus).

Wenn Du brauchbaren Quellcode sucht, kann man den von SUN aus der Freely Distributable LIBM 5.3 netlib/fdlibm oder den aus der Cephes-Bibliothek empfehlen. Im Gegensatz zu Delphi und FPU arbeiten diese auch richtig über den ganzen Definitionsbereich und liefern nicht sin(1e22)=1e22.

[Wolfgang Mix]

@gammatester:

danke für dir Links, werde ich mir 'mal zu Gemüte führen.

[OT] Die hast hin und wieder Blödsinn in Delphi-Funktionen geschildert.
Könntest Du die 'mal bitte gesammelt posten? Wäre nett [/OT]

Gruß

Wolfgang

[gammatester]

Zitat:

Zitat von Wolfgang Mix

[OT] Die hast hin und wieder Blödsinn in Delphi-Funktionen geschildert.
Könntest Du die 'mal bitte gesammelt posten? Wäre nett [/OT]

Habe rein zufällig gestern die erste Ausgabe meiner AMath-Bibliothek veröffentlicht. Amath enthält zZ verbesserte elementare Funktionen als Ersatz für Delphis Math-Funktionen, außerdem eine über den ganzen Bereich gute exp-Funktion (Fehler AMath: Max=0.74 eps, RMS=0.24 eps, Delphi Max=4563.61 eps, RMS=1574.67, nit eps=1.084e-19).

In der nächsten Ausgabe kommen dann die trigonometrischen Funktionen und eine nochmals verbesserte Power-Funktion.