- Timestamp:
- Aug 25, 2018, 6:30:55 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/libs/libmath/k_rem_pio2.c
r469 r511 168 168 /* set up f[0] to f[jx+jk] where f[jx+jk] = ipio2[jv+jk] */ 169 169 j = jv-jx; m = jx+jk; 170 for (i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j];170 for (i=0;i<=m;i++,j++) { f[i] = (j<0)? zero : (double) ipio2[j]; } 171 171 172 172 /* compute q[0],q[1],...q[jk] */ 173 for (i=0;i<=jk;i++) { 174 for(j=0,fw=0.0;j<=jx;j++) fw += x[j]*f[jx+i-j]; q[i] = fw; 173 for ( i = 0; i <= jk; i++ ) { 174 for ( j = 0, fw = 0.0; j <= jx; j++) { 175 fw += x[j] * f[jx + i - j]; 176 } 177 q[i] = fw; 175 178 } 176 179
Note: See TracChangeset
for help on using the changeset viewer.