protected override void configureNewParticle(SSParticle p)
{
base.configureNewParticle(p);
float r = Interpolate.Lerp(radiusOffsetMin, radiusOffsetMax, nextFloat());
float theta = Interpolate.Lerp(thetaMin, thetaMax, nextFloat());
float phi = Interpolate.Lerp(phiMin, phiMax, nextFloat());
float xy = (float)Math.Cos(phi);
float x = xy * (float)Math.Cos(theta);
float y = xy * (float)Math.Sin(theta);
float z = (float)Math.Sin(phi);
Vector3 xAxis, yAxis;
OpenTKHelper.TwoPerpAxes(up, out xAxis, out yAxis);
Vector3 xyz = x * xAxis + y * yAxis + z * up;
p.pos = center + r * xyz;
float velocityMag = Interpolate.Lerp(velocityFromCenterMagnitudeMin, velocityFromCenterMagnitudeMax, nextFloat());
p.vel += velocityMag * xyz;
if (orientAwayFromCenter)
{
p.orientation.Z = -theta;
p.orientation.Y = phi;
}
//Console.WriteLine("particle emitted with vel = " + p.vel.Length);
}
protected void ComputeProjections(
List <SSObject> objects,
Matrix4 cameraView,
Matrix4 cameraProj,
float fov, float aspect, float cameraNearZ, float cameraFarZ)
{
if (m_light.GetType() != typeof(SSDirectionalLight))
{
throw new NotSupportedException();
}
SSDirectionalLight dirLight = (SSDirectionalLight)m_light;
// light-aligned unit vectors
Vector3 lightZ = dirLight.Direction.Normalized();
Vector3 lightX, lightY;
OpenTKHelper.TwoPerpAxes(lightZ, out lightX, out lightY);
// transform matrix from regular space into light aligned space
Matrix4 lightTransform = new Matrix4(
lightX.X, lightX.Y, lightX.Z, 0f,
lightY.X, lightY.Y, lightY.Z, 0f,
lightZ.X, lightZ.Y, lightZ.Z, 0f,
0f, 0f, 0f, 0f
);
// Step 0: camera projection matrix (nearZ and farZ modified) for each frustum split
float prevFarZ = cameraNearZ;
for (int i = 0; i < c_numberOfSplits; ++i)
{
// generate frustum splits using Practical Split Scheme (GPU Gems 3, 10.2.1)
float iRatio = (float)(i + 1) / (float)c_numberOfSplits;
float cLog = cameraNearZ * (float)Math.Pow(cameraFarZ / cameraNearZ, iRatio);
float cUni = cameraNearZ + (cameraFarZ - cameraNearZ) * iRatio;
float nextFarZ = LogVsLinearSplitFactor * cLog + (1f - LogVsLinearSplitFactor) * cUni;
float nextNearZ = prevFarZ;
// exported to the shader
m_viewSplits [i] = nextFarZ;
// create a view proj matrix with the nearZ, farZ values for the current split
m_frustumViewProjMatrices[i] = cameraView
* Matrix4.CreatePerspectiveFieldOfView(fov, aspect, nextNearZ, nextFarZ);
// create light-aligned AABBs of frustums
m_frustumLightBB [i] = SSAABB.FromFrustum(ref lightTransform, ref m_frustumViewProjMatrices [i]);
prevFarZ = nextFarZ;
}
#if true
// Optional scene-dependent optimization
for (int i = 0; i < c_numberOfSplits; ++i)
{
m_objsLightBB[i] = new SSAABB(float.PositiveInfinity, float.NegativeInfinity);
m_splitFrustums[i] = new SSFrustumCuller(ref m_frustumViewProjMatrices[i]);
m_shrink[i] = false;
}
foreach (var obj in objects)
{
// pass through all shadow casters and receivers
if (obj.renderState.toBeDeleted || obj.localBoundingSphereRadius <= 0f ||
!obj.renderState.visible || !obj.renderState.receivesShadows)
{
continue;
}
else
{
for (int i = 0; i < c_numberOfSplits; ++i)
{
if (m_splitFrustums[i].isSphereInsideFrustum(obj.worldBoundingSphere))
{
// determine AABB in light coordinates of the objects so far
m_shrink[i] = true;
Vector3 lightAlignedPos = Vector3.Transform(obj.worldBoundingSphereCenter, lightTransform);
Vector3 rad = new Vector3(obj.worldBoundingSphereRadius);
Vector3 localMin = lightAlignedPos - rad;
Vector3 localMax = lightAlignedPos + rad;
m_objsLightBB[i].UpdateMin(localMin);
m_objsLightBB[i].UpdateMax(localMax);
}
}
}
}
#endif
for (int i = 0; i < c_numberOfSplits; ++i)
{
if (m_shrink [i])
{
m_resultLightBB[i].Min = Vector3.ComponentMax(m_frustumLightBB [i].Min,
m_objsLightBB [i].Min);
m_resultLightBB [i].Max = Vector3.ComponentMin(m_frustumLightBB [i].Max,
m_objsLightBB [i].Max);
}
else
{
m_resultLightBB [i] = m_frustumLightBB [i];
}
}
for (int i = 0; i < c_numberOfSplits; ++i)
{
// Obtain view + projection + crop matrix, need it later
Matrix4 shadowView, shadowProj;
viewProjFromLightAlignedBB(ref m_resultLightBB [i], ref lightTransform, ref lightY,
out shadowView, out shadowProj);
m_shadowViewProjMatrices[i] = shadowView * shadowProj * c_cropMatrices[i];
// obtain view + projection + clio + bias
m_shadowViewProjBiasMatrices[i] = m_shadowViewProjMatrices[i] * c_biasMatrix;
// There is, currently, no mathematically derived solution to how much Poisson spread scaling
// you need for each split. Current improvisation combines 1) increasing spread for the near
// splits; reducing spread for the far splits and 2) reducing spread for splits with larger
// light-aligned areas; increasing spread for splits with smaller light-aligned areas
m_poissonScaling [i] = m_resultLightBB [i].Diff().Xy / (100f * (float)Math.Pow(3.0, i - 1));
}
// Combine all splits' BB into one and extend it to include shadow casters closer to light
SSAABB castersLightBB = new SSAABB(float.PositiveInfinity, float.NegativeInfinity);
for (int i = 0; i < c_numberOfSplits; ++i)
{
castersLightBB.Combine(ref m_resultLightBB [i]);
}
// extend Z of the AABB to cover shadow-casters closer to the light
foreach (var obj in objects)
{
if (obj.renderState.toBeDeleted || obj.localBoundingSphereRadius <= 0f ||
!obj.renderState.visible || !obj.renderState.castsShadow)
{
continue;
}
else
{
Vector3 lightAlignedPos = Vector3.Transform(obj.worldBoundingSphereCenter, lightTransform);
Vector3 rad = new Vector3(obj.worldBoundingSphereRadius);
Vector3 localMin = lightAlignedPos - rad;
Vector3 localMax = lightAlignedPos + rad;
if (localMin.Z < castersLightBB.Min.Z)
{
if (OpenTKHelper.RectsOverlap(castersLightBB.Min.Xy,
castersLightBB.Max.Xy,
localMin.Xy,
localMax.Xy))
{
castersLightBB.Min.Z = localMin.Z;
}
}
}
}
// Generate frustum culler from the BB extended towards light to include shadow casters
Matrix4 frustumView, frustumProj;
viewProjFromLightAlignedBB(ref castersLightBB, ref lightTransform, ref lightY,
out frustumView, out frustumProj);
Matrix4 frustumMatrix = frustumView * frustumProj;
FrustumCuller = new SSFrustumCuller(ref frustumMatrix);
}
private void ComputeProjections(List <SSObject> objects,
SSLightBase light,
Matrix4 cameraView, Matrix4 cameraProj)
{
if (light.GetType() != typeof(SSDirectionalLight))
{
throw new NotSupportedException();
}
SSDirectionalLight dirLight = (SSDirectionalLight)light;
// light-aligned unit vectors
Vector3 lightZ = dirLight.Direction.Normalized();
Vector3 lightX, lightY;
OpenTKHelper.TwoPerpAxes(lightZ, out lightX, out lightY);
// transform matrix from regular space into light aligned space
Matrix4 lightTransform = new Matrix4(
lightX.X, lightX.Y, lightX.Z, 0f,
lightY.X, lightY.Y, lightY.Z, 0f,
lightZ.X, lightZ.Y, lightZ.Z, 0f,
0f, 0f, 0f, 0f
);
// Find AABB of frustum corners in light coordinates
Matrix4 cameraViewProj = cameraView * cameraProj;
SSAABB frustumLightBB = SSAABB.FromFrustum(ref lightTransform, ref cameraViewProj);
bool shrink = false;
SSAABB objsLightBB = new SSAABB(float.PositiveInfinity, float.NegativeInfinity);
#if true
// (optional) scene dependent optimization
// Trim the light-bounding box by the shadow receivers (only in light-space x,y,maxz)
SSFrustumCuller cameraFrustum = new SSFrustumCuller(ref cameraViewProj);
foreach (var obj in objects)
{
// pass through all shadow casters and receivers
if (obj.renderState.toBeDeleted || !obj.renderState.visible ||
!obj.renderState.receivesShadows || obj.localBoundingSphereRadius <= 0f)
{
continue;
}
else if (cameraFrustum.isSphereInsideFrustum(obj.worldBoundingSphere))
{
// determine AABB in light coordinates of the objects so far
shrink = true;
Vector3 lightAlignedPos = Vector3.Transform(obj.worldBoundingSphereCenter, lightTransform);
Vector3 rad = new Vector3(obj.worldBoundingSphereRadius);
Vector3 localMin = lightAlignedPos - rad;
Vector3 localMax = lightAlignedPos + rad;
objsLightBB.UpdateMin(localMin);
objsLightBB.UpdateMax(localMax);
}
}
#endif
// Optimize the light-frustum-projection bounding box by the object-bounding-box
SSAABB resultLightBB = new SSAABB(float.PositiveInfinity, float.NegativeInfinity);
if (shrink)
{
// shrink the XY & far-Z coordinates..
resultLightBB.Min.Xy = Vector2.ComponentMax(frustumLightBB.Min.Xy, objsLightBB.Min.Xy);
resultLightBB.Min.Z = objsLightBB.Min.Z;
resultLightBB.Max = Vector3.ComponentMin(frustumLightBB.Max, objsLightBB.Max);
}
else
{
resultLightBB = frustumLightBB;
}
// View and projection matrices, used by the scene later
viewProjFromLightAlignedBB(ref resultLightBB, ref lightTransform, ref lightY,
out m_shadowViewMatrix, out m_shadowProjMatrix);
// Now extend Z of the result AABB to cover shadow-casters closer to the light inside the
// original box
foreach (var obj in objects)
{
if (obj.renderState.toBeDeleted || !obj.renderState.visible ||
!obj.renderState.castsShadow || obj.localBoundingSphereRadius <= 0f)
{
continue;
}
Vector3 lightAlignedPos = Vector3.Transform(obj.worldBoundingSphereCenter, lightTransform);
Vector3 rad = new Vector3(obj.worldBoundingSphereRadius);
Vector3 localMin = lightAlignedPos - rad;
if (localMin.Z < resultLightBB.Min.Z)
{
Vector3 localMax = lightAlignedPos + rad;
if (OpenTKHelper.RectsOverlap(resultLightBB.Min.Xy,
resultLightBB.Max.Xy,
localMin.Xy,
localMax.Xy))
{
resultLightBB.Min.Z = localMin.Z;
}
}
}
// Generate frustum culler from the BB extended towards light to include shadow casters
Matrix4 frustumView, frustumProj;
viewProjFromLightAlignedBB(ref resultLightBB, ref lightTransform, ref lightY,
out frustumView, out frustumProj);
Matrix4 frustumMatrix = frustumView * frustumProj;
FrustumCuller = new SSFrustumCuller(ref frustumMatrix);
}