Lines Matching refs:J
5074 #--LET J BE THE LEAST SIG. BIT OF D0, LET SGN := (-1)**J.
6737 # 2.3 Calculate J = N mod 64; so J = 0,1,2,..., #
6739 # 2.4 Calculate M = (N - J)/64; so N = 64M + J. #
6741 # 2^(J/64). #
6799 # Step 5. Compute 2^(J/64)*exp(R) = 2^(J/64)*(1+p) by #
6801 # where T and t are the stored values for 2^(J/64). #
6802 # Notes: 2^(J/64) is stored as T and t where T+t approximates #
6803 # 2^(J/64) to roughly 85 bits; T is in extended precision #
6812 # exp(X) = 2^M * 2^(J/64) * exp(R). #
6848 # 8.3 Calculate J = N mod 64, J = 0,1,...,63 #
6849 # 8.4 K := (N-J)/64, M1 := truncate(K/2), M = K-M1, #
6852 # 2^(J/64). #
7132 and.l &0x3F,%d1 # D0 is J = N mod 64
7134 add.l %d1,%a1 # address of 2^(J/64)
7143 #--a0 points to 2^(J/64), D0 is biased expo. of 2^(M)
7181 fmov.x (%a1)+,%fp1 # fp1 is lead. pt. of 2^(J/64)
7186 #--EXP(X) = 2^M * ( 2^(J/64) + 2^(J/64)*(EXP(R)-1) )
7188 fmul.x %fp1,%fp0 # 2^(J/64)*(Exp(R)-1)
7190 fadd.s (%a1),%fp0 # accurate 2^(J/64)
7192 fadd.x %fp1,%fp0 # 2^(J/64) + 2^(J/64)*...
7228 and.l &0x3F,%d1 # D0 is J = N mod 64
7230 add.l %d1,%a1 # address of 2^(J/64)
7296 and.l &0x3F,%d1 # D0 is J = N mod 64
7298 add.l %d1,%a1 # address of 2^(J/64)
7305 #--a0 points to 2^(J/64), D0 and a1 both contain M
7357 #--Compute 2^(J/64)*p
7359 fmul.x (%a1),%fp0 # 2^(J/64)*(Exp(R)-1)
9009 lea TEXPTBL(%pc),%a1 # LOAD ADDRESS OF TABLE OF 2^(J/64)
9013 and.l &0x3F,%d1 # D0 IS J
9014 asl.l &4,%d1 # DISPLACEMENT FOR 2^(J/64)
9015 add.l %d1,%a1 # ADDRESS FOR 2^(J/64)
9016 asr.l &6,%d2 # d2 IS L, N = 64L + J
9019 sub.l %d1,%d2 # d2 IS M', N = 64(M+M') + J
9022 #--SUMMARY: a1 IS ADDRESS FOR THE LEADING PORTION OF 2^(J/64),
9106 lea TEXPTBL(%pc),%a1 # LOAD ADDRESS OF TABLE OF 2^(J/64)
9110 and.l &0x3F,%d1 # D0 IS J
9111 asl.l &4,%d1 # DISPLACEMENT FOR 2^(J/64)
9112 add.l %d1,%a1 # ADDRESS FOR 2^(J/64)
9113 asr.l &6,%d2 # d2 IS L, N = 64L + J
9116 sub.l %d1,%d2 # d2 IS M', N = 64(M+M') + J
9119 #--SUMMARY: a1 IS ADDRESS FOR THE LEADING PORTION OF 2^(J/64),
9151 #-- 2**(M'+M) * 2**(J/64) * EXP(R)
9178 #--EXP(X) = 2^M*2^(J/64) + 2^M*2^(J/64)*(EXP(R)-1) - (1 OR 0)