public static void AorThread(Mesh mesh, float[] squares, int offset, int count)
{
if (offset + count > mesh.Verteces.Count)
{
count = mesh.Verteces.Count - offset;
}
for (int i = offset; i < offset + count; i += 3)
{
var v0 = mesh.Verteces[i];
var v1 = mesh.Verteces[i + 1];
var v2 = mesh.Verteces[i + 2];
//var center1 = (v0.Position + v1.Position + v2.Position) / 3;
v0.Ao = v1.Ao = v2.Ao = 1.0f;
float res = 0;
for (int j = 0; j < mesh.Verteces.Count; j += 3)
{
if (j < i || j > i + 2)
{
var vv0 = mesh.Verteces[j];
var vv1 = mesh.Verteces[j + 1];
var vv2 = mesh.Verteces[j + 2];
//var center2 = (vv0.Position + vv1.Position + vv2.Position) / 3;
var v = vv0.Position - v0.Position;
var d2 = Vector3.Dot(v, v) + 1e-16;
v *= 1.0f / SstlHelper.Sqrt((float)d2);
//var rsq = (center1 - center2).LengthSquared;
//if (float.IsNaN(res))
//{
// continue;
//}
var value = ElementShadow(v, (float)d2, v0.Normal, vv0.Normal, squares[j]);
if (float.IsNaN(value))
{
value = 1;
}
res += value;
if (res >= 1)
{
break;
} // http://http.developer.nvidia.com/GPUGems2/gpugems2_chapter14.html
}
}
//res = SstlHelper.Clamp(res, 1, 0);
v0.Ao = v1.Ao = v2.Ao = 1.0f - res;
mesh.Verteces[i] = v0;
mesh.Verteces[i + 1] = v1;
mesh.Verteces[i + 2] = v2;
AorResult += 3;
}
}
static float ElementShadow(Vector3 v, float rSquared, Vector3 receiverNormal,
Vector3 emitterNormal, float emitterArea)
{
// we assume that emitterArea has already been divided by PI
return((1.0f - 1.0f / SstlHelper.Sqrt(emitterArea / rSquared + 1.0f)) *
SstlHelper.Clamp(Vector3.Dot(emitterNormal, v), 0f, 1.0f) *
SstlHelper.Clamp(4 * Vector3.Dot(receiverNormal, v), 0f, 1.0f));
}
