unsigned int julery_isqrt(unsigned int val) {
	unsigned int temp;
	unsigned int g;
	g=0;
	unsigned int b;
	b = 0x8000;
	unsigned int bshft;
	bshft = 15;
	int i;
	for (i = 0; i < 16; i++) {
		temp = g << 1;
		temp = temp + b;
		temp = temp << bshft;
		if ( val >= temp) {
			g += b;
			val -= temp;
		}
		bshft--;
		b = b >> 1;
	}
	return g;
}