98     } mat;  in SetLutTable()  local
128     mat.F_dr = - 1.440f / (refi * refi) + 0.71f / refi + 0.668f + 0.0636f * refi;  in SetLutTable()
129     A = (1.0f + mat.F_dr) / (1.0f - mat.F_dr);  in SetLutTable()
130     mat.Bi = 1.0f + (4.0f / 3.0f) * A;  in SetLutTable()
131     mat.refi = refi;  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()
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()
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()
183         qlut[co][127 + 256] = pow(mat.dif_refl[co], display_gamma) - qlut[co][127];  in SetLutTable()