piyush01123 · March 16, 2023 01:36 · piyush01123 · Aug 11, 2021
diff --git a/fisr.cpp b/fisr.cpp
 float fisr(float n){
 	long x = 0x5f3759df - (*(long *)&n >> 1); //magic number = 1.5*2^23*(127-mu); mu represents log(1+x)=x+mu; mu=0.0450465
 	float y = *(float *)&x;  // Solution from bit manipulation
 	y = 1.5F*y - 0.5F*n*y*y*y; // Newton improvement 1
        y = 1.5F*y - 0.5F*n*y*y*y; // Newton improvement 2
 	return y;
 }
	float fisr(float n){
	long x = 0x5f3759df - ((long )&n >> 1); //magic number = 1.52^23(127-mu); mu represents log(1+x)=x+mu; mu=0.0450465
	float y = (float )&x; // Solution from bit manipulation
	y = 1.5Fy - 0.5Fnyy*y; // Newton improvement 1
	y = 1.5Fy - 0.5Fnyy*y; // Newton improvement 2
	return y;
	}