From 1756c4a269c72d48977e05dabc3091c1190fc24c Mon Sep 17 00:00:00 2001 From: Alex Ionescu Date: Mon, 28 Nov 2005 22:24:37 +0000 Subject: [PATCH] - Remove inlined gas assembly + c code mix in math.c and write the functions in real assembly code. - Fix pch usage and qsort callback definition. RTL can now be built with msvc. svn path=/trunk/; revision=19729 --- reactos/lib/rtl/i386/math_asm.S | 223 ++++++++++++++++++++++++++++++-- reactos/lib/rtl/largeint.c | 1 - reactos/lib/rtl/math.c | 215 ------------------------------ reactos/lib/rtl/qsort.c | 5 +- reactos/lib/rtl/rtl.h | 3 +- reactos/lib/rtl/rtl.xml | 1 - 6 files changed, 218 insertions(+), 230 deletions(-) delete mode 100644 reactos/lib/rtl/math.c diff --git a/reactos/lib/rtl/i386/math_asm.S b/reactos/lib/rtl/i386/math_asm.S index ca1d6caf6eb..3ec7c072409 100644 --- a/reactos/lib/rtl/i386/math_asm.S +++ b/reactos/lib/rtl/i386/math_asm.S @@ -28,7 +28,7 @@ * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS - * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) + * OR SERVICES// LOSS OF USE, DATA, OR PROFITS// OR BUSINESS INTERRUPTION) * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF @@ -48,12 +48,32 @@ .globl __alldiv .globl __aulldvrm .globl __alldvrm +.globl _atan +.globl _ceil +.globl _cos +.globl _fabs +.globl _floor +.globl _log +.globl _pow +.globl _sin +.globl _sqrt +.globl _tan +.globl __fltused + +/* DATA ********************************************************************/ + +fzero: + .long 0 // Floating point zero + .long 0 // Floating point zero + +__fltused: + .long 0x9875 /* FUNCTIONS ***************************************************************/ /* * long long - * __alldiv(long long Dividend, long long Divisor); + * __alldiv(long long Dividend, long long Divisor)// * * Parameters: * [ESP+04h] - long long Dividend @@ -71,7 +91,7 @@ __alldiv: /* * long long - * __allmul(long long Multiplier, long long Multiplicand); + * __allmul(long long Multiplier, long long Multiplicand)// * * Parameters: * [ESP+04h] - long long Multiplier @@ -114,7 +134,7 @@ __allmul: /* * unsigned long long - * __aullrem(unsigned long long Dividend, unsigned long long Divisor); + * __aullrem(unsigned long long Dividend, unsigned long long Divisor)// * * Parameters: * [ESP+04h] - unsigned long long Dividend @@ -132,7 +152,7 @@ __aullrem: /* * long long - * __allshl(long long Value, unsigned char Shift); + * __allshl(long long Value, unsigned char Shift)// * * Parameters: * EDX:EAX - signed long long value to be shifted left @@ -154,7 +174,7 @@ __allshl: /* * long long - * __allshr(long long Value, unsigned char Shift); + * __allshr(long long Value, unsigned char Shift)// * * Parameters: * EDX:EAX - signed long long value to be shifted right @@ -176,7 +196,7 @@ __allshr: /* * unsigned long long - * __aulldiv(unsigned long long Dividend, unsigned long long Divisor); + * __aulldiv(unsigned long long Dividend, unsigned long long Divisor)// * * Parameters: * [ESP+04h] - unsigned long long Dividend @@ -194,7 +214,7 @@ __aulldiv: /* * unsigned long long - * __aullshr(unsigned long long Value, unsigned char Shift); + * __aullshr(unsigned long long Value, unsigned char Shift)// * * Parameters: * EDX:EAX - unsigned long long value to be shifted right @@ -215,7 +235,7 @@ __aullshr: /* * long long - * __allrem(long long Dividend, long long Divisor); + * __allrem(long long Dividend, long long Divisor)// * * Parameters: * [ESP+04h] - long long Dividend @@ -629,3 +649,188 @@ __aulldvrm: ret 16 +_atan: + push ebp + mov ebp,esp + fld qword ptr [ebp+8] // Load real from stack + fld1 // Load constant 1 + fpatan // Take the arctangent + pop ebp + ret + +_ceil: + push ebp + mov ebp,esp + sub esp,4 // Allocate temporary space + fld qword ptr [ebp+8] // Load real from stack + fstcw [ebp-2] // Save control word + fclex // Clear exceptions + mov word ptr [ebp-4],0xb63 // Rounding control word + fldcw [ebp-4] // Set new rounding control + frndint // Round to integer + fclex // Clear exceptions + fldcw [ebp-2] // Restore control word + mov esp,ebp // Deallocate temporary space + pop ebp + ret + +_cos: + push ebp + mov ebp,esp // Point to the stack frame + fld qword ptr [ebp+8] // Load real from stack + fcos // Take the cosine + pop ebp + ret + +_fabs: + push ebp + mov ebp,esp + fld qword ptr [ebp+8] // Load real from stack + fabs // Take the absolute value + pop ebp + ret + +_floor: + push ebp + mov ebp,esp + sub esp,4 // Allocate temporary space + fld qword ptr [ebp+8] // Load real from stack + fstcw [ebp-2] // Save control word + fclex // Clear exceptions + mov word ptr [ebp-4],0x763 // Rounding control word + fldcw [ebp-4] // Set new rounding control + frndint // Round to integer + fclex // Clear exceptions + fldcw [ebp-2] // Restore control word + mov esp,ebp + pop ebp + ret + +_log: + push ebp + mov ebp,esp + fld qword ptr [ebp+8] // Load real from stack + fldln2 // Load log base e of 2 + fxch st(1) // Exchange st, st(1) + fyl2x // Compute the natural log(x) + pop ebp + ret + +_pow: + push ebp + mov ebp,esp + sub esp,12 // Allocate temporary space + push edi // Save register edi + push eax // Save register eax + mov dword ptr [ebp-12],0 // Set negation flag to zero + fld qword ptr [ebp+16] // Load real from stack + fld qword ptr [ebp+8] // Load real from stack + mov edi,offset flat:fzero // Point to real zero + fcom qword ptr [edi] // Compare x with zero + fstsw ax // Get the FPU status word + mov al,ah // Move condition flags to AL + lahf // Load Flags into AH + and al, 0b01000101 // Isolate C0, C2 and C3 + and ah,not 0b01000101 // Turn off CF, PF and ZF + or ah,al // Set new CF, PF and ZF + sahf // Store AH into Flags + jb __fpow1 // Re-direct if x < 0 + ja __fpow3 // Re-direct if x > 0 + fxch // Swap st, st(1) + fcom qword ptr [edi] // Compare y with zero + fxch // Restore x as top of stack + fstsw ax // Get the FPU status word + mov al,ah // Move condition flags to AL + lahf // Load Flags into AH + and al, 0b01000101 // Isolate C0, C2 and C3 + and ah,not 0b01000101 // Turn off CF, PF and ZF + or ah,al // Set new CF, PF and ZF + sahf // Store AH into Flags + ja __fpow3 // Re-direct if y > 0 + fstp st(1) // Set new stack top and pop + mov eax,1 // Set domain error (EDOM) + jmp __fpow5 // End of case +__fpow1: fxch // Put y on top of stack + fld st // Duplicate y as st(1) + frndint // Round to integer + fxch // Put y on top of stack + fcomp // y = int(y) ? + fstsw ax // Get the FPU status word + mov al,ah // Move condition flags to AL + lahf // Load Flags into AH + and al, 0b01000101 // Isolate C0, C2 and C3 + and ah,not 0b01000101 // Turn off CF, PF and ZF + or ah,al // Set new CF, PF and ZF + sahf // Store AH into Flags + je __fpow2 // Proceed if y = int(y) + fstp st(1) // Set new stack top and pop + fldz // Set result to zero + fstp st(1) // Set new stack top and pop + mov eax,1 // Set domain error (EDOM) + jmp __fpow5 // End of case +__fpow2: fist dword ptr [ebp-12] // Store y as integer + and dword ptr [ebp-12],1 // Set bit if y is odd + fxch // Put x on top of stack + fabs // x = |x| +__fpow3: fldln2 // Load log base e of 2 + fxch st(1) // Exchange st, st(1) + fyl2x // Compute the natural log(x) + fmulp // Compute y * ln(x) + fldl2e // Load log base 2(e) + fmulp st(1),st // Multiply x * log base 2(e) + fst st(1) // Push result + frndint // Round to integer + fsub st(1),st // Subtract + fxch // Exchange st, st(1) + f2xm1 // Compute 2 to the (x - 1) + fld1 // Load real number 1 + faddp // 2 to the x + fscale // Scale by power of 2 + fstp st(1) // Set new stack top and pop + test dword ptr [ebp-12],1 // Negation required ? + jz __fpow4 // No, re-direct + fchs // Negate the result +__fpow4: fstp qword ptr [ebp-8] // Save (double)pow(x, y) + fld qword ptr [ebp-8] // Load (double)pow(x, y) + fxam // Examine st + fstsw ax // Get the FPU status word + cmp ah,5 // Infinity ? + jne __fpow6 // No, end of case + mov eax,2 // Set range error (ERANGE) + // Get errno pointer offset +__fpow5: int 3 + mov edi,0 // TODO: offset flat:__crt_errno + mov edi,[edi] // Get C errno variable pointer + mov dword ptr [edi],eax // Set errno +__fpow6: pop eax // Restore register eax + pop edi // Restore register edi + mov esp,ebp // Deallocate temporary space + pop ebp + ret + +_sin: + push ebp // Save register bp + mov ebp,esp // Point to the stack frame + fld qword ptr [ebp+8] // Load real from stack + fsin // Take the sine + pop ebp // Restore register bp + ret + +_sqrt: + push ebp + mov ebp,esp + fld qword ptr [ebp+8] // Load real from stack + fsqrt // Take the square root + pop ebp + ret + +_tan: + push ebp + mov ebp,esp + sub esp,4 // Allocate temporary space + fld qword ptr [ebp+8] // Load real from stack + fptan // Take the tangent + fstp dword ptr [ebp-4] // Throw away the constant 1 + mov esp,ebp // Deallocate temporary space + pop ebp + ret diff --git a/reactos/lib/rtl/largeint.c b/reactos/lib/rtl/largeint.c index af9f3061048..e5ef9125c02 100644 --- a/reactos/lib/rtl/largeint.c +++ b/reactos/lib/rtl/largeint.c @@ -7,7 +7,6 @@ /* INCLUDES *****************************************************************/ -#define NO_RTL_INLINES #include #define NDEBUG diff --git a/reactos/lib/rtl/math.c b/reactos/lib/rtl/math.c deleted file mode 100644 index f292e5c4e9b..00000000000 --- a/reactos/lib/rtl/math.c +++ /dev/null @@ -1,215 +0,0 @@ -/* COPYRIGHT: See COPYING in the top level directory - * Copyright (C) 1995, 1996, 1997 Free Software Foundation, Inc. - * PROJECT: ReactOS system libraries - * FILE: lib/rtl/math.c - * PURPOSE: Math functions for i387. - * PROGRAMMER: John C. Bowman - */ - -/* INCLUDES *****************************************************************/ - -#include - -#define NDEBUG -#include - -/* FUNCTIONS *****************************************************************/ - -double atan (double __x); -double ceil (double __x); -double cos (double __x); -double fabs (double __x); -double floor (double __x); -double log (double __x); -double __log2 (double __x); -double pow (double __x, double __y); -double sin (double __x); -double sqrt (double __x); -double tan (double __x); - -int _fltused = 0x9875; - -double atan (double __x) -{ - register double __val; - __asm __volatile__ - ("fld1\n\t" - "fpatan" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double ceil (double __x) -{ - register double __val; - __volatile unsigned short int __cw, __cwtmp; - - __asm __volatile ("fnstcw %0" : "=m" (__cw)); - __cwtmp = (__cw & 0xf3ff) | 0x0800; /* rounding up */ - __asm __volatile ("fldcw %0" : : "m" (__cwtmp)); - __asm __volatile ("frndint" : "=t" (__val) : "0" (__x)); - __asm __volatile ("fldcw %0" : : "m" (__cw)); - - return __val; -} - -/* - * @implemented - */ -double cos (double __x) -{ - register double __val; - __asm __volatile__ - ("fcos" - : "=t" (__val): "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double fabs (double __x) -{ - register double __val; - __asm __volatile__ - ("fabs" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double floor (double __x) -{ - register double __val; - __volatile unsigned short int __cw, __cwtmp; - - __asm __volatile ("fnstcw %0" : "=m" (__cw)); - __cwtmp = (__cw & 0xf3ff) | 0x0400; /* rounding down */ - __asm __volatile ("fldcw %0" : : "m" (__cwtmp)); - __asm __volatile ("frndint" : "=t" (__val) : "0" (__x)); - __asm __volatile ("fldcw %0" : : "m" (__cw)); - - return __val; -} - -/* - * @implemented - */ -double log (double __x) -{ - register double __val; - __asm __volatile__ - ("fldln2\n\t" - "fxch\n\t" - "fyl2x" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -double __log2 (double __x) -{ - register double __val; - __asm __volatile__ - ("fld1\n\t" - "fxch\n\t" - "fyl2x" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double pow (double __x, double __y) -{ - register double __val, __exponent; - long __p = (long) __y; - - if (__x == 0.0 && __y > 0.0) - return 0.0; - if (__y == (double) __p) - { - double __r = 1.0; - if (__p == 0) - return 1.0; - if (__p < 0) - { - __p = -__p; - __x = 1.0 / __x; - } - while (1) - { - if (__p & 1) - __r *= __x; - __p >>= 1; - if (__p == 0) - return __r; - __x *= __x; - } - /* NOTREACHED */ - } - __asm __volatile__ - ("fmul %%st(1) # y * log2(x)\n\t" - "fst %%st(1)\n\t" - "frndint # int(y * log2(x))\n\t" - "fxch\n\t" - "fsub %%st(1) # fract(y * log2(x))\n\t" - "f2xm1 # 2^(fract(y * log2(x))) - 1\n\t" - : "=t" (__val), "=u" (__exponent) : "0" (__log2 (__x)), "1" (__y)); - __val += 1.0; - __asm __volatile__ - ("fscale" - : "=t" (__val) : "0" (__val), "u" (__exponent)); - - return __val; -} - -/* - * @implemented - */ -double sin (double __x) -{ - register double __val; - __asm __volatile__ - ("fsin" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double sqrt (double __x) -{ - register double __val; - __asm __volatile__ - ("fsqrt" - : "=t" (__val) : "0" (__x)); - - return __val; -} - -/* - * @implemented - */ -double tan (double __x) -{ - register double __val; - register double __val2 __attribute__ ((unused)); - __asm __volatile__ - ("fptan" - : "=t" (__val2), "=u" (__val) : "0" (__x)); - - return __val; -} diff --git a/reactos/lib/rtl/qsort.c b/reactos/lib/rtl/qsort.c index c9951420c52..d3fdc1321fe 100644 --- a/reactos/lib/rtl/qsort.c +++ b/reactos/lib/rtl/qsort.c @@ -37,8 +37,7 @@ * FILE: lib/rtl/qsort.c * PROGRAMMER: Adapted from CygWin newlib 2000-03-12. */ - -#include +#include #ifndef __GNUC__ #define inline @@ -46,7 +45,7 @@ /* FIXME: these types should be from the default includes */ -typedef int (* _pfunccmp_t) (char *, char *); +typedef int (* _pfunccmp_t) (const void *, const void *); #define min(a,b) ((a)<(b)?(a):(b)) diff --git a/reactos/lib/rtl/rtl.h b/reactos/lib/rtl/rtl.h index e80b3483baf..78b9390f06d 100644 --- a/reactos/lib/rtl/rtl.h +++ b/reactos/lib/rtl/rtl.h @@ -12,8 +12,9 @@ #define _NTSYSTEM_ #define _NTDLLBUILD_ #define WIN32_NO_STATUS -#define NTOS_MODE_USER #define _INC_SWPRINTF_INL_ +#define _CRT_SECURE_NO_DEPRECATE +#define NO_RTL_INLINES /* C Headers */ #include diff --git a/reactos/lib/rtl/rtl.xml b/reactos/lib/rtl/rtl.xml index 5b797755c21..50ecc6a6c52 100644 --- a/reactos/lib/rtl/rtl.xml +++ b/reactos/lib/rtl/rtl.xml @@ -34,7 +34,6 @@ handle.c heap.c image.c - math.c mem.c message.c largeint.c