Home | History | Annotate | Download | only in libm

Lines Matching defs:jz

62 	__int32_t jz,jx,jv,jp,jk,carry,n,iq[20],i,j,k,m,q0,ih;
83 jz = jk;
86 for(i=0,j=jz,z=q[jz];j>0;i++,j--) {
98 if(q0>0) { /* need iq[jz-1] to determine n */
99 i = (iq[jz-1]>>(8-q0)); n += i;
100 iq[jz-1] -= i<<(8-q0);
101 ih = iq[jz-1]>>(7-q0);
103 else if(q0==0) ih = iq[jz-1]>>8;
108 for(i=0;i<jz ;i++) { /* compute 1-q */
119 iq[jz-1] &= 0x7f; break;
121 iq[jz-1] &= 0x3f; break;
133 for (i=jz-1;i>=jk;i--) j |= iq[i];
137 for(i=jz+1;i<=jz+k;i++) { /* add q[jz+1] to q[jz+k] */
142 jz += k;
149 jz -= 1; q0 -= 8;
150 while(iq[jz]==0) { jz--; q0-=8;}
155 iq[jz] = (__int32_t)(z-two8*fw);
156 jz += 1; q0 += 8;
157 iq[jz] = (__int32_t) fw;
158 } else iq[jz] = (__int32_t) z ;
163 for(i=jz;i>=0;i--) {
167 /* compute PIo2[0,...,jp]*q[jz,...,0] */
168 for(i=jz;i>=0;i--) {
169 for(fw=0.0,k=0;k<=jp&&k<=jz-i;k++) fw += PIo2[k]*q[i+k];
170 fq[jz-i] = fw;
177 for (i=jz;i>=0;i--) fw += fq[i];
183 for (i=jz;i>=0;i--) fw += fq[i];
186 for (i=1;i<=jz;i++) fw += fq[i];
190 for (i=jz;i>0;i--) {
195 for (i=jz;i>1;i--) {
200 for (fw=0.0,i=jz;i>=2;i--) fw += fq[i];