| 32 |
|
static LUTAB mat_tab = LU_SINIT(free,free); /* material lookup table */ |
| 33 |
|
static LUTAB vtx_tab = LU_SINIT(free,free); /* vertex lookup table */ |
| 34 |
|
|
| 35 |
< |
static int setspectrum(C_COLOR *, double wlmin, double wlmax, |
| 36 |
< |
int ac, char **av); |
| 35 |
> |
static int setspectrum(C_COLOR *, int ac, char **av); |
| 36 |
|
|
| 37 |
|
|
| 38 |
|
int |
| 39 |
|
c_hcolor(int ac, char **av) /* handle color entity */ |
| 40 |
|
{ |
| 41 |
|
double w, wsum; |
| 42 |
< |
register int i; |
| 43 |
< |
register LUENT *lp; |
| 42 |
> |
int i; |
| 43 |
> |
LUENT *lp; |
| 44 |
|
|
| 45 |
|
switch (mg_entity(av[0])) { |
| 46 |
|
case MG_E_COLOR: /* get/set color context */ |
| 107 |
|
c_ccolor->clock++; |
| 108 |
|
return(MG_OK); |
| 109 |
|
case MG_E_CSPEC: /* assign spectral values */ |
| 110 |
< |
if (ac < 5) |
| 112 |
< |
return(MG_EARGC); |
| 113 |
< |
if (!isflt(av[1]) | !isflt(av[2])) |
| 114 |
< |
return(MG_ETYPE); |
| 115 |
< |
return(setspectrum(c_ccolor, atof(av[1]), atof(av[2]), |
| 116 |
< |
ac-3, av+3)); |
| 110 |
> |
return(setspectrum(c_ccolor, ac-1, av+1)); |
| 111 |
|
case MG_E_CCT: /* assign black body spectrum */ |
| 112 |
|
if (ac != 2) |
| 113 |
|
return(MG_EARGC); |
| 153 |
|
c_hmaterial(int ac, char **av) /* handle material entity */ |
| 154 |
|
{ |
| 155 |
|
int i; |
| 156 |
< |
register LUENT *lp; |
| 156 |
> |
LUENT *lp; |
| 157 |
|
|
| 158 |
|
switch (mg_entity(av[0])) { |
| 159 |
|
case MG_E_MATERIAL: /* get/set material context */ |
| 299 |
|
c_hvertex(int ac, char **av) /* handle a vertex entity */ |
| 300 |
|
{ |
| 301 |
|
int i; |
| 302 |
< |
register LUENT *lp; |
| 302 |
> |
LUENT *lp; |
| 303 |
|
|
| 304 |
|
switch (mg_entity(av[0])) { |
| 305 |
|
case MG_E_VERTEX: /* get/set vertex context */ |
| 399 |
|
C_MATERIAL * |
| 400 |
|
c_getmaterial(char *name) /* get a named material */ |
| 401 |
|
{ |
| 402 |
< |
register LUENT *lp; |
| 402 |
> |
LUENT *lp; |
| 403 |
|
|
| 404 |
|
if ((lp = lu_find(&mat_tab, name)) == NULL) |
| 405 |
|
return(NULL); |
| 410 |
|
C_VERTEX * |
| 411 |
|
c_getvert(char *name) /* get a named vertex */ |
| 412 |
|
{ |
| 413 |
< |
register LUENT *lp; |
| 413 |
> |
LUENT *lp; |
| 414 |
|
|
| 415 |
|
if ((lp = lu_find(&vtx_tab, name)) == NULL) |
| 416 |
|
return(NULL); |
| 421 |
|
C_COLOR * |
| 422 |
|
c_getcolor(char *name) /* get a named color */ |
| 423 |
|
{ |
| 424 |
< |
register LUENT *lp; |
| 424 |
> |
LUENT *lp; |
| 425 |
|
|
| 426 |
|
if ((lp = lu_find(&clr_tab, name)) == NULL) |
| 427 |
|
return(NULL); |
| 430 |
|
|
| 431 |
|
|
| 432 |
|
static int |
| 433 |
< |
setspectrum( /* convert a spectrum */ |
| 434 |
< |
register C_COLOR *clr, |
| 441 |
< |
double wlmin, |
| 442 |
< |
double wlmax, |
| 433 |
> |
setspectrum( /* convert a spectrum */ |
| 434 |
> |
C_COLOR *clr, |
| 435 |
|
int ac, |
| 436 |
|
char **av |
| 437 |
|
) |
| 438 |
|
{ |
| 439 |
< |
double scale; |
| 440 |
< |
float va[C_CNSS]; |
| 441 |
< |
register int i, pos; |
| 442 |
< |
int n, imax; |
| 443 |
< |
int wl; |
| 444 |
< |
double wl0, wlstep; |
| 445 |
< |
double boxpos, boxstep; |
| 446 |
< |
/* check bounds */ |
| 447 |
< |
if ((wlmax <= C_CMINWL) | (wlmax <= wlmin) | (wlmin >= C_CMAXWL)) |
| 448 |
< |
return(MG_EILL); |
| 449 |
< |
wlstep = (wlmax - wlmin)/(ac-1); |
| 450 |
< |
while (wlmin < C_CMINWL) { |
| 451 |
< |
wlmin += wlstep; |
| 452 |
< |
ac--; av++; |
| 439 |
> |
double wlmin, wlmax; |
| 440 |
> |
int i; |
| 441 |
> |
float *va; |
| 442 |
> |
/* get bounds */ |
| 443 |
> |
if (ac < 4) |
| 444 |
> |
return(MG_EARGC); |
| 445 |
> |
if (!isflt(av[0]) || !isflt(av[1])) |
| 446 |
> |
return(MG_ETYPE); |
| 447 |
> |
wlmin = atof(av[0]); |
| 448 |
> |
wlmax = atof(av[1]); |
| 449 |
> |
ac -= 2; av += 2; |
| 450 |
> |
va = (float *)malloc(sizeof(float)*ac); |
| 451 |
> |
if (va == NULL) |
| 452 |
> |
return(MG_EMEM); |
| 453 |
> |
for (i = ac; i--; ) { |
| 454 |
> |
if (!isflt(av[i])) |
| 455 |
> |
return(MG_ETYPE); |
| 456 |
> |
va[i] = atof(av[i]); |
| 457 |
|
} |
| 458 |
< |
while (wlmax > C_CMAXWL) { |
| 459 |
< |
wlmax -= wlstep; |
| 464 |
< |
ac--; |
| 465 |
< |
} |
| 466 |
< |
imax = ac; /* box filter if necessary */ |
| 467 |
< |
boxpos = 0; |
| 468 |
< |
boxstep = 1; |
| 469 |
< |
if (wlstep < C_CWLI) { |
| 470 |
< |
imax = (wlmax - wlmin)/C_CWLI + (1-FTINY); |
| 471 |
< |
boxpos = (wlmin - C_CMINWL)/C_CWLI; |
| 472 |
< |
boxstep = wlstep/C_CWLI; |
| 473 |
< |
wlstep = C_CWLI; |
| 474 |
< |
} |
| 475 |
< |
scale = 0.; /* get values and maximum */ |
| 476 |
< |
pos = 0; |
| 477 |
< |
for (i = 0; i < imax; i++) { |
| 478 |
< |
va[i] = 0.; n = 0; |
| 479 |
< |
while (boxpos < i+.5 && pos < ac) { |
| 480 |
< |
if (!isflt(av[pos])) |
| 481 |
< |
return(MG_ETYPE); |
| 482 |
< |
va[i] += atof(av[pos++]); |
| 483 |
< |
n++; |
| 484 |
< |
boxpos += boxstep; |
| 485 |
< |
} |
| 486 |
< |
if (n > 1) |
| 487 |
< |
va[i] /= (double)n; |
| 488 |
< |
if (va[i] > scale) |
| 489 |
< |
scale = va[i]; |
| 490 |
< |
else if (va[i] < -scale) |
| 491 |
< |
scale = -va[i]; |
| 492 |
< |
} |
| 493 |
< |
if (scale <= FTINY) |
| 458 |
> |
if (c_sset(clr, wlmin, wlmax, va, ac) <= FTINY) { |
| 459 |
> |
free(va); |
| 460 |
|
return(MG_EILL); |
| 461 |
< |
scale = C_CMAXV / scale; |
| 462 |
< |
clr->ssum = 0; /* convert to our spacing */ |
| 497 |
< |
wl0 = wlmin; |
| 498 |
< |
pos = 0; |
| 499 |
< |
for (i = 0, wl = C_CMINWL; i < C_CNSS; i++, wl += C_CWLI) |
| 500 |
< |
if ((wl < wlmin) | (wl > wlmax)) |
| 501 |
< |
clr->ssamp[i] = 0; |
| 502 |
< |
else { |
| 503 |
< |
while (wl0 + wlstep < wl+FTINY) { |
| 504 |
< |
wl0 += wlstep; |
| 505 |
< |
pos++; |
| 506 |
< |
} |
| 507 |
< |
if ((wl+FTINY >= wl0) & (wl-FTINY <= wl0)) |
| 508 |
< |
clr->ssamp[i] = scale*va[pos] + .5; |
| 509 |
< |
else /* interpolate if necessary */ |
| 510 |
< |
clr->ssamp[i] = .5 + scale / wlstep * |
| 511 |
< |
( va[pos]*(wl0+wlstep - wl) + |
| 512 |
< |
va[pos+1]*(wl - wl0) ); |
| 513 |
< |
clr->ssum += clr->ssamp[i]; |
| 514 |
< |
} |
| 515 |
< |
clr->flags = C_CDSPEC|C_CSSPEC; |
| 461 |
> |
} |
| 462 |
> |
free(va); |
| 463 |
|
clr->clock++; |
| 464 |
|
return(MG_OK); |
| 465 |
|
} |
