csqrt32.cpp

这是我本人编写的一个32位定点小数运算的函数库。对于没有浮点运算器的场合

	3846924535u, 3847509841u, 3848095058u, 3848680186u, 3849265225u, 3849850175u, 
	3850435037u, 3851019809u, 3851604493u, 3852189088u, 3852773595u, 3853358012u, 
	3853942341u, 3854526582u, 3855110734u, 3855694797u, 3856278772u, 3856862659u, 
	3857446457u, 3858030167u, 3858613788u, 3859197322u, 3859780767u, 3860364124u, 
	3860947392u, 3861530573u, 3862113666u, 3862696670u, 3863279587u, 3863862415u, 
	3864445156u, 3865027809u, 3865610374u, 3866192851u, 3866775241u, 3867357543u, 
	3867939757u, 3868521883u, 3869103922u, 3869685874u, 3870267738u, 3870849514u, 
	3871431203u, 3872012805u, 3872594319u, 3873175746u, 3873757086u, 3874338338u, 
	3874919503u, 3875500582u, 3876081573u, 3876662476u, 3877243293u, 3877824023u, 
	3878404666u, 3878985222u, 3879565692u, 3880146074u, 3880726369u, 3881306578u, 
	3881886700u, 3882466736u, 3883046684u, 3883626546u, 3884206322u, 3884786011u, 
	3885365614u, 3885945130u, 3886524559u, 3887103903u, 3887683160u, 3888262330u, 
	3888841415u, 3889420413u, 3889999325u, 3890578151u, 3891156891u, 3891735544u, 
	3892314112u, 3892892594u, 3893470989u, 3894049299u, 3894627523u, 3895205661u, 
	3895783714u, 3896361680u, 3896939561u, 3897517356u, 3898095066u, 3898672690u, 
	3899250228u, 3899827681u, 3900405048u, 3900982330u, 3901559527u, 3902136638u, 
	3902713664u, 3903290604u, 3903867459u, 3904444229u, 3905020914u, 3905597514u, 
	3906174028u, 3906750457u, 3907326802u, 3907903061u, 3908479236u, 3909055325u, 
	3909631330u, 3910207250u, 3910783085u, 3911358835u, 3911934500u, 3912510081u, 
	3913085577u, 3913660988u, 3914236315u, 3914811558u, 3915386715u, 3915961789u, 
	3916536778u, 3917111682u, 3917686502u, 3918261238u, 3918835890u, 3919410457u, 
	3919984940u, 3920559339u, 3921133653u, 3921707884u, 3922282031u, 3922856093u, 
	3923430072u, 3924003966u, 3924577777u, 3925151503u, 3925725146u, 3926298705u, 
	3926872181u, 3927445572u, 3928018880u, 3928592104u, 3929165245u, 3929738302u, 
	3930311275u, 3930884165u, 3931456972u, 3932029695u, 3932602334u, 3933174891u, 
	3933747363u, 3934319753u, 3934892059u, 3935464282u, 3936036422u, 3936608479u, 
	3937180453u, 3937752343u, 3938324151u, 3938895875u, 3939467517u, 3940039075u, 
	3940610551u, 3941181944u, 3941753254u, 3942324481u, 3942895626u, 3943466687u, 
	3944037666u, 3944608563u, 3945179376u, 3945750108u, 3946320756u, 3946891322u, 
	3947461806u, 3948032207u, 3948602526u, 3949172763u, 3949742917u, 3950312989u, 
	3950882978u, 3951452886u, 3952022711u, 3952592454u, 3953162115u, 3953731694u, 
	3954301191u, 3954870605u, 3955439938u, 3956009189u, 3956578358u, 3957147445u, 
	3957716451u, 3958285374u, 3958854216u, 3959422976u, 3959991654u, 3960560251u, 
	3961128766u, 3961697200u, 3962265552u, 3962833822u, 3963402011u, 3963970118u, 
	3964538145u, 3965106089u, 3965673953u, 3966241735u, 3966809436u, 3967377055u, 
	3967944594u, 3968512051u, 3969079427u, 3969646722u, 3970213936u, 3970781069u, 
	3971348121u, 3971915092u, 3972481982u, 3973048791u, 3973615519u, 3974182167u, 
	3974748733u, 3975315219u, 3975881625u, 3976447949u, 3977014193u, 3977580356u, 
	3978146439u, 3978712441u, 3979278363u, 3979844204u, 3980409965u, 3980975645u, 
	3981541245u, 3982106765u, 3982672204u, 3983237563u, 3983802842u, 3984368041u, 
	3984933160u, 3985498198u, 3986063156u, 3986628034u, 3987192833u, 3987757551u, 
	3988322189u, 3988886747u, 3989451226u, 3990015624u, 3990579943u, 3991144182u, 
	3991708341u, 3992272421u, 3992836421u, 3993400341u, 3993964181u, 3994527942u, 
	3995091623u, 3995655225u, 3996218748u, 3996782190u, 3997345554u, 3997908838u, 
	3998472043u, 3999035168u, 3999598214u, 4000161181u, 4000724069u, 4001286877u, 
	4001849607u, 4002412257u, 4002974828u, 4003537320u, 4004099733u, 4004662067u, 
	4005224322u, 4005786499u, 4006348596u, 4006910614u, 4007472554u, 4008034415u, 
	4008596197u, 4009157900u, 4009719525u, 4010281071u, 4010842539u, 4011403927u, 
	4011965238u, 4012526469u, 4013087623u, 4013648698u, 4014209694u, 4014770612u, 
	4015331452u, 4015892213u, 4016452896u, 4017013501u, 4017574027u, 4018134476u, 
	4018694846u, 4019255138u, 4019815352u, 4020375488u, 4020935546u, 4021495526u, 
	4022055427u, 4022615251u, 4023174998u, 4023734666u, 4024294256u, 4024853769u, 
	4025413204u, 4025972561u, 4026531840u, 4027091042u, 4027650166u, 4028209212u, 
	4028768181u, 4029327072u, 4029885886u, 4030444623u, 4031003282u, 4031561863u, 
	4032120367u, 4032678794u, 4033237143u, 4033795415u, 4034353610u, 4034911728u, 
	4035469769u, 4036027732u, 4036585618u, 4037143427u, 4037701159u, 4038258814u, 
	4038816392u, 4039373893u, 4039931317u, 4040488665u, 4041045935u, 4041603129u, 
	4042160245u, 4042717285u, 4043274248u, 4043831135u, 4044387945u, 4044944678u, 
	4045501334u, 4046057914u, 4046614418u, 4047170845u, 4047727195u, 4048283469u, 
	4048839666u, 4049395787u, 4049951832u, 4050507801u, 4051063693u, 4051619509u, 
	4052175248u, 4052730911u, 4053286499u, 4053842010u, 4054397445u, 4054952804u, 
	4055508086u, 4056063293u, 4056618424u, 4057173479u, 4057728458u, 4058283361u, 
	4058838188u, 4059392939u, 4059947615u, 4060502215u, 4061056739u, 4061611187u, 
	4062165560u, 4062719857u, 4063274078u, 4063828224u, 4064382294u, 4064936289u, 
	4065490208u, 4066044052u, 4066597820u, 4067151513u, 4067705131u, 4068258673u, 
	4068812140u, 4069365532u, 4069918848u, 4070472089u, 4071025255u, 4071578346u, 
	4072131362u, 4072684302u, 4073237168u, 4073789959u, 4074342674u, 4074895315u, 
	4075447880u, 4076000371u, 4076552787u, 4077105128u, 4077657394u, 4078209586u, 
	4078761702u, 4079313744u, 4079865711u, 4080417604u, 4080969422u, 4081521165u, 
	4082072834u, 4082624428u, 4083175948u, 4083727393u, 4084278764u, 4084830060u, 
	4085381282u, 4085932430u, 4086483503u, 4087034502u, 4087585426u, 4088136277u, 
	4088687053u, 4089237755u, 4089788383u, 4090338937u, 4090889416u, 4091439822u, 
	4091990153u, 4092540411u, 4093090594u, 4093640704u, 4094190740u, 4094740702u, 
	4095290589u, 4095840404u, 4096390144u, 4096939811u, 4097489403u, 4098038923u, 
	4098588368u, 4099137740u, 4099687038u, 4100236263u, 4100785414u, 4101334491u, 
	4101883495u, 4102432426u, 4102981283u, 4103530067u, 4104078777u, 4104627414u, 
	4105175977u, 4105724468u, 4106272885u, 4106821229u, 4107369499u, 4107917697u, 
	4108465821u, 4109013872u, 4109561850u, 4110109756u, 4110657588u, 4111205347u, 
	4111753033u, 4112300646u, 4112848186u, 4113395653u, 4113943048u, 4114490369u, 
	4115037618u, 4115584794u, 4116131898u, 4116678928u, 4117225886u, 4117772771u, 
	4118319584u, 4118866324u, 4119412992u, 4119959587u, 4120506109u, 4121052559u, 
	4121598937u, 4122145242u, 4122691474u, 4123237635u, 4123783723u, 4124329739u, 
	4124875682u, 4125421553u, 4125967352u, 4126513079u, 4127058734u, 4127604316u, 
	4128149826u, 4128695265u, 4129240631u, 4129785925u, 4130331148u, 4130876298u, 
	4131421376u, 4131966383u, 4132511317u, 4133056180u, 4133600971u, 4134145690u, 
	4134690337u, 4135234913u, 4135779417u, 4136323849u, 4136868210u, 4137412499u, 
	4137956716u, 4138500862u, 4139044936u, 4139588939u, 4140132870u, 4140676730u, 
	4141220518u, 4141764235u, 4142307881u, 4142851455u, 4143394958u, 4143938390u, 
	4144481750u, 4145025040u, 4145568258u, 4146111404u, 4146654480u, 4147197485u, 
	4147740418u, 4148283281u, 4148826072u, 4149368792u, 4149911442u, 4150454020u, 
	4150996528u, 4151538965u, 4152081330u, 4152623625u, 4153165849u, 4153708003u, 
	4154250085u, 4154792097u, 4155334038u, 4155875909u, 4156417709u, 4156959438u, 
	4157501097u, 4158042685u, 4158584202u, 4159125650u, 4159667026u, 4160208332u, 
	4160749568u, 4161290733u, 4161831828u, 4162372853u, 4162913807u, 4163454691u, 
	4163995505u, 4164536248u, 4165076922u, 4165617525u, 4166158058u, 4166698521u, 
	4167238914u, 4167779236u, 4168319489u, 4168859672u, 4169399784u, 4169939827u, 
	4170479800u, 4171019703u, 4171559536u, 4172099299u, 4172638992u, 4173178616u, 
	4173718170u, 4174257654u, 4174797068u, 4175336413u, 4175875687u, 4176414893u, 
	4176954029u, 4177493095u, 4178032091u, 4178571018u, 4179109876u, 4179648664u, 
	4180187383u, 4180726032u, 4181264612u, 4181803122u, 4182341563u, 4182879935u, 
	4183418238u, 4183956471u, 4184494635u, 4185032730u, 4185570756u, 4186108712u, 
	4186646600u, 4187184418u, 4187722167u, 4188259847u, 4188797458u, 4189335001u, 
	4189872474u, 4190409878u, 4190947214u, 4191484480u, 4192021678u, 4192558806u, 
	4193095866u, 4193632858u, 4194169780u, 4194706634u, 4195243419u, 4195780135u, 
	4196316783u, 4196853362u, 4197389873u, 4197926314u, 4198462688u, 4198998993u, 
	4199535229u, 4200071397u, 4200607497u, 4201143528u, 4201679490u, 4202215385u, 
	4202751211u, 4203286968u, 4203822658u, 4204358279u, 4204893832u, 4205429317u, 
	4205964733u, 4206500082u, 4207035362u, 4207570574u, 4208105718u, 4208640794u, 
	4209175802u, 4209710742u, 4210245614u, 4210780419u, 4211315155u, 4211849823u, 
	4212384424u, 4212918956u, 4213453421u, 4213987818u, 4214522147u, 4215056409u, 
	4215590603u, 4216124729u, 4216658787u, 4217192778u, 4217726701u, 4218260557u, 
	4218794345u, 4219328066u, 4219861719u, 4220395305u, 4220928823u, 4221462273u, 
	4221995657u, 4222528973u, 4223062221u, 4223595402u, 4224128516u, 4224661563u, 
	4225194543u, 4225727455u, 4226260300u, 4226793078u, 4227325788u, 4227858432u, 
	4228391008u, 4228923518u, 4229455960u, 4229988336u, 4230520644u, 4231052885u, 
	4231585060u, 4232117167u, 4232649208u, 4233181182u, 4233713089u, 4234244929u, 
	4234776702u, 4235308409u, 4235840048u, 4236371621u, 4236903128u, 4237434568u, 
	4237965941u, 4238497247u, 4239028487u, 4239559660u, 4240090767u, 4240621807u, 
	4241152781u, 4241683688u, 4242214529u, 4242745304u, 4243276012u, 4243806653u, 
	4244337229u, 4244867738u, 4245398180u, 4245928557u, 4246458867u, 4246989111u, 
	4247519289u, 4248049400u, 4248579446u, 4249109425u, 4249639338u, 4250169186u, 
	4250698967u, 4251228682u, 4251758331u, 4252287914u, 4252817431u, 4253346883u, 
	4253876268u, 4254405587u, 4254934841u, 4255464029u, 4255993151u, 4256522207u, 
	4257051198u, 4257580123u, 4258108982u, 4258637775u, 4259166503u, 4259695165u, 
	4260223762u, 4260752292u, 4261280758u, 4261809158u, 4262337492u, 4262865761u, 
	4263393964u, 4263922102u, 4264450175u, 4264978182u, 4265506124u, 4266034001u, 
	4266561812u, 4267089558u, 4267617238u, 4268144854u, 4268672404u, 4269199889u, 
	4269727309u, 4270254664u, 4270781953u, 4271309178u, 4271836337u, 4272363431u, 
	4272890461u, 4273417425u, 4273944324u, 4274471159u, 4274997928u, 4275524633u, 
	4276051273u, 4276577848u, 4277104358u, 4277630803u, 4278157183u, 4278683499u, 
	4279209750u, 4279735936u, 4280262058u, 4280788115u, 4281314107u, 4281840035u, 
	4282365898u, 4282891696u, 4283417430u, 4283943100u, 4284468705u, 4284994245u, 
	4285519721u, 4286045133u, 4286570480u, 4287095763u, 4287620981u, 4288146135u, 
	4288671225u, 4289196251u, 4289721212u, 4290246109u, 4290770942u, 4291295711u, 
	4291820415u, 4292345056u, 4292869632u, 4293394144u, 4293918592u, 4294442976u,
	4294967295u
};

