ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/Development/ray/src/cv/bsdfmesh.c
(Generate patch)

Comparing ray/src/cv/bsdfmesh.c (file contents):
Revision 2.31 by greg, Thu Aug 21 13:44:05 2014 UTC vs.
Revision 2.34 by greg, Fri Jan 29 16:21:55 2016 UTC

# Line 51 | Line 51 | eval_DSFsurround(const RBFNODE *rbf, const FVECT outve
51                  if (i) spinvector(tvec, tvec, outvec, phinc);
52                  if (tvec[2] > 0 ^ output_orient > 0)
53                          continue;
54 <                sum += eval_rbfrep(rbf, tvec) * output_orient*tvec[2];
54 >                sum += eval_rbfrep(rbf, tvec) * COSF(tvec[2]);
55                  ++n;
56          }
57          if (n < 2)                              /* should never happen! */
# Line 64 | Line 64 | static double
64   est_DSFrad(const RBFNODE *rbf, const FVECT outvec)
65   {
66          const double    rad_epsilon = 0.03;
67 <        const double    DSFtarget = 0.60653066 * eval_rbfrep(rbf,outvec)
68 <                                                * output_orient*outvec[2];
67 >        const double    DSFtarget = 0.60653066 * eval_rbfrep(rbf,outvec) *
68 >                                                        COSF(outvec[2]);
69          double          inside_rad = rad_epsilon;
70          double          outside_rad = 0.5;
71          double          DSFinside = eval_DSFsurround(rbf, outvec, inside_rad);
72          double          DSFoutside = eval_DSFsurround(rbf, outvec, outside_rad);
73   #define interp_rad      inside_rad + (outside_rad-inside_rad) * \
74                                  (DSFtarget-DSFinside) / (DSFoutside-DSFinside)
75 <                                                /* interpolation search */
76 <        while (outside_rad-inside_rad > rad_epsilon) {
75 >                                                /* Newton's method (sort of) */
76 >        do {
77                  double  test_rad = interp_rad;
78 <                double  DSFtest = eval_DSFsurround(rbf, outvec, test_rad);
79 <                if (DSFtarget < DSFtest) {
78 >                double  DSFtest;
79 >                if (test_rad >= outside_rad)
80 >                        return(test_rad);
81 >                if (test_rad <= inside_rad)
82 >                        return(test_rad*(test_rad>0));
83 >                DSFtest = eval_DSFsurround(rbf, outvec, test_rad);
84 >                if (DSFtest > DSFtarget) {
85                          inside_rad = test_rad;
86                          DSFinside = DSFtest;
87                  } else {
88                          outside_rad = test_rad;
89                          DSFoutside = DSFtest;
90                  }
91 <        }
91 >                if (DSFoutside >= DSFinside)
92 >                        return(test_rad);
93 >        } while (outside_rad-inside_rad > rad_epsilon);
94          return(interp_rad);
95   #undef interp_rad
96   }
# Line 92 | Line 99 | est_DSFrad(const RBFNODE *rbf, const FVECT outvec)
99   static void
100   comp_bsdf_spec(void)
101   {
102 +        const double    max_hemi = 0.9;
103          double          peak_sum = 0;
104          double          rad_sum = 0;
105          int             n = 0;
# Line 113 | Line 121 | comp_bsdf_spec(void)
121          }
122          bsdf_spec_peak = peak_sum/(double)n;
123          bsdf_spec_rad = rad_sum/(double)n;
124 +        if ((2.*M_PI)*bsdf_spec_peak*bsdf_spec_rad*bsdf_spec_rad > max_hemi)
125 +                bsdf_spec_peak = max_hemi/((2.*M_PI)*bsdf_spec_rad*bsdf_spec_rad);
126   }
127  
128   /* Create a new migration holder (sharing memory for multiprocessing) */

Diff Legend

– Removed lines
+ Added lines
< Changed lines (old)
> Changed lines (new)