求C#大神帮忙用C#重写下面代码，要求编译正确。

Question

float InvSqrt(float x) {
float xhalf = 0.5f*x;
int i = *(int*)&x; // get bits for floating VALUE
i = 0x5f375a86- (i>>1); // gives initial guess y0
x = *(float*)&i; // convert bits BACK to
float x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
return x;
}
备注：函数的输入与返回值要求是decimal类的（不知道在原则上是否可行），请指教！

Answer 1

快速平方根倒数
decimal是128bit浮点数，原则上可行。
1、要找一个类似0x5f375a86的魔术数字R就可以了。
2、这个R也要是128bit的整型，C#也不支持。可以用long[2]来实现减法和位移操作，
C#/.net甚至PC/x86先天不支持128bit的整型，要先实现整型的指针、位移、减法操作，
就意味这个"快速"算法作用于128bit时失去了快速的优势。
而用普通算法，这个函数值f(x)等于1除以x的平方根，得到一样的值。
所以作用于128bit时是否真快速，有待研究吧

对于32bit和64bit，可以改成：
//32bit
static public unsafe float InvSqrt(float x) //加unsafe关键字
{
float xhalf = 0.5f * x;
int i = *(int*)&x;
i = 0x5f375a86 - (i >> 1); // 魔术数字
x = *(float*)&i;
x *= (1.5f - xhalf * x * x);
return x;
}
//64bit
static public unsafe double InvSqrt(double x)
{
double xhalf = 0.5 * x;
long i = *(long*)&x;
i = 0x5fe6eb50c7aa19f9 - (i >> 1);
x = *(double*)&i;
x *= (1.5 - xhalf * x * x);
return x;
}
这2个快速算法90年代时值得称赞。在21世纪的PC已经失去了快速优势，x86的cpu集成了SSE的机器指令 rsqrtss 能更快速更精确地求得该值。

Answer 2

decimal InvSqrt(float x)
{
float xhalf = 0.5f * x;
int i = (int)x; // get bits for floating VALUE
i = 0x5f375a86 - (i >> 1); // gives initial guess y0
x = i; // convert bits BACK to
return Convert.ToDecimal(x * (1.5f - xhalf * x * x)); // Newton step, repeating increases accuracy
}

追问

请问能不能把函数入口的类型也设为decimal，而且函数中间运算时也最好都是decimal的，因为我希望得到尽可能高的精度，（本人初学C#，其中数据类型的原委不是很清楚，如果思路有错，请指正）