X-Git-Url: https://git.ucc.asn.au/?a=blobdiff_plain;f=Usermode%2FLibraries%2Flibgcc.so_src%2Flibgcc.c;h=8f6021f20e892677e31d65f926b4f227f552a1f8;hb=25f7e9ab0f31ca486c0c981a406d381e160637a4;hp=010cc35a9ea2cd5370dc87efc61ccdb0e2393777;hpb=a2210987109ab5a6337c72b45f7e52cfc9092f8f;p=tpg%2Facess2.git

diff --git a/Usermode/Libraries/libgcc.so_src/libgcc.c b/Usermode/Libraries/libgcc.so_src/libgcc.c
index 010cc35a..8f6021f2 100644
--- a/Usermode/Libraries/libgcc.so_src/libgcc.c
+++ b/Usermode/Libraries/libgcc.so_src/libgcc.c
@@ -25,14 +25,66 @@ void __stack_chk_fail()
  */
 uint64_t __udivdi3(uint64_t Num, uint64_t Den)
 {
+	#if 0
 	uint64_t	ret = 0;
 	if(Den == 0)	// Call Div by Zero Error
 		__asm__ __volatile__ ("int $0");
+	
+	if(Den == 1)	return Num;
+	if(Den == 2)	return Num >> 1;
+	if(Den == 4)	return Num >> 2;
+	if(Den == 8)	return Num >> 3;
+	if(Den == 16)	return Num >> 4;
+	if(Den == 32)	return Num >> 5;
+	if(Den == 64)	return Num >> 6;
+	if(Den == 128)	return Num >> 7;
+	if(Den == 256)	return Num >> 8;
+	
 	while(Num > Den) {
 		ret ++;
 		Num -= Den;
 	}
 	return ret;
+	#else
+	uint64_t	P[2];
+	uint64_t	q;
+	 int	i;
+	
+	if(Den == 0)	__asm__ __volatile__ ("int $0x0");
+	// Common speedups
+	if(Num <= 0xFFFFFFFF && Den <= 0xFFFFFFFF)
+		return (uint32_t)Num / (uint32_t)Den;
+	if(Den == 1)	return Num;
+	if(Den == 2)	return Num >> 1;
+	if(Den == 16)	return Num >> 4;
+	if(Num < Den)	return 0;
+	if(Num < Den*2)	return 1;
+	if(Num == Den*2)	return 2;
+	
+	// Restoring division, from wikipedia
+	// http://en.wikipedia.org/wiki/Division_(digital)
+	P[0] = Num;	P[1] = 0;
+	for( i = 64; i--; )
+	{
+		// P <<= 1;
+		P[1] = (P[1] << 1) | (P[0] >> 63);
+		P[0] = P[0] << 1;
+		
+		// P -= Den << 64
+		P[1] -= Den;
+		
+		// P >= 0
+		if( !(P[1] & (1ULL<<63)) ) {
+			q |= (uint64_t)1 << (63-i);
+		}
+		else {
+			//q |= 0 << (63-i);
+			P[1] += Den;
+		}
+	}
+	
+	return q;
+	#endif
 }
 
 /**
@@ -41,9 +93,36 @@ uint64_t __udivdi3(uint64_t Num, uint64_t Den)
  */
 uint64_t __umoddi3(uint64_t Num, uint64_t Den)
 {
-	if(Den == 0)	// Call Div by Zero Error
-		__asm__ __volatile__ ("int $0");
-	while(Num > Den)
-		Num -= Den;
+	#if 0
+	if(Den == 0)	__asm__ __volatile__ ("int $0");	// Call Div by Zero Error
+	
+	if(Den == 1)	return 0;
+	if(Den == 2)	return Num & 0x01;
+	if(Den == 4)	return Num & 0x03;
+	if(Den == 8)	return Num & 0x07;
+	if(Den == 16)	return Num & 0x0F;
+	if(Den == 32)	return Num & 0x1F;
+	if(Den == 64)	return Num & 0x3F;
+	if(Den == 128)	return Num & 0x3F;
+	if(Den == 256)	return Num & 0x7F;
+	
+	while(Num >= Den)	Num -= Den;
+	
 	return Num;
+	#else
+	if(Den == 0)	__asm__ __volatile__ ("int $0");	// Call Div by Zero Error
+	
+	// Speedups
+	if(Num < Den)	return Num;
+	if(Num == Den)	return 0;
+	if(Num <= 0xFFFFFFFF && Den <= 0xFFFFFFFF)
+		return (uint32_t)Num % (uint32_t)Den;
+	
+	// Speedups for common operations
+	if(Den == 1)	return 0;
+	if(Den == 2)	return Num & 0x01;
+	if(Den == 8)	return Num & 0x07;
+	if(Den == 16)	return Num & 0x0F;
+	return Num - __udivdi3(Num, Den) * Den;
+	#endif
 }