src/libm/k_rem_pio2.c
changeset 11838 5ef6e4e70103
parent 11683 48bcba563d9c
child 12068 ce88faaf8bd2
equal deleted inserted replaced
11837:25872a19809a 11838:5ef6e4e70103
   165 	j = jv-jx; m = jx+jk;
   165 	j = jv-jx; m = jx+jk;
   166 	for(i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j];
   166 	for(i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j];
   167 
   167 
   168     /* compute q[0],q[1],...q[jk] */
   168     /* compute q[0],q[1],...q[jk] */
   169 	for (i=0;i<=jk;i++) {
   169 	for (i=0;i<=jk;i++) {
   170 	    for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j]; q[i] = fw;
   170 	    for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j];
       
   171 	    q[i] = fw;
   171 	}
   172 	}
   172 
   173 
   173 	jz = jk;
   174 	jz = jk;
   174 recompute:
   175 recompute:
   175     /* distill q[] into iq[] reversingly */
   176     /* distill q[] into iq[] reversingly */