Differential D16655 Diff 58189 intern/cycles/kernel/integrator/mnee.h

Changeset View

Standalone View

intern/cycles/kernel/integrator/mnee.h

Show First 20 Lines • Show All 102 Lines • ▼ Show 20 Lines
{		{
float det = mat22_determinant(m);		float det = mat22_determinant(m);
if (fabsf(det) < MNEE_MIN_DETERMINANT)		if (fabsf(det) < MNEE_MIN_DETERMINANT)
return 0.f;		return 0.f;
m_inverse = make_float4(m.w, -m.y, -m.z, m.x) / det;		m_inverse = make_float4(m.w, -m.y, -m.z, m.x) / det;
return det;		return det;
}		}

/* Update light sample */
ccl_device_forceinline void mnee_update_light_sample(KernelGlobals kg,
const float3 P,
ccl_private LightSample *ls)
{
/* correct light sample position/direction and pdf
* NOTE: preserve pdf in area measure */
const ccl_global KernelLight *klight = &kernel_data_fetch(lights, ls->lamp);

if (ls->type == LIGHT_POINT \|\| ls->type == LIGHT_SPOT) {
ls->D = normalize_len(ls->P - P, &ls->t);
ls->Ng = -ls->D;

float2 uv = map_to_sphere(ls->Ng);
ls->u = uv.x;
ls->v = uv.y;

float invarea = klight->spot.invarea;
ls->eval_fac = (0.25f * M_1_PI_F) * invarea;
ls->pdf = invarea;

if (ls->type == LIGHT_SPOT) {
/* spot light attenuation */
float3 dir = make_float3(klight->spot.dir[0], klight->spot.dir[1], klight->spot.dir[2]);
ls->eval_fac *= spot_light_attenuation(
dir, klight->spot.spot_angle, klight->spot.spot_smooth, ls->Ng);
}
}
else if (ls->type == LIGHT_AREA) {
float invarea = fabsf(klight->area.invarea);
ls->D = normalize_len(ls->P - P, &ls->t);
ls->pdf = invarea;
if (klight->area.tan_spread > 0.f) {
ls->eval_fac = 0.25f * invarea;
ls->eval_fac *= light_spread_attenuation(
ls->D, ls->Ng, klight->area.tan_spread, klight->area.normalize_spread);
}
}

ls->pdf *= kernel_data.integrator.pdf_lights;
}

/* Manifold vertex setup from ray and intersection data */		/* Manifold vertex setup from ray and intersection data */
ccl_device_forceinline void mnee_setup_manifold_vertex(KernelGlobals kg,		ccl_device_forceinline void mnee_setup_manifold_vertex(KernelGlobals kg,
ccl_private ManifoldVertex *vtx,		ccl_private ManifoldVertex *vtx,
ccl_private ShaderClosure *bsdf,		ccl_private ShaderClosure *bsdf,
const float eta,		const float eta,
const float2 n_offset,		const float2 n_offset,
ccl_private const Ray *ray,		ccl_private const Ray *ray,
ccl_private const Intersection *isect,		ccl_private const Intersection *isect,
▲ Show 20 Lines • Show All 653 Lines • ▼ Show 20 Lines	ccl_device_forceinline bool mnee_path_contribution(KernelGlobals kg,
float wo_len;		float wo_len;
float3 wo = normalize_len(vertices[0].p - sd->P, &wo_len);		float3 wo = normalize_len(vertices[0].p - sd->P, &wo_len);

/* Initialize throughput and evaluate receiver bsdf * \|n.wo\|. */		/* Initialize throughput and evaluate receiver bsdf * \|n.wo\|. */
surface_shader_bsdf_eval(kg, state, sd, wo, throughput, ls->shader);		surface_shader_bsdf_eval(kg, state, sd, wo, throughput, ls->shader);

/* Update light sample with new position / direct.ion		/* Update light sample with new position / direct.ion
* and keep pdf in vertex area measure */		* and keep pdf in vertex area measure */
mnee_update_light_sample(kg, vertices[vertex_count - 1].p, ls);		light_sample_update_position(kg, ls, vertices[vertex_count - 1].p);

/* Save state path bounce info in case a light path node is used in the refractive interface or		/* Save state path bounce info in case a light path node is used in the refractive interface or
* light shader graph. */		* light shader graph. */
const int transmission_bounce = INTEGRATOR_STATE(state, path, transmission_bounce);		const int transmission_bounce = INTEGRATOR_STATE(state, path, transmission_bounce);
const int diffuse_bounce = INTEGRATOR_STATE(state, path, diffuse_bounce);		const int diffuse_bounce = INTEGRATOR_STATE(state, path, diffuse_bounce);
const int bounce = INTEGRATOR_STATE(state, path, bounce);		const int bounce = INTEGRATOR_STATE(state, path, bounce);

/* Set diffuse bounce info . */		/* Set diffuse bounce info . */
▲ Show 20 Lines • Show All 278 Lines • Show Last 20 Lines