Back to index

glibc  2.9
Classes | Defines | Functions | Variables
e_sqrt.c File Reference
#include <features.h>
#include <sysdeps/ieee754/dbl-64/e_sqrt.c>

Go to the source code of this file.

Classes

struct  sqrt_data_struct

Defines

#define __ieee754_sqrt   __full_ieee754_sqrt

Functions

static const struct
sqrt_data_struct 
__attribute__ ((used))
 asm ("\ /* Define offsets into the structure defined in C above. */ \n\ $DN = 0*8 \n\ $UP = 1*8 \n\ $HALF = 2*8 \n\ $ALMOST_THREE_HALF = 3*8 \n\ $NAN = 7*8 \n\ $T2 = 8*8 \n\ \n\ /* Stack variables. */ \n\ $K = 0 \n\ $Y = 8 \n\ \n\ .text \n\ .align 5 \n\ .globl __ieee754_sqrt \n\ .ent __ieee754_sqrt \n\ __ieee754_sqrt: \n\ ldgp $29, 0($27) \n\ subq $sp, 16, $sp \n\ .frame $sp, 16, $26, 0\n"" .prologue 1 \n\ \n\ .align 4 \n\ stt $f16, $K($sp) # e0 : \n\ mult $f31, $f31, $f31 # .. fm : \n\ lda $4, sqrt_data # e0 : \n\ fblt $f16, $fixup # .. fa : \n\ \n\ ldah $2, 0x5fe8 # e0 : \n\ ldq $3, $K($sp) # .. e1 : \n\ ldt $f12, $HALF($4) # e0 : \n\ ldt $f18, $ALMOST_THREE_HALF($4) # .. e1 : \n\ \n\ sll $3, 52, $5 # e0 : \n\ lda $6, 0x7fd # .. e1 : \n\ fnop # .. fa : \n\ fnop # .. fm : \n\ \n\ subq $5, 1, $5 # e1 : \n\ srl $3, 33, $1 # .. e0 : \n\ cmpule $5, $6, $5 # e0 : \n\ beq $5, $fixup # .. e1 : \n\ \n\ mult $f16, $f12, $f11 # fm : $f11 = x * 0.5 \n\ subl $2, $1, $2 # .. e0 : \n\ addt $f12, $f12, $f17 # .. fa : $f17 = 1.0 \n\ srl $2, 12, $1 # e0 : \n\ \n\ and $1, 0xfc, $1 # e0 : \n\ addq $1, $4, $1 # e1 : \n\ ldl $1, $T2($1) # e0 : \n\ addt $f12, $f17, $f15 # .. fa : $f15 = 1.5 \n\ \n\ subl $2, $1, $2 # e0 : \n\ ldt $f14, $DN($4) # .. e1 : \n\ sll $2, 32, $2 # e0 : \n\ stq $2, $Y($sp) # e0 : \n\ \n\ ldt $f13, $Y($sp) # e0 : \n\ mult/su $f11, $f13, $f10 # fm 2: $f10 = (x * 0.5) * y \n\ mult $f10, $f13, $f10 # fm 4: $f10 = ((x*0.5)*y)*y \n\ subt $f15, $f10, $f1 # fa 4: $f1 = (1.5-0.5*x*y*y) \n\ \n\ mult $f13, $f1, $f13 # fm 4: yp = y*(1.5-0.5*x*y^2)\n\ mult/su $f11, $f13, $f1 # fm 4: $f11 = x * 0.5 * yp \n\ mult $f1, $f13, $f11 # fm 4: $f11 = (x*0.5*yp)*yp \n\ subt $f18, $f11, $f1 # fa 4: $f1=(1.5-2^-30)-x/2*yp^2\n\ \n\ mult $f13, $f1, $f13 # fm 4: ypp = $f13 = yp*$f1 \n\ subt $f15, $f12, $f1 # .. fa : $f1 = (1.5 - 0.5) \n\ ldt $f15, $UP($4) # .. e0 : \n\ mult/su $f16, $f13, $f10 # fm 4: z = $f10 = x * ypp \n\ \n\ mult $f10, $f13, $f11 # fm 4: $f11 = z*ypp \n\ mult $f10, $f12, $f12 # fm : $f12 = z*0.5 \n\ subt $f1, $f11, $f1 # fa 4: $f1 = 1 - z*ypp \n\ mult $f12, $f1, $f12 # fm 4: $f12 = z/2*(1 - z*ypp)\n\ \n\ addt $f10, $f12, $f0 # fa 4: zp=res= z+z/2*(1-z*ypp)\n\ mult/c $f0, $f14, $f12 # fm 4: zmi = zp * DN \n\ mult/c $f0, $f15, $f11 # fm : zpl = zp * UP \n\ mult/c $f0, $f12, $f1 # fm : $f1 = zp * zmi \n\ \n\ mult/c $f0, $f11, $f15 # fm : $f15 = zp * zpl \n\ subt/su $f1, $f16, $f13 # .. fa : y1 = zp*zmi - x \n\ subt/su $f15, $f16, $f14 # fa 4: y2 = zp*zpl - x \n\ fcmovge $f13, $f12, $f0 # fa 3: res = (y1>=0)?zmi:res \n\ \n\ fcmovlt $f14, $f11, $f0 # fa 4: res = (y2<0)?zpl:res \n\ addq $sp, 16, $sp # .. e0 : \n\ ret # .. e1 : \n\ \n\ .align 4 \n\ $fixup: \n\ addq $sp, 16, $sp \n\ br __full_ieee754_sqrt !samegp \n\ \n\ .end __ieee754_sqrt")

