| 241 |
|
register int x; |
| 242 |
|
/* compute foveal y position */ |
| 243 |
|
iy = dy = (y+.5)/numscans(&inpres)*fvyr - .5; |
| 244 |
< |
if (iy >= fvyr-1) iy--; |
| 244 |
> |
while (iy >= fvyr-1) iy--; |
| 245 |
|
dy -= (double)iy; |
| 246 |
|
for (x = 0; x < scanlen(&inpres); x++) { |
| 247 |
|
/* compute foveal x position */ |
| 248 |
|
ix = dx = (x+.5)/scanlen(&inpres)*fvxr - .5; |
| 249 |
< |
if (ix >= fvxr-1) ix--; |
| 249 |
> |
while (ix >= fvxr-1) ix--; |
| 250 |
|
dx -= (double)ix; |
| 251 |
|
/* interpolate sample rate */ |
| 252 |
|
sr = (1.-dy)*((1.-dx)*tsampr(ix,iy) + dx*tsampr(ix+1,iy)) + |
| 302 |
|
dy -= (double)iy; |
| 303 |
|
/* get scanlines */ |
| 304 |
|
sl0 = getascan(sb, iy); |
| 305 |
+ |
#ifdef DEBUG |
| 306 |
|
if (sl0 == NULL) { |
| 307 |
|
fprintf(stderr, "%s: internal - cannot backspace in ascanval\n", |
| 308 |
|
progname); |
| 309 |
< |
exit(1); |
| 309 |
> |
abort(); |
| 310 |
|
} |
| 311 |
+ |
#endif |
| 312 |
|
sl1 = getascan(sb, iy+1); |
| 313 |
|
/* 2D linear interpolation */ |
| 314 |
|
copycolor(col, sl0[ix]); |
| 324 |
|
scalecolor(col, 1.-dy); |
| 325 |
|
scalecolor(c1y, dy); |
| 326 |
|
addcolor(col, c1y); |
| 327 |
+ |
for (ix = 0; ix < 3; ix++) /* make sure no negative */ |
| 328 |
+ |
if (colval(col,ix) < 0.) |
| 329 |
+ |
colval(col,ix) = 0.; |
| 330 |
|
} |
| 331 |
|
|
| 332 |
|
|
| 333 |
|
SCANBAR * |
| 334 |
< |
sballoc(sr, ns, sl) /* allocate scanbar */ |
| 335 |
< |
int sr; /* sampling rate */ |
| 334 |
> |
sballoc(se, ns, sl) /* allocate scanbar */ |
| 335 |
> |
int se; /* sampling rate exponent */ |
| 336 |
|
int ns; /* number of scanlines */ |
| 337 |
|
int sl; /* original scanline length */ |
| 338 |
|
{ |
| 339 |
|
SCANBAR *sbarr; |
| 340 |
|
register SCANBAR *sb; |
| 341 |
|
|
| 342 |
< |
sbarr = sb = (SCANBAR *)malloc((sr+1)*sizeof(SCANBAR)); |
| 342 |
> |
sbarr = sb = (SCANBAR *)malloc((se+1)*sizeof(SCANBAR)); |
| 343 |
|
if (sb == NULL) |
| 344 |
|
syserror("malloc"); |
| 345 |
|
do { |
| 346 |
< |
sb->sdata = (COLOR *)malloc((sl>>sr)*ns*sizeof(COLOR)); |
| 346 |
> |
sb->sampe = se; |
| 347 |
> |
sb->len = sl>>se; |
| 348 |
> |
sb->nscans = ns; |
| 349 |
> |
sb->sdata = (COLOR *)malloc(sb->len*ns*sizeof(COLOR)); |
| 350 |
|
if (sb->sdata == NULL) |
| 351 |
|
syserror("malloc"); |
| 344 |
– |
sb->sampe = sr; |
| 345 |
– |
sb->nscans = ns; |
| 346 |
– |
sb->len = sl>>sr; |
| 352 |
|
sb->nread = 0; |
| 353 |
|
ns <<= 1; |
| 354 |
|
sb++; |
| 355 |
< |
} while (--sr >= 0); |
| 355 |
> |
} while (--se >= 0); |
| 356 |
|
return(sbarr); |
| 357 |
|
} |
| 358 |
|
|
| 381 |
|
} |
| 382 |
|
fcross(cp, diffx, diffy); |
| 383 |
|
omega = 0.5 * sqrt(DOT(cp,cp)); |
| 384 |
< |
if (omega <= FTINY) |
| 384 |
> |
if (omega <= FTINY*FTINY) |
| 385 |
|
tsampr(x,y) = 1.; |
| 386 |
|
else if ((tsampr(x,y) = PI/180. / sqrt(omega) / |
| 387 |
|
hacuity(plum(fovscan(y)[x]))) > maxsr) |
| 393 |
|
tsampr(x,fvyr-1) = tsampr(x,fvyr-2); |
| 394 |
|
} |
| 395 |
|
for (y = 0; y < fvyr; y++) { |
| 396 |
< |
tsampr(y,0) = tsampr(y,1); |
| 397 |
< |
tsampr(y,fvxr-1) = tsampr(y,fvxr-2); |
| 396 |
> |
tsampr(0,y) = tsampr(1,y); |
| 397 |
> |
tsampr(fvxr-1,y) = tsampr(fvxr-2,y); |
| 398 |
|
} |
| 399 |
|
/* initialize with next power of two */ |
| 400 |
|
rootbar = sballoc((int)(log(maxsr)/log(2.))+1, 2, scanlen(&inpres)); |
