Intravision/reactos
1
0
Fork 0
mirror of https://github.com/reactos/reactos.git synced 2025-05-16 15:50:24 +00:00
Code Issues Projects Releases Packages Wiki Activity

- Remove inlined gas assembly + c code mix in math.c and write the functions in real assembly code.

Browse source 
- Fix pch usage and qsort callback definition.
RTL can now be built with msvc.

svn path=/trunk/; revision=19729
This commit is contained in:
Alex Ionescu 2005-11-28 22:24:37 +00:00
parent 0586868b52
commit 1756c4a269
6 changed files with 218 additions and 230 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

223
reactos/lib/rtl/i386/math_asm.S
View file

@ -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

1
reactos/lib/rtl/largeint.c
View file

@ -7,7 +7,6 @@
/* INCLUDES *****************************************************************/
#define NO_RTL_INLINES
#include <rtl.h>
#define NDEBUG

215
reactos/lib/rtl/math.c
View file

@ -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 <bowman@ipp-garching.mpg.de>
*/
/* INCLUDES *****************************************************************/
#include <rtl.h>
#define NDEBUG
#include <debug.h>
/* 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;
}

5
reactos/lib/rtl/qsort.c
View file

@ -37,8 +37,7 @@
* FILE: lib/rtl/qsort.c
* PROGRAMMER: Adapted from CygWin newlib 2000-03-12.
*/
#include <stddef.h>
#include <rtl.h>
#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))

3
reactos/lib/rtl/rtl.h
View file

@ -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 <stdio.h>

1
reactos/lib/rtl/rtl.xml
View file

@ -34,7 +34,6 @@
<file>handle.c</file>
<file>heap.c</file>
<file>image.c</file>
<file>math.c</file>
<file>mem.c</file>
<file>message.c</file>
<file>largeint.c</file>