Variables

unsigned long dn
unsigned long up
unsigned long half
unsigned long almost_three_half
unsigned long one_and_a_half
unsigned long two_to_minus_30
unsigned long one
unsigned long nan
const int T2 [64]

Class Documentation

struct sqrt_data_struct

Definition at line 32 of file e_sqrt.c.

Class Members
unsigned long almost_three_half
unsigned long dn
unsigned long half
unsigned long nan
unsigned long one
unsigned long one_and_a_half
const int T2
unsigned long two_to_minus_30
unsigned long up

Define Documentation

#define __ieee754_sqrt   __full_ieee754_sqrt

Function Documentation

static const struct sqrt_data_struct __attribute__ ( (used)  ) [static]
asm ( "\ /* Define offsets into the structure defined in C above. */ \n\ $  DN = 0*8 \n\ $UP = 1*8 \n\ $HALF = 2*8 \n\ $ALMOST_THREE_HALF = 3*8 \n\ $NAN = 7*8 \n\ $T2 = 8*8 \n\ \n\  /* Stack variables.  */ \n\ $K = 0 \n\ $Y = 8 \n\ \n\ .text \n\ .align 5 \n\ .globl __ieee754_sqrt \n\ .ent __ieee754_sqrt \n\__ieee754_sqrt: \n\ ldgp $29,
0($27)\n\subq $  sp,
16  ,
$sp\n\.frame $  sp,
16  ,
26,
0\n"".prologue 1\n\\n\.align 4\n\stt $  f16,
$K($sp)#e0:\n\mult $  f31,
f31,
$f31#..fm:\n\lda $  4,
sqrt_data#e0:\n\fblt $  f16,
$fixup#..fa:\n\\n\ldah $  2,
0x5fe8#e0:\n\ldq $  3,
$K($sp)#..e1:\n\ldt $  f12,
$HALF($4)#e0:\n\ldt $  f18,
$ALMOST_THREE_HALF($4)#..e1:\n\\n\sll $  3,
52  ,
$5#e0:\n\lda $  6,
0x7fd#..e1:\n\fnop#..fa:\n\fnop#..fm:\n\\n\subq $  5,
,
$5#e1:\n\srl $  3,
33  ,
$1#..e0:\n\cmpule $  5,
6,
$5#e0:\n\beq $  5,
$fixup#..e1:\n\\n\mult $  f16,
f12,
$f11#fm:$  f11 = x * 0.5 \n\ subl $2,
1,
$2#..e0:\n\addt $  f12,
f12,
$f17#..fa:$  f17 = 1.0 \n\ srl $2,
12  ,
$1#e0:\n\\n\and $  1,
0xfc  ,
$1#e0:\n\addq $  1,
4,
$1#e1:\n\ldl $  1,
$T2($1)#e0:\n\addt $  f12,
f17,
$f15#..fa:$  f15 = 1.5 \n\ \n\ subl $2,
1,
$2#e0:\n\ldt $  f14,
$DN($4)#..e1:\n\sll $  2,
32  ,
$2#e0:\n\stq $  2,
$Y($sp)#e0:\n\\n\ldt $  f13,
$Y($sp)#e0:\n\mult/su $  f11,
f13,
$f10#fm 2:$  f10 = (x * 0.5) * y \n\ mult $f10,
f13,
$f10#fm 4:$  f10 = ((x*0.5)*y)*y \n\ subt $f15,
f10,
$f1#fa 4:$  f1 = (1.5-0.5*x*y*y) \n\ \n\ mult $f13,
f1,
$f13#fm 4:yp  = y*(1.5-0.5*x*y^2)\n\  mult/su $f11,
f13,
$f1#fm 4:$  f11 = x * 0.5 * yp \n\ mult $f1,
f13,
$f11#fm 4:$  f11 = (x*0.5*yp)*yp \n\ subt $f18,
f11,
$f1#fa 4:$  f1 = (1.5-2^-30)-x/2*yp^2\n\ \n\ mult $f13,
f1,
$f13#fm 4:ypp  = $f13 = yp*$f1 \n\ subt $f15,
f12,
$f1#..fa:$  f1 = (1.5 - 0.5) \n\ ldt $f15,
$UP($4)#..e0:\n\mult/su $  f16,
f13,
$f10#fm 4:z  = $f10 = x * ypp \n\ \n\ mult $f10,
f13,
$f11#fm 4:$  f11 = z*ypp \n\ mult $f10,
f12,
$f12#fm:$  f12 = z*0.5 \n\ subt $f1,
f11,
$f1#fa 4:$  f1 = 1 - z*ypp \n\ mult $f12,
f1,
$f12#fm 4:$  f12 = z/2*(1 - z*ypp)\n\ \n\ addt $f10,
f12,
$f0#fa 4:zp  = resz+z/2*(1-z*ypp)\n\ mult/c $f0,
f14,
$f12#fm 4:zmi  = zp * DN \n\ mult/c $f0,
f15,
$f11#fm:zpl  = zp * UP \n\ mult/c $f0,
f12,
$f1#fm:$  f1 = zp * zmi \n\ \n\ mult/c $f0,
f11,
$f15#fm:$  f15 = zp * zpl \n\ subt/su $f1,
f16,
$f13#..fa:y1  = zp*zmi - x \n\ subt/su $f15,
f16,
$f14#fa 4:y2  = zp*zpl - x \n\ fcmovge $f13,
f12,
$f0#fa 3:res  = (y1>=0)?zmi:res \n\ \n\ fcmovlt $f14,
f11,
$f0#fa 4:res  = (y2<0)?zpl:res \n\ addq $sp,
16  ,
$sp#..e0:\n\ret#..e1:\n\\n\.align 4\n\$fixup:\n\addq $  sp,
16  ,
$sp\n\br __full_ieee754_sqrt!samegp\n\\n\.end __ieee754_sqrt  
)

Variable Documentation

unsigned long almost_three_half

Definition at line 55 of file e_sqrt.c.

unsigned long dn

Definition at line 55 of file e_sqrt.c.

unsigned long half

Definition at line 55 of file e_sqrt.c.

unsigned long nan

Definition at line 56 of file e_sqrt.c.

unsigned long one

Definition at line 56 of file e_sqrt.c.

unsigned long one_and_a_half

Definition at line 56 of file e_sqrt.c.

const int T2[64]

Definition at line 57 of file e_sqrt.c.

unsigned long two_to_minus_30

Definition at line 56 of file e_sqrt.c.

unsigned long up

Definition at line 55 of file e_sqrt.c.