////////////////////////////////////////////////////////////

#pragma warning(disable:4554)
long FixSqrt32_C(long x,long y)
{
	short x1=*((short*)(&x)+1),y1=*((short*)(&y)+1);
	unsigned short x0=*((unsigned short*)(&x)),y0=*((unsigned short*)(&y));

	long z1=(long)x1*(long)x1, c1=(long)y1*(long)y1;
	long z2=(long)x1*(long)x0, c2=(long)y1*(long)y0;
	long z3=(long)x0*(long)x0, c3=(long)y0*(long)y0;
	long z=z1+c1+(*((short*)(&z2)+1)+*((short*)(&c2)+1))*2;
	long z0=*((unsigned short*)(&z3)+1)+*((unsigned short*)(&c3)+1)+
			(*(unsigned short*)(&z2)+*(unsigned short*)(&c2))*2;

	z+=(z0>>16); z0&=(1<<16)-1;

	int e=0;
	unsigned long zz=z,s,t;

	if (zz!=0) {
		if (zz<(1<<16)) { e+=16; zz<<=16; }
		if (zz<(1<<24)) { e+=8; zz<<=8; }
		if (zz<(1<<28)) { e+=4; zz<<=4; }
		if (zz<(1<<30)) { e+=2; zz<<=2; }
		//ASSERT(e<=20);

		if (e<=20) {
			s=(zz>>20)-1024;
			t=(e>8)?(((zz<<12)+(unsigned long)(z0<<e-4))>>20):((zz<<12)>>20);
		}
		else {
			s=(zz>>20)-1024+((unsigned long)z0>>(36-e));
			t=(e<=24)?((unsigned long)z0>>(24-e)):((unsigned long)z0<<(e-24));
			t&=(1<<12)-1;
		}
	}
	else if (z0!=0) {
		zz=(unsigned)z0; e=32;
		if (zz<(1<<8)) { e+=8; zz<<=8; }
		if (zz<(1<<12)) { e+=4; zz<<=4; }
		if (zz<(1<<14)) { e+=2; zz<<=2; }
		s=(zz>>4)-1024; t=zz&((1<<4)-1);
	}
	else
		return abs(x)>abs(y)?abs(x):abs(y);

	return (unsigned long)((SqrtTable[s]+(((unsigned long)t*(unsigned long)(SqrtTable[s+1]-SqrtTable[s]))>>12)))>>(e/2);
}
#pragma warning(default:4554)

////////////////////////////////////////////////////////////

}