ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/bsdf_m.c
(Generate patch)

Comparing ray/src/common/bsdf_m.c (file contents):
Revision 3.34 by greg, Sun Aug 23 00:56:00 2015 UTC vs.
Revision 3.44 by greg, Fri Dec 17 20:51:55 2021 UTC

# Line 171 | Line 171 | fo_getndx(const FVECT v, void *p)
171  
172          if (v == NULL)
173                  return -1;
174 <        if ((v[2] < 0) | (v[2] > 1.))
174 >        if ((v[2] < 0) | (v[2] > 1.00001))
175                  return -1;
176          pol = 180.0/M_PI*Acos(v[2]);
177          azi = 180.0/M_PI*atan2(v[1], v[0]);
# Line 191 | Line 191 | fo_getndx(const FVECT v, void *p)
191   double
192   io_getohm(int ndx, void *p)
193   {
194 +        static void     *last_p = NULL;
195          static int      last_li = -1;
196          static double   last_ohm;
197          ANGLE_BASIS     *ab = (ANGLE_BASIS *)p;
# Line 201 | Line 202 | io_getohm(int ndx, void *p)
202                  return -1.;
203          for (li = 0; ndx >= ab->lat[li].nphis; li++)
204                  ndx -= ab->lat[li].nphis;
205 <        if (li == last_li)                      /* cached latitude? */
205 >        if ((p == last_p) & (li == last_li))            /* cached latitude? */
206                  return last_ohm;
207 +        last_p = p;
208          last_li = li;
209          theta = M_PI/180. * ab->lat[li].tmin;
210          theta1 = M_PI/180. * ab->lat[li+1].tmin;
# Line 294 | Line 296 | mBSDF_color(float coef[], const SDMat *dp, int i, int
296   {
297          C_COLOR cxy;
298  
299 <        coef[0] = mBSDF_value(dp, i, o);
299 >        coef[0] = mBSDF_value(dp, o, i);
300          if (dp->chroma == NULL)
301                  return 1;       /* grayscale */
302  
303 <        c_decodeChroma(&cxy, dp->chroma[o*dp->ninc + i]);
303 >        c_decodeChroma(&cxy, mBSDF_chroma(dp,o,i));
304          c_toSharpRGB(&cxy, coef[0], coef);
305          coef[0] *= mtx_RGB_coef[0];
306          coef[1] *= mtx_RGB_coef[1];
# Line 381 | Line 383 | get_extrema(SDSpectralDF *df)
383          for (i = dp->ninc; i--; ) {
384                  double  hemi = .0;
385                  for (o = dp->nout; o--; )
386 <                        hemi += ohma[o] * mBSDF_value(dp, i, o);
386 >                        hemi += ohma[o] * mBSDF_value(dp, o, i);
387                  if (hemi > df->maxHemi)
388                          df->maxHemi = hemi;
389          }
# Line 517 | Line 519 | load_bsdf_data(SDData *sd, ezxml_t wdb, int ct, int ro
519                  if (rowinc) {
520                          int     r = i/dp->nout;
521                          int     c = i - r*dp->nout;
522 <                        mBSDF_value(dp,r,c) = val;
522 >                        mBSDF_value(dp,c,r) = val;
523                  } else
524                          dp->bsdf[i] = val;
525                  sdata = sdnext;
# Line 622 | Line 624 | subtract_min(C_COLOR *cs, SDMat *sm)
624                                  coef[c] = (coef[c] - min_coef[c]) /
625                                                  mtx_RGB_coef[c];
626                          if (c_fromSharpRGB(coef, &cxy) > 1e-5)
627 <                                sm->chroma[o*sm->ninc + i] = c_encodeChroma(&cxy);
628 <                        mBSDF_value(sm,i,o) -= ymin;
627 >                                mBSDF_chroma(sm,o,i) = c_encodeChroma(&cxy);
628 >                        mBSDF_value(sm,o,i) -= ymin;
629                  }
630                                          /* return colored minimum */
631          for (i = 3; i--; )
# Line 647 | Line 649 | extract_diffuse(SDValue *dv, SDSpectralDF *df)
649                                          /* subtract minimum value */
650          dv->cieY = subtract_min(&dv->spec, (SDMat *)df->comp[0].dist);
651          df->maxHemi -= dv->cieY;        /* adjust maximum hemispherical */
652 <                                        /* make sure everything is set */
653 <        c_ccvt(&dv->spec, C_CSXY+C_CSSPEC);
652 >                                
653 >        c_ccvt(&dv->spec, C_CSXY);      /* make sure (x,y) is set */
654          return df;
655   }
656  
# Line 707 | Line 709 | SDloadMtx(SDData *sd, ezxml_t wtl)
709                                          /* separate diffuse components */
710          sd->rf = extract_diffuse(&sd->rLambFront, sd->rf);
711          sd->rb = extract_diffuse(&sd->rLambBack, sd->rb);
712 <        if (sd->tf != NULL)
713 <                sd->tf = extract_diffuse(&sd->tLamb, sd->tf);
714 <        if (sd->tb != NULL)
715 <                sd->tb = extract_diffuse(&sd->tLamb, sd->tb);
712 >        sd->tf = extract_diffuse(&sd->tLambFront, sd->tf);
713 >        if (sd->tb != NULL) {
714 >                sd->tb = extract_diffuse(&sd->tLambBack, sd->tb);
715 >                if (sd->tf == NULL)
716 >                        sd->tLambFront = sd->tLambBack;
717 >        } else if (sd->tf != NULL)
718 >                sd->tLambBack = sd->tLambFront;
719                                          /* return success */
720          return SDEnone;
721   }
722  
723   /* Get Matrix BSDF value */
724   static int
725 < SDgetMtxBSDF(float coef[SDmaxCh], const FVECT outVec,
726 <                                const FVECT inVec, SDComponent *sdc)
725 > SDgetMtxBSDF(float coef[SDmaxCh], const FVECT inVec,
726 >                                const FVECT outVec, SDComponent *sdc)
727   {
728          const SDMat     *dp;
729          int             i_ndx, o_ndx;
# Line 803 | Line 808 | make_cdist(SDMatCDst *cd, const FVECT inVec, SDMat *dp
808          cmtab[0] = .0;
809          for (o = 0; o < cd->calen; o++) {
810                  if (rev)
811 <                        cmtab[o+1] = mBSDF_value(dp, o, cd->indx) *
811 >                        cmtab[o+1] = mBSDF_value(dp, cd->indx, o) *
812                                          (*dp->ib_ohm)(o, dp->ib_priv);
813                  else
814 <                        cmtab[o+1] = mBSDF_value(dp, cd->indx, o) *
814 >                        cmtab[o+1] = mBSDF_value(dp, o, cd->indx) *
815                                          (*dp->ob_ohm)(o, dp->ob_priv);
816                  cmtab[o+1] += cmtab[o];
817          }
# Line 906 | Line 911 | SDsampMtxCDist(FVECT ioVec, double randX, const SDCDst
911   }
912  
913   /* Fixed resolution BSDF methods */
914 < SDFunc                  SDhandleMtx = {
914 > const SDFunc            SDhandleMtx = {
915                                  &SDgetMtxBSDF,
916                                  &SDqueryMtxProjSA,
917                                  &SDgetMtxCDist,

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines