VirtualBox

Ignore:
Timestamp:
Aug 5, 2022 12:22:06 PM (3 years ago)
Author:
vboxsync
svn:sync-xref-src-repo-rev:
152859
Message:

IPRT/nocrt: Make use of existing floating point assembly code, morphing it into double and float variants where those were missing. bugref:10261

Location:
trunk/src/VBox/Runtime/common/math
Files:
12 copied

Legend:

Unmodified
Added
Removed

  • trunk/src/VBox/Runtime/common/math/ceil.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT ceill - AMD64 & X86.
     3; IPRT - No-CRT ceil - AMD64 & X86.
    44;
    55
    … …  
    3232;;
    3333; Compute the smallest integral value not less than lrd.
    34 ; @returns st(0)
    35 ; @param    lrd     [rbp + 8]
    36 RT_NOCRT_BEGINPROC ceill
     34; @returns st(0) / xmm0
     35; @param    rd      [rbp + 8] / xmm0
     36RT_NOCRT_BEGINPROC ceil
    3737    push    xBP
    3838    mov     xBP, xSP
    3939    sub     xSP, 10h
    4040
    41     fld     tword [xBP + xCB*2]
     41%ifdef RT_ARCH_AMD64 ;; @todo there is probably some sse instruction for this.
     42    movsd   [xSP], xmm0
     43    fld     qword [xSP]
     44%else
     45    fld     qword [xBP + xCB*2]
     46%endif
    4247
    4348    ; Make it round up by modifying the fpu control word.
    … …  
    5560    fldcw   [xBP - 10h]
    5661
     62%ifdef RT_ARCH_AMD64
     63    fstp    qword [xSP]
     64    movsd   xmm0, [xSP]
     65%endif
    5766    leave
    5867    ret
    59 ENDPROC   RT_NOCRT(ceill)
     68ENDPROC   RT_NOCRT(ceil)
    6069

  • trunk/src/VBox/Runtime/common/math/ceilf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT ceill - AMD64 & X86.
     3; IPRT - No-CRT ceilf - AMD64 & X86.
    44;
    55
    … …  
    3131
    3232;;
    33 ; Compute the smallest integral value not less than lrd.
    34 ; @returns st(0)
    35 ; @param    lrd     [rbp + 8]
    36 RT_NOCRT_BEGINPROC ceill
     33; Compute the smallest integral value not less than r32.
     34; @returns st(0) / xmm0
     35; @param    r32      [rbp + 8] / xmm0
     36RT_NOCRT_BEGINPROC ceilf
    3737    push    xBP
    3838    mov     xBP, xSP
    3939    sub     xSP, 10h
    4040
    41     fld     tword [xBP + xCB*2]
     41%ifdef RT_ARCH_AMD64 ;; @todo there is probably some sse instruction for this.
     42    movss   [xSP], xmm0
     43    fld     dword [xSP]
     44%else
     45    fld     dword [xBP + xCB*2]
     46%endif
    4247
    4348    ; Make it round up by modifying the fpu control word.
    … …  
    5560    fldcw   [xBP - 10h]
    5661
     62%ifdef RT_ARCH_AMD64
     63    fstp    dword [xSP]
     64    movss   xmm0, [xSP]
     65%endif
    5766    leave
    5867    ret
    59 ENDPROC   RT_NOCRT(ceill)
     68ENDPROC   RT_NOCRT(ceilf)
    6069

  • trunk/src/VBox/Runtime/common/math/cos.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT cosl - AMD64 & X86.
     3; IPRT - No-CRT cos - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; compute the cosine of ldr, measured in radians.
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 RT_NOCRT_BEGINPROC cosl
     32; compute the cosine of dr, measured in radians.
     33; @returns st(0) / xmm0
     34; @param    rd      [rbp + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC cos
    3636    push    xBP
    3737    mov     xBP, xSP
     38%ifdef RT_ARCH_AMD64
    3839    sub     xSP, 10h
    3940
    40     fld     tword [xBP + xCB*2]
     41    movsd   [xSP], xmm0
     42    fld     qword [xSP]
     43%else
     44    fld     qword [xBP + xCB*2]
     45%endif
    4146    fcos
    4247    fnstsw  ax
    … …  
    5762
    5863.done:
     64%ifdef RT_ARCH_AMD64
     65    fstp    qword [xSP]
     66    movsd   xmm0, [xSP]
     67%endif
    5968    leave
    6069    ret
    61 ENDPROC   RT_NOCRT(cosl)
     70ENDPROC   RT_NOCRT(cos)
    6271

  • trunk/src/VBox/Runtime/common/math/cosf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT cosl - AMD64 & X86.
     3; IPRT - No-CRT cosf - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; compute the cosine of ldr, measured in radians.
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 RT_NOCRT_BEGINPROC cosl
     32; compute the cosine of r32, measured in radians.
     33; @returns st(0) / xmm0
     34; @param    r32     [rbp + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC cosf
    3636    push    xBP
    3737    mov     xBP, xSP
     38%ifdef RT_ARCH_AMD64
    3839    sub     xSP, 10h
    3940
    40     fld     tword [xBP + xCB*2]
     41    movss   [xSP], xmm0
     42    fld     dword [xSP]
     43%else
     44    fld     dword [xBP + xCB*2]
     45%endif
    4146    fcos
    4247    fnstsw  ax
    … …  
    5762
    5863.done:
     64%ifdef RT_ARCH_AMD64
     65    fstp    dword [xSP]
     66    movss   xmm0, [xSP]
     67%endif
    5968    leave
    6069    ret
    61 ENDPROC   RT_NOCRT(cosl)
     70ENDPROC   RT_NOCRT(cosf)
    6271

  • trunk/src/VBox/Runtime/common/math/ldexp.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT ldexpl - AMD64 & X86.
     3; IPRT - No-CRT ldexp - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Computes lrd * 2^exp
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 ; @param    exp     [ebp + 14h]  gcc:edi  msc:ecx
    36 RT_NOCRT_BEGINPROC ldexpl
     32; Computes rd * 2^exp
     33; @returns st(0) / xmm0
     34; @param    rd      [rbp + xCB*2] / xmm0
     35; @param    exp     [ebp + 10h]  gcc:edi  msc:ecx
     36RT_NOCRT_BEGINPROC ldexp
    3737    push    xBP
    3838    mov     xBP, xSP
     39%ifdef RT_ARCH_AMD64
     40    sub     xSP, 20h
     41%else
    3942    sub     xSP, 10h
     43%endif
    4044
    4145    ; load exp
    42 %ifdef RT_ARCH_AMD64 ; ASSUMES ONLY GCC HERE!
     46%ifdef RT_ARCH_AMD64
     47 %ifdef ASM_CALL64_GCC
    4348    mov     [rsp], edi
     49 %else
     50    mov     [rsp], ecx
     51 %endif
    4452    fild    dword [rsp]
     53    movsd   [rsp + 8], xmm0
     54    fld     qword [rsp + 8]
    4555%else
    46     fild    dword [ebp + xCB*2 + RTLRD_CB]
     56    fild    dword [ebp + xCB*2 + 8]
     57    fld     qword [xBP + xCB*2]
    4758%endif
    48     fld     tword [xBP + xCB*2]
    4959    fscale
    5060    fstp    st1
    5161
     62%ifdef RT_ARCH_AMD64
     63    fstp    qword [xSP]
     64    movsd   xmm0, [xSP]
     65%endif
    5266    leave
    5367    ret
    54 ENDPROC   RT_NOCRT(ldexpl)
     68ENDPROC   RT_NOCRT(ldexp)
    5569

  • trunk/src/VBox/Runtime/common/math/ldexpf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT ldexpl - AMD64 & X86.
     3; IPRT - No-CRT ldexpf - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Computes lrd * 2^exp
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 ; @param    exp     [ebp + 14h]  gcc:edi  msc:ecx
    36 RT_NOCRT_BEGINPROC ldexpl
     32; Computes r32 * 2^exp
     33; @returns st(0) / xmm0
     34; @param    r32     [rbp + xCB*2] / xmm0
     35; @param    exp     [ebp + 10h]  gcc:edi  msc:ecx
     36RT_NOCRT_BEGINPROC ldexpf
    3737    push    xBP
    3838    mov     xBP, xSP
     39%ifdef RT_ARCH_AMD64
     40    sub     xSP, 20h
     41%else
    3942    sub     xSP, 10h
     43%endif
    4044
    4145    ; load exp
    42 %ifdef RT_ARCH_AMD64 ; ASSUMES ONLY GCC HERE!
     46%ifdef RT_ARCH_AMD64
     47 %ifdef ASM_CALL64_GCC
    4348    mov     [rsp], edi
     49 %else
     50    mov     [rsp], ecx
     51 %endif
    4452    fild    dword [rsp]
     53    movss   [rsp + 8], xmm0
     54    fld     dword [rsp + 8]
    4555%else
    46     fild    dword [ebp + xCB*2 + RTLRD_CB]
     56    fild    dword [ebp + xCB*2 + 8]
     57    fld     dword [xBP + xCB*2]
    4758%endif
    48     fld     tword [xBP + xCB*2]
    4959    fscale
    5060    fstp    st1
    5161
     62%ifdef RT_ARCH_AMD64
     63    fstp    dword [xSP]
     64    movss   xmm0, [xSP]
     65%endif
    5266    leave
    5367    ret
    54 ENDPROC   RT_NOCRT(ldexpl)
     68ENDPROC   RT_NOCRT(ldexpf)
    5569

  • trunk/src/VBox/Runtime/common/math/log.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT logl - AMD64 & X86.
     3; IPRT - No-CRT log - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; compute the natural logarithm of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 RT_NOCRT_BEGINPROC logl
     32; compute the natural logarithm of rd
     33; @returns st(0) / xmm0
     34; @param    rd      [rbp + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC log
    3636    push    xBP
    3737    mov     xBP, xSP
     38%ifdef RT_ARCH_AMD64
    3839    sub     xSP, 10h
     40%endif
    3941
    4042    fldln2                              ; st0=log(2)
    41     fld     tword [xBP + xCB*2]         ; st1=log(2) st0=lrd
     43%ifdef RT_ARCH_AMD64
     44    movsd   [xSP], xmm0
     45    fld     qword [xSP]
     46%else
     47    fld     qword [xBP + xCB*2]         ; st1=log(2) st0=lrd
     48%endif
    4249    fld     st0                         ; st1=log(2) st0=lrd st0=lrd
    4350    fsub    qword [.one xWrtRIP]        ; st2=log(2) st1=lrd st0=lrd-1.0
    … …  
    5865
    5966.done:
     67%ifdef RT_ARCH_AMD64
     68    fstp    qword [xSP]
     69    movsd   xmm0, [xSP]
     70%endif
    6071    leave
    6172    ret
    6273.one:   dq  1.0
    6374.limit: dq  0.29
    64 ENDPROC   RT_NOCRT(logl)
     75ENDPROC   RT_NOCRT(log)
    6576

  • trunk/src/VBox/Runtime/common/math/logf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT logl - AMD64 & X86.
     3; IPRT - No-CRT logf - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; compute the natural logarithm of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [rbp + xCB*2]
    35 RT_NOCRT_BEGINPROC logl
     32; compute the natural logarithm of r32
     33; @returns st(0) / xmm0
     34; @param    r32     [rbp + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC logf
    3636    push    xBP
    3737    mov     xBP, xSP
     38%ifdef RT_ARCH_AMD64
    3839    sub     xSP, 10h
     40%endif
    3941
    4042    fldln2                              ; st0=log(2)
    41     fld     tword [xBP + xCB*2]         ; st1=log(2) st0=lrd
     43%ifdef RT_ARCH_AMD64
     44    movss   [xSP], xmm0
     45    fld     dword [xSP]
     46%else
     47    fld     dword [xBP + xCB*2]         ; st1=log(2) st0=lrd
     48%endif
    4249    fld     st0                         ; st1=log(2) st0=lrd st0=lrd
    4350    fsub    qword [.one xWrtRIP]        ; st2=log(2) st1=lrd st0=lrd-1.0
    … …  
    5865
    5966.done:
     67%ifdef RT_ARCH_AMD64
     68    fstp    dword [xSP]
     69    movss   xmm0, [xSP]
     70%endif
    6071    leave
    6172    ret
    6273.one:   dq  1.0
    6374.limit: dq  0.29
    64 ENDPROC   RT_NOCRT(logl)
     75ENDPROC   RT_NOCRT(logf)
    6576

  • trunk/src/VBox/Runtime/common/math/sin.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT sinl - AMD64 & X86.
     3; IPRT - No-CRT sin - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Compute the sine of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [xSP + xCB*2]
    35 RT_NOCRT_BEGINPROC sinl
     32; Compute the sine of rd
     33; @returns st(0)/xmm0
     34; @param    rd      [xSP + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC sin
    3636    push    xBP
    3737    mov     xBP, xSP
     38
     39%ifdef RT_ARCH_AMD64
    3840    sub     xSP, 10h
    3941
    40     fld     tword [xBP + xCB*2]
     42    movsd   [xSP], xmm0
     43    fld     qword [xSP]
     44%else
     45    fld     qword [xBP + xCB*2]
     46%endif
    4147    fsin
    4248    fnstsw  ax
    … …  
    5662
    5763.done:
     64%ifdef RT_ARCH_AMD64
     65    fstp    qword [xSP]
     66    movsd   xmm0, [xSP]
     67%endif
    5868    leave
    5969    ret
    60 ENDPROC   RT_NOCRT(sinl)
     70ENDPROC   RT_NOCRT(sin)
    6171

  • trunk/src/VBox/Runtime/common/math/sinf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT sinl - AMD64 & X86.
     3; IPRT - No-CRT sinf - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Compute the sine of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [xSP + xCB*2]
    35 RT_NOCRT_BEGINPROC sinl
     32; Compute the sine of r32
     33; @returns st(0)/xmm0
     34; @param    r32     [xSP + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC sinf
    3636    push    xBP
    3737    mov     xBP, xSP
     38
     39%ifdef RT_ARCH_AMD64
    3840    sub     xSP, 10h
    3941
    40     fld     tword [xBP + xCB*2]
     42    movss   [xSP], xmm0
     43    fld     dword [xSP]
     44%else
     45    fld     dword [xBP + xCB*2]
     46%endif
    4147    fsin
    4248    fnstsw  ax
    … …  
    5662
    5763.done:
     64%ifdef RT_ARCH_AMD64
     65    fstp    dword [xSP]
     66    movss   xmm0, [xSP]
     67%endif
    5868    leave
    5969    ret
    60 ENDPROC   RT_NOCRT(sinl)
     70ENDPROC   RT_NOCRT(sinf)
    6171

  • trunk/src/VBox/Runtime/common/math/tan.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT tanl - AMD64 & X86.
     3; IPRT - No-CRT tan - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Compute the sine of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [xSP + xCB*2]
    35 RT_NOCRT_BEGINPROC tanl
     32; Compute the sine of rd
     33; @returns st(0) / xmm0
     34; @param    rd      [xSP + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC tan
    3636    push    xBP
    3737    mov     xBP, xSP
     38
     39%ifdef RT_ARCH_AMD64
    3840    sub     xSP, 10h
    3941
    40     fld     tword [xBP + xCB*2]
     42    movsd   [xSP], xmm0
     43    fld     qword [xSP]
     44%else
     45    fld     qword [xBP + xCB*2]
     46%endif
    4147    fptan
    4248    fnstsw  ax
    … …  
    5763.done:
    5864    fstp    st0
     65%ifdef RT_ARCH_AMD64
     66    fstp    qword [xSP]
     67    movsd   xmm0, [xSP]
     68%endif
    5969    leave
    6070    ret
    61 ENDPROC   RT_NOCRT(tanl)
     71ENDPROC   RT_NOCRT(tan)
    6272

  • trunk/src/VBox/Runtime/common/math/tanf.asm

    r96039 r96060  
    11; $Id$
    22;; @file
    3 ; IPRT - No-CRT tanl - AMD64 & X86.
     3; IPRT - No-CRT tanf - AMD64 & X86.
    44;
    55
    … …  
    3030
    3131;;
    32 ; Compute the sine of lrd
    33 ; @returns st(0)
    34 ; @param    lrd     [xSP + xCB*2]
    35 RT_NOCRT_BEGINPROC tanl
     32; Compute the sine of r32
     33; @returns st(0) / xmm0
     34; @param    r32     [xSP + xCB*2] / xmm0
     35RT_NOCRT_BEGINPROC tanf
    3636    push    xBP
    3737    mov     xBP, xSP
     38
     39%ifdef RT_ARCH_AMD64
    3840    sub     xSP, 10h
    3941
    40     fld     tword [xBP + xCB*2]
     42    movss   [xSP], xmm0
     43    fld     dword [xSP]
     44%else
     45    fld     dword [xBP + xCB*2]
     46%endif
    4147    fptan
    4248    fnstsw  ax
    … …  
    5763.done:
    5864    fstp    st0
     65%ifdef RT_ARCH_AMD64
     66    fstp    dword [xSP]
     67    movss   xmm0, [xSP]
     68%endif
    5969    leave
    6070    ret
    61 ENDPROC   RT_NOCRT(tanl)
     71ENDPROC   RT_NOCRT(tanf)
    6272
Note: See TracChangeset for help on using the changeset viewer.

© 2024 Oracle Support Privacy / Do Not Sell My Info Terms of Use Trademark Policy Automated Access Etiquette