|
Server IP : 2a02:4780:11:1359:0:1d43:a566:2 / Your IP : 216.73.216.161 Web Server : LiteSpeed System : Linux in-mum-web1259.main-hosting.eu 4.18.0-553.37.1.lve.el8.x86_64 #1 SMP Mon Feb 10 22:45:17 UTC 2025 x86_64 User : u490972518 ( 490972518) PHP Version : 5.6.40 Disable Function : system, exec, shell_exec, passthru, mysql_list_dbs, ini_alter, dl, symlink, link, chgrp, leak, popen, apache_child_terminate, virtual, mb_send_mail MySQL : ON | cURL : ON | WGET : ON | Perl : OFF | Python : OFF Directory (0755) : /home/../opt/golang/1.22.0/src/maps/../math/ |
| [ Home ] | [ C0mmand ] | [ Upload File ] |
|---|
// Copyright 2010 The Go Authors. All rights reserved. // Use of this source code is governed by a BSD-style // license that can be found in the LICENSE file. #include "textflag.h" // func archHypot(p, q float64) float64 TEXT ·archHypot(SB),NOSPLIT,$0 // test bits for not-finite MOVL p_hi+4(FP), AX // high word p ANDL $0x7ff00000, AX CMPL AX, $0x7ff00000 JEQ not_finite MOVL q_hi+12(FP), AX // high word q ANDL $0x7ff00000, AX CMPL AX, $0x7ff00000 JEQ not_finite FMOVD p+0(FP), F0 // F0=p FABS // F0=|p| FMOVD q+8(FP), F0 // F0=q, F1=|p| FABS // F0=|q|, F1=|p| FUCOMI F0, F1 // compare F0 to F1 JCC 2(PC) // jump if F0 >= F1 FXCHD F0, F1 // F0=|p| (larger), F1=|q| (smaller) FTST // compare F0 to 0 FSTSW AX ANDW $0x4000, AX JNE 10(PC) // jump if F0 = 0 FXCHD F0, F1 // F0=q (smaller), F1=p (larger) FDIVD F1, F0 // F0=q(=q/p), F1=p FMULD F0, F0 // F0=q*q, F1=p FLD1 // F0=1, F1=q*q, F2=p FADDDP F0, F1 // F0=1+q*q, F1=p FSQRT // F0=sqrt(1+q*q), F1=p FMULDP F0, F1 // F0=p*sqrt(1+q*q) FMOVDP F0, ret+16(FP) RET FMOVDP F0, F1 // F0=0 FMOVDP F0, ret+16(FP) RET not_finite: // test bits for -Inf or +Inf MOVL p_hi+4(FP), AX // high word p ORL p_lo+0(FP), AX // low word p ANDL $0x7fffffff, AX CMPL AX, $0x7ff00000 JEQ is_inf MOVL q_hi+12(FP), AX // high word q ORL q_lo+8(FP), AX // low word q ANDL $0x7fffffff, AX CMPL AX, $0x7ff00000 JEQ is_inf MOVL $0x7ff80000, ret_hi+20(FP) // return NaN = 0x7FF8000000000001 MOVL $0x00000001, ret_lo+16(FP) RET is_inf: MOVL AX, ret_hi+20(FP) // return +Inf = 0x7FF0000000000000 MOVL $0x00000000, ret_lo+16(FP) RET