101     int j, co;  in SetLutTable()  local
106     for (co = 0; co < 3; co++)  in SetLutTable()
107         memset(qlut[co], 0, sizeof(qlut[0]));  in SetLutTable()
132     for (co = 0; co < 3; co++)  in SetLutTable()
134         sigma_t_prim[co] = sigma_s_prim[co] + sigma_a[co];  in SetLutTable()
135         mat.zr[co] = geom_scale / sigma_t_prim[co];  in SetLutTable()
136         mat.albedo[co] = sigma_s_prim[co] / sigma_t_prim[co];  in SetLutTable()
137         mat.beta[co] = sqrt(3.0f * (1.0f - mat.albedo[co]));  in SetLutTable()
138         mat.beta_prime[co] = mat.beta[co] * mat.Bi;  in SetLutTable()
139         mat.dif_refl[co] = 0.5f * mat.albedo[co] * (exp(-mat.beta[co]) + exp(-mat.beta_prime[co]));  in SetLutTable()
140         mat.i0[co] = exp(-mat.beta[co]) * sqrt(2.0f * DMP_PI / mat.beta[co]) +  in SetLutTable()
141           mat.Bi * exp(-mat.beta_prime[co]) * sqrt(2.0f * DMP_PI / mat.beta_prime[co]);  in SetLutTable()
142         mat.i1[co] = DMP_PI * (exp(-mat.beta[co]) + exp(-mat.beta_prime[co])) / mat.i0[co];  in SetLutTable()
150         for (co = 0; co < 3; co++)  in SetLutTable()
154                 qlut[co][j] = mat.dif_refl[co] * LN;  in SetLutTable()
155             gamma = mat.zr[co] * sqrt(kappa);  in SetLutTable()
157 …qlut[co][j] += mat.i0[co] * (1.0f / (1.0f + mat.i1[co] * h)) * gamma * mat.albedo[co] * 0.25f * RE…  in SetLutTable()
158             qlut[co][j] = pow(qlut[co][j], display_gamma);  in SetLutTable()
166         for (co = 0; co < 3; co++)  in SetLutTable()
169                 qlut[co][j] = mat.dif_refl[co] * LN;  in SetLutTable()
170             gamma = mat.zr[co] * sqrt(kappa);  in SetLutTable()
172             qlut[co][j] += mat.i0[co] * (1.0f / (1.0f + mat.i1[co] * h)) *  in SetLutTable()
173                 gamma * mat.albedo[co] * 0.25f * REV_PI;  in SetLutTable()
174             qlut[co][j] = pow(qlut[co][j], display_gamma);  in SetLutTable()
179         for (co = 0; co < 3; co++)  in SetLutTable()
180             qlut[co][j + 256] = qlut[co][j + 1] - qlut[co][j];  in SetLutTable()
182     for (co = 0; co < 3; co++)  in SetLutTable()
183         qlut[co][127 + 256] = pow(mat.dif_refl[co], display_gamma) - qlut[co][127];  in SetLutTable()
186         for (co = 0; co < 3; co++)  in SetLutTable()
187             qlut[co][j + 256] = qlut[co][j + 1] - qlut[co][j];  in SetLutTable()
189     for (co = 0; co < 3; co++)  in SetLutTable()
190         qlut[co][255 + 256] = qlut[co][0] - qlut[co][255];  in SetLutTable()