Intressant artikel om Quake3:s InvSqrt()|Programmering / Grafisk formgivning|Forum|Nordichardware

Intressant nog så fungerar den koden inte med modernare C-kompilatorer (till exempel gcc-4.1). Problemet är de roliga castarna mellan int och float; de följer inte riktigt C-standarden och gcc anser därför att den har rätt att optimera hej vilt. Om man kan assembler kan man titta på en objdump och se att kompilatorn mer eller mindre har slängt bort en massa kod.

#include <stdio.h>

	

	static float invsqrt2(float x)

	{

	        float xhalf = 0.5f * x;

	        union {

	                int i;

	                float f;

	        } u = { .f = x };

	

	        u.i = 0x5f3759df - (u.i >> 1);

	        x = u.f;

	        x = x * (1.5f - xhalf * x * x);

	        return x;

	}

	

	

	static float invsqrt1(float x)

	{

	        float xhalf = 0.5f * x;

	        int i = *(int*)&x;

	        i = 0x5f3759df - (i >> 1);

	        x = *(float*)&i;

	        x = x * (1.5f - xhalf * x * x);

	        return x;

	}

	

	int main(int argc, char **argv)

	{

	        printf("version 1: 1/sqrt(25) = %fn", invsqrt1(25.0f));

	        printf("version 2: 1/sqrt(25) = %fn", invsqrt2(25.0f));

	}

	

	$ gcc --version

	gcc (GCC) 4.1.1

	Copyright (C) 2006 Free Software Foundation, Inc.

	This is free software; see the source for copying conditions.  There is NO

	warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

	$ gcc -march=athlon64 -msse3  -O3  -o foo foo.c

	$ ./foo

	version 1: 1/sqrt(25) = 0.000000

	version 2: 1/sqrt(25) = 0.199690