| 11 |
|
* |
| 12 |
|
*/ |
| 13 |
|
|
| 14 |
+ |
#define _USE_MATH_DEFINES |
| 15 |
|
#include "rtio.h" |
| 16 |
|
#include <stdlib.h> |
| 17 |
|
#include <math.h> |
| 87 |
|
return (a <= 1e-6) & (a >= -1e-6); |
| 88 |
|
} |
| 89 |
|
|
| 89 |
– |
/* Returns the name of the given tag */ |
| 90 |
– |
#ifdef ezxml_name |
| 91 |
– |
#undef ezxml_name |
| 92 |
– |
static char * |
| 93 |
– |
ezxml_name(ezxml_t xml) |
| 94 |
– |
{ |
| 95 |
– |
if (xml == NULL) |
| 96 |
– |
return NULL; |
| 97 |
– |
return xml->name; |
| 98 |
– |
} |
| 99 |
– |
#endif |
| 100 |
– |
|
| 90 |
|
/* Returns the given tag's character content or empty string if none */ |
| 91 |
|
#ifdef ezxml_txt |
| 92 |
|
#undef ezxml_txt |
| 178 |
|
{ |
| 179 |
|
ANGLE_BASIS *ab = (ANGLE_BASIS *)p; |
| 180 |
|
int li, ndx; |
| 181 |
< |
double pol, azi, d; |
| 181 |
> |
double pol, azi; |
| 182 |
|
|
| 183 |
|
if (v == NULL) |
| 184 |
|
return -1; |
| 218 |
|
if (li == last_li) /* cached latitude? */ |
| 219 |
|
return last_ohm; |
| 220 |
|
last_li = li; |
| 232 |
– |
theta1 = M_PI/180. * ab->lat[li+1].tmin; |
| 233 |
– |
if (ab->lat[li].nphis == 1) /* special case */ |
| 234 |
– |
return last_ohm = M_PI*(1. - sq(cos(theta1))); |
| 221 |
|
theta = M_PI/180. * ab->lat[li].tmin; |
| 222 |
+ |
theta1 = M_PI/180. * ab->lat[li+1].tmin; |
| 223 |
|
return last_ohm = M_PI*(sq(cos(theta)) - sq(cos(theta1))) / |
| 224 |
|
(double)ab->lat[li].nphis; |
| 225 |
|
} |
| 333 |
|
ezxml_child(ezxml_child(wbb, |
| 334 |
|
"ThetaBounds"), "UpperTheta"))); |
| 335 |
|
if (!i) |
| 336 |
< |
abase_list[nabases].lat[i].tmin = |
| 350 |
< |
-abase_list[nabases].lat[i+1].tmin; |
| 336 |
> |
abase_list[nabases].lat[0].tmin = 0; |
| 337 |
|
else if (!fequal(atof(ezxml_txt(ezxml_child(ezxml_child(wbb, |
| 338 |
|
"ThetaBounds"), "LowerTheta"))), |
| 339 |
|
abase_list[nabases].lat[i].tmin)) { |
| 514 |
|
if (*sdnext == ',') sdnext++; |
| 515 |
|
if (rowinc) { |
| 516 |
|
int r = i/dp->nout; |
| 517 |
< |
int c = i - c*dp->nout; |
| 517 |
> |
int c = i - r*dp->nout; |
| 518 |
|
mBSDF_value(dp,r,c) = atof(sdata); |
| 519 |
|
} else |
| 520 |
|
dp->bsdf[i] = atof(sdata); |
| 595 |
|
sd->name); |
| 596 |
|
return SDEsupport; |
| 597 |
|
} |
| 598 |
< |
/* get angle basis */ |
| 599 |
< |
rval = load_angle_basis(ezxml_child(ezxml_child(wtl, |
| 600 |
< |
"DataDefinition"), "AngleBasis")); |
| 601 |
< |
if (rval < 0) |
| 602 |
< |
return convert_errcode(rval); |
| 598 |
> |
/* get angle bases */ |
| 599 |
> |
for (wld = ezxml_child(ezxml_child(wtl, "DataDefinition"), "AngleBasis"); |
| 600 |
> |
wld != NULL; wld = wld->next) { |
| 601 |
> |
rval = load_angle_basis(wld); |
| 602 |
> |
if (rval < 0) |
| 603 |
> |
return convert_errcode(rval); |
| 604 |
> |
} |
| 605 |
|
/* load BSDF components */ |
| 606 |
|
for (wld = ezxml_child(wtl, "WavelengthData"); |
| 607 |
|
wld != NULL; wld = wld->next) { |
| 756 |
|
} |
| 757 |
|
cdlast = NULL; /* check for it in cache list */ |
| 758 |
|
for (cd = (SDMatCDst *)sdc->cdList; cd != NULL; |
| 759 |
< |
cdlast = cd, cd = (SDMatCDst *)cd->next) |
| 759 |
> |
cdlast = cd, cd = cd->next) |
| 760 |
|
if (cd->indx == myCD.indx && (cd->calen == myCD.calen) & |
| 761 |
|
(cd->ob_priv == myCD.ob_priv) & |
| 762 |
|
(cd->ob_vec == myCD.ob_vec)) |
| 775 |
|
} |
| 776 |
|
if (cdlast != NULL) { /* move entry to head of cache list */ |
| 777 |
|
cdlast->next = cd->next; |
| 778 |
< |
cd->next = sdc->cdList; |
| 778 |
> |
cd->next = (SDMatCDst *)sdc->cdList; |
| 779 |
|
sdc->cdList = (SDCDst *)cd; |
| 780 |
|
} |
| 781 |
|
return (SDCDst *)cd; /* ready to go */ |
| 795 |
|
/* binary search to find index */ |
| 796 |
|
ilower = 0; iupper = mcd->calen; |
| 797 |
|
while ((i = (iupper + ilower) >> 1) != ilower) |
| 798 |
< |
if ((long)target >= (long)mcd->carr[i]) |
| 798 |
> |
if (target >= mcd->carr[i]) |
| 799 |
|
ilower = i; |
| 800 |
|
else |
| 801 |
|
iupper = i; |
