Schneller Wurzeln ziehen als sqrtf das kann?

[jaque]

Hallo,

habe mich mal an eine eigene sqrt-Funktion für floats gemacht. Bin allerdings auf Dauer nicht annähernd an die Geschwindigkeit von sqrtf aus der Standardlib rangekommen. Gibts da nen schön schnellen Algorithmus?

Gruß
jaque

[Gast]

Ich verschiebe das mal nach "Rund um die Programmierung", da ein Algorithmus selbst recht wenig mit der zur Implementierung benutzten Sprache zu tun hat.

Wenn dein Prozessor eine eigene FPU hat, würde ich mir eher weniger Hoffnung auf einen schnelleren Algorithmus machen

[C++ Forumbot]

Dieser Thread wurde von Moderator/in TactX aus dem Forum ANSI C in das Forum Rund um die Programmierung verschoben.

Im Zweifelsfall bitte auch folgende Hinweise beachten:
C/C++ Forum :: FAQ - Sonstiges :: Wohin mit meiner Frage?

Dieses Posting wurde automatisch erzeugt.

[rapso]

jaque schrieb:

Hallo,

habe mich mal an eine eigene sqrt-Funktion für floats gemacht. Bin allerdings auf Dauer nicht annähernd an die Geschwindigkeit von sqrtf aus der Standardlib rangekommen. Gibts da nen schön schnellen Algorithmus?

Gruß
jaque

du wirst, wie schon gesagt, nichts schnelleres als die fpu macht machen koennen bei gleicher genauigkeit. ansonsten bastel dir ein LUT.

[geh heim]

float sqrt (float f)
{
	if (f == 0.0f)
		return 0.0f;

	float f1 = f;
	unsigned long *ptr = (unsigned long*)&f;
	*ptr = (0xbe6f0000-*ptr)>>1;

	return f1*f;
}

[rapso]

geh heim schrieb:

float sqrt (float f)
{
	if (f == 0.0f)
		return 0.0f;

	float f1 = f;
	unsigned long *ptr = (unsigned long*)&f;
	*ptr = (0xbe6f0000-*ptr)>>1;

	return f1*f;
}

ist langsammer als die fpu.

[jaque]

geh heim schrieb:

float sqrt (float f)
{
	if (f == 0.0f)
		return 0.0f;

	float f1 = f;
	unsigned long *ptr = (unsigned long*)&f;
	*ptr = (0xbe6f0000-*ptr)>>1;

	return f1*f;
}

Nett, danke. Gibts ne Erklärung dazu?