1 | ///////////////////////////////////////////////////////////////////////////////
|
---|
2 | //
|
---|
3 | // ## ######
|
---|
4 | // ###### ###
|
---|
5 | // ## ############### Shark 3D Engine (www.shark3d.com)
|
---|
6 | // ########## # # #
|
---|
7 | // ######## Copyright (c) 1996-2006 Spinor GmbH.
|
---|
8 | // ######### # # # All rights reserved.
|
---|
9 | // ## ##########
|
---|
10 | // ##
|
---|
11 | //
|
---|
12 | ///////////////////////////////////////////////////////////////////////////////
|
---|
13 |
|
---|
14 | //@cpp
|
---|
15 |
|
---|
16 | #include "gtp_shader_opticsutil.h"
|
---|
17 | #include "../../../util/math/util_geo_types.h"
|
---|
18 |
|
---|
19 | ///////////////////////////////////////////////////////////////////////////////
|
---|
20 |
|
---|
21 | bool gtp_ShaderOpticsUtil::CalcViewNormalDist(
|
---|
22 | s3d_CUtilVec3f &ViewNormal, float &ViewDist,
|
---|
23 | s3d_CUtilTranslQuatf_cr ViewTransf)
|
---|
24 | {
|
---|
25 | if(ViewTransf.m_Quat == 0)
|
---|
26 | return false;
|
---|
27 |
|
---|
28 | // Clip plane is local z-plane.
|
---|
29 | s3d_CUtilVec3f LocalNormal = s3d_CUtilVec3f(0, 0, 1);
|
---|
30 |
|
---|
31 | // Clip plane in view coordinates. Plane equation is:
|
---|
32 | // Point * ViewNormal + ViewDist = 0.
|
---|
33 | s3d_CUtilMat3x3f TransfMat = ViewTransf.m_Quat.GetMat();
|
---|
34 | ViewNormal = (TransfMat * LocalNormal).GetNormalized();
|
---|
35 | ViewDist = - ViewNormal * ViewTransf.m_Transl;
|
---|
36 |
|
---|
37 | // Front side visible?:
|
---|
38 | return ViewDist > 0;
|
---|
39 | }
|
---|
40 |
|
---|
41 | bool gtp_ShaderOpticsUtil::CalcChildCamFull(
|
---|
42 | s3d_CEngGfxCam &ChildCam,
|
---|
43 | void *Exclude,
|
---|
44 | int MaxRecursion, float Coeff, float Precision,
|
---|
45 | bool UseCompactFrustum, s3d_CEngShaderScope *Scope,
|
---|
46 | s3d_CUtilTranslQuatf_cr ViewTransf,
|
---|
47 | s3d_CUtilVec3f_cr ViewNormal, float ViewDist)
|
---|
48 | {
|
---|
49 | s3d_CEngGfxPortion *Portion = Scope->m_Portion;
|
---|
50 |
|
---|
51 | // Check recursion level:
|
---|
52 | if(Portion->m_Cam.m_Recursion >= MaxRecursion)
|
---|
53 | return false;
|
---|
54 |
|
---|
55 | s3d_CUtilTranslQuatf RecTransf = Portion->m_Cam.m_Transf;
|
---|
56 | bool RecReverse = Portion->m_Cam.m_Reverse;
|
---|
57 |
|
---|
58 | s3d_CUtilMat4x4f RecProj = Portion->m_Cam.m_Proj;
|
---|
59 |
|
---|
60 | s3d_CUtilVec4f RecClipPlane;
|
---|
61 | RecClipPlane.m_x = - ViewNormal.m_x;
|
---|
62 | RecClipPlane.m_y = - ViewNormal.m_y;
|
---|
63 | RecClipPlane.m_z = - ViewNormal.m_z;
|
---|
64 | RecClipPlane.m_w = - ViewDist;
|
---|
65 |
|
---|
66 | if(Coeff != 0 && Coeff != 1)
|
---|
67 | {
|
---|
68 | // Calculate extra transformation in view coordinates:
|
---|
69 |
|
---|
70 | float Change = Coeff - 1;
|
---|
71 | // (Assuming that normal vector is normalized:)
|
---|
72 | // NewPoint
|
---|
73 | // = Point + Point-To-Plane-Distance * (ViewNormal * Change)
|
---|
74 | // = Point + (Point * ViewNormal + ViewDist) * (ViewNormal * Change)
|
---|
75 | // = (I + ViewNormal (x) ViewNormal * Change) * Point
|
---|
76 | // + ViewNormal * (ViewDist * Change)
|
---|
77 | // =: ExtraMat * Point + ExtraTransl
|
---|
78 | s3d_CUtilMat3x3f ExtraMat(s3d_CUtilMat3x3f(1)
|
---|
79 | + s3d_CUtilMat3x3f::ByTensorProd(
|
---|
80 | ViewNormal, ViewNormal * Change));
|
---|
81 | s3d_CUtilVec3f ExtraTransl(ViewNormal * (ViewDist * Change));
|
---|
82 |
|
---|
83 | // PartQuat maps the view z into the reflected view z:
|
---|
84 | s3d_CUtilVec3f BeforeDirZ(0, 0, 1);
|
---|
85 | s3d_CUtilVec3f AppliedDirZ = ExtraMat.GetColZ().GetNormalized();
|
---|
86 | s3d_CUtilVec3f DirY(0, 1, 0);
|
---|
87 | s3d_CUtilQuatf PartQuat = s3d_CUtilQuatf::ByPathAlso(
|
---|
88 | AppliedDirZ, BeforeDirZ, DirY, DirY);
|
---|
89 | // (ExtraTransl, ExtraMat)
|
---|
90 | // := (0, RestMat) x (PartTransl, PartQuat)
|
---|
91 | // = (RestMat * PartTransl, RestMat * PartQuat)
|
---|
92 | s3d_CUtilMat3x3f RestMat = ExtraMat * PartQuat.GetInvMat();
|
---|
93 | s3d_CUtilVec3f PartTransl = RestMat.GetInv() * ExtraTransl;
|
---|
94 | s3d_CUtilTranslQuatf PartTransf(PartTransl, PartQuat);
|
---|
95 |
|
---|
96 | // Apply:
|
---|
97 | RecTransf = RecTransf.ConcatInv(PartTransf);
|
---|
98 | RecProj = s3d_CUtilGeof::ConcatMapTranslMat3(RecProj, 0, RestMat);
|
---|
99 | RecClipPlane = RecClipPlane
|
---|
100 | * s3d_CUtilGeof::GetMapOfInvTransf3(PartTransf);
|
---|
101 | if(Coeff < 0)
|
---|
102 | {
|
---|
103 | RecClipPlane = - RecClipPlane;
|
---|
104 | RecReverse = !RecReverse;
|
---|
105 | }
|
---|
106 | }
|
---|
107 |
|
---|
108 | if(UseCompactFrustum)
|
---|
109 | {
|
---|
110 | // Use the near clipping plane as mirror clipping plane.
|
---|
111 | //
|
---|
112 | // Clipping plane condition:
|
---|
113 | // -> (NewProj * P).z >= - (NewProj * P).w
|
---|
114 | // -> (NewProj.z? + NewProj.w?) * P >= 0
|
---|
115 | // Only modify Proj.z?
|
---|
116 | // -> (NewProj.z? + Proj.w?) * P >= 0
|
---|
117 | //
|
---|
118 | // Plane:
|
---|
119 | // -> Plane * P >= 0
|
---|
120 | // -> NewProj.z? + NewProj.w? = Plane * Fac
|
---|
121 | // -> NewProj.z? = Plane * Fac - Proj.w?
|
---|
122 | //
|
---|
123 | // Choose far clipping plane to include FarPoint
|
---|
124 | // (typically choose plane dir infinity):
|
---|
125 | // NewProj * FarPoint = ~ (0, 0, 1, 1)
|
---|
126 | // -> NewProj.z? * FarPoint = NewProj.w? * FarPoint
|
---|
127 | // -> [Plane * Fac - Proj.w?] * FarPoint = Proj.w? * FarPoint
|
---|
128 | // -> Fac = (Proj.w? * FarPoint * 2) / (Plane * FarPoint)
|
---|
129 | //
|
---|
130 | // Choose FarPoint at infinity in that way
|
---|
131 | // that the complete old frustum is included.
|
---|
132 | // For this, select one frustum vertex at infinity
|
---|
133 | // depending on the plane slope signs.
|
---|
134 |
|
---|
135 | s3d_CUtilVec4f FarPoint(
|
---|
136 | s3d_SysSign(RecClipPlane.m_x),
|
---|
137 | s3d_SysSign(RecClipPlane.m_y),
|
---|
138 | 1,
|
---|
139 | 0);
|
---|
140 |
|
---|
141 | s3d_CUtilVec4f RecProjRowW = RecProj.GetRowW();
|
---|
142 | float PlaneFac = ((RecProjRowW * FarPoint) * 2)
|
---|
143 | / (RecClipPlane * FarPoint);
|
---|
144 | RecProj.SetRowZ(RecClipPlane * PlaneFac - RecProjRowW);
|
---|
145 | RecClipPlane = 0;
|
---|
146 | }
|
---|
147 |
|
---|
148 | s3d_CUtilVec3f WorldTransl
|
---|
149 | = Portion->m_Cam.m_Transf.ApplyToPoint(ViewTransf.m_Transl);
|
---|
150 |
|
---|
151 | // New camera structure:
|
---|
152 | ChildCam = Portion->m_Cam;
|
---|
153 | ChildCam.m_Exclude = Exclude;
|
---|
154 | ChildCam.m_Recursion = Portion->m_Cam.m_Recursion + 1;
|
---|
155 | ChildCam.m_VisibPoint = WorldTransl;
|
---|
156 | ChildCam.m_Transf = RecTransf;
|
---|
157 | ChildCam.m_Proj = RecProj;
|
---|
158 | ChildCam.m_Reverse = RecReverse;
|
---|
159 |
|
---|
160 | if(Precision > 0)
|
---|
161 | {
|
---|
162 | ChildCam.m_CullPix *= Precision;
|
---|
163 | ChildCam.m_MeshPix *= Precision;
|
---|
164 | ChildCam.m_VisibMeshPix *= Precision;
|
---|
165 | }
|
---|
166 |
|
---|
167 | return true;
|
---|
168 | }
|
---|
169 |
|
---|
170 | bool gtp_ShaderOpticsUtil::CalcChildRange(
|
---|
171 | s3d_CUtilVec3f &SubNeg, s3d_CUtilVec3f &SubPos,
|
---|
172 | s3d_CUtilTranslQuatf_cr ViewTransf,
|
---|
173 | s3d_CUtilVec3f_cr BoundCen, s3d_CUtilVec3f_cr BoundExt,
|
---|
174 | s3d_CEngShaderScope *Scope)
|
---|
175 | {
|
---|
176 | s3d_CEngGfxPortion *Portion = Scope->m_Portion;
|
---|
177 |
|
---|
178 | // Calculate bounding rectangle in screen coordinates:
|
---|
179 | s3d_CUtilMat3x3f TransfMat = ViewTransf.m_Quat.GetMat();
|
---|
180 | s3d_CUtilVec3f ViewBoundCen = ViewTransf.m_Transl + TransfMat * BoundCen;
|
---|
181 | s3d_CUtilMat3x3f ViewBoundExt = TransfMat.ColGetScaled(BoundExt);
|
---|
182 |
|
---|
183 | s3d_CUtilMat4x4f TotalBoundMap =
|
---|
184 | s3d_CUtilGeof::ConcatMapTranslMat3(
|
---|
185 | Portion->m_Cam.m_Proj, ViewBoundCen, ViewBoundExt);
|
---|
186 | s3d_CUtilVec3f MapMin, MapMax;
|
---|
187 | s3d_CUtilGeof::CalcRangeMapMinusPlusOne3(
|
---|
188 | MapMin, MapMax, -1, +1, TotalBoundMap);
|
---|
189 | MapMin.m_z = -1; // Keep near plane to avoid problems with linked cam.
|
---|
190 | MapMax.m_z = 1; // Look to infinity.
|
---|
191 | if(MapMin.IsAnyCompGe(MapMax))
|
---|
192 | return false;
|
---|
193 |
|
---|
194 | // Calculate new rendering target rectangle:
|
---|
195 | SubNeg = (MapMin + s3d_CUtilVec3f(1)) * 0.5f;
|
---|
196 | SubPos = (MapMax + s3d_CUtilVec3f(1)) * 0.5f;
|
---|
197 |
|
---|
198 | return true;
|
---|
199 | }
|
---|
200 |
|
---|
201 | void gtp_ShaderOpticsUtil::SnapChildRange(
|
---|
202 | s3d_CUtilVec3f &SubNeg, s3d_CUtilVec3f &SubPos,
|
---|
203 | int SnapWidth, int SnapHeight)
|
---|
204 | {
|
---|
205 | S3D_SYS_ASSERT(SnapWidth > 0 && SnapHeight > 0);
|
---|
206 | float RcpSnapWidth = 1.0f / SnapWidth;
|
---|
207 | float RcpSnapHeight = 1.0f / SnapHeight;
|
---|
208 | SubNeg.m_x = s3d_SysRoundFloor(SubNeg.m_x * SnapWidth) * RcpSnapWidth;
|
---|
209 | SubPos.m_x = s3d_SysRoundCeil(SubPos.m_x * SnapWidth) * RcpSnapWidth;
|
---|
210 | SubNeg.m_y = s3d_SysRoundFloor(SubNeg.m_y * SnapHeight) * RcpSnapHeight;
|
---|
211 | SubPos.m_y = s3d_SysRoundCeil(SubPos.m_y * SnapHeight) * RcpSnapHeight;
|
---|
212 | }
|
---|
213 |
|
---|
214 | void gtp_ShaderOpticsUtil::ApplyRangeToChildCam(
|
---|
215 | s3d_CEngGfxCam &ChildCam, s3d_CUtilMat4x4f &BaseProj,
|
---|
216 | s3d_CUtilVec3f_cr SubNeg, s3d_CUtilVec3f_cr SubPos)
|
---|
217 | {
|
---|
218 | s3d_CUtilVec3f SubDelta = SubPos - SubNeg;
|
---|
219 | s3d_CUtilVec3f SubSum = SubPos + SubNeg;
|
---|
220 |
|
---|
221 | // Calculate new target matrix:
|
---|
222 | s3d_CUtilVec2f SubNegXY(SubNeg.m_x, SubNeg.m_y);
|
---|
223 | s3d_CUtilVec2f SubDeltaXY(SubDelta.m_x, SubDelta.m_y);
|
---|
224 | ChildCam.m_ViewPane.m_Transl
|
---|
225 | = ChildCam.m_ViewPane.m_Transl
|
---|
226 | + SubNegXY.CompMul(ChildCam.m_ViewPane.m_Scale);
|
---|
227 | ChildCam.m_ViewPane.m_Scale
|
---|
228 | = SubDeltaXY.CompMul(ChildCam.m_ViewPane.m_Scale);
|
---|
229 | float ViewDepthDelta = ChildCam.m_ViewDepthEnd - ChildCam.m_ViewDepthStart;
|
---|
230 | ChildCam.m_ViewDepthStart
|
---|
231 | = ChildCam.m_ViewDepthStart + ViewDepthDelta * SubNeg.m_z;
|
---|
232 |
|
---|
233 | // Calculate new projection matrix:
|
---|
234 | s3d_CUtilVec3f ScrScal = s3d_CUtilVec3f(1).CompDiv(SubDelta);
|
---|
235 | s3d_CUtilVec3f ScrOffs = (s3d_CUtilVec3f(1) - SubSum).CompMul(ScrScal);
|
---|
236 | s3d_CUtilMat4x4f ExtraProj = 1;
|
---|
237 | ExtraProj.m_xx = ScrScal.m_x;
|
---|
238 | ExtraProj.m_yy = ScrScal.m_y;
|
---|
239 | ExtraProj.m_zz = ScrScal.m_z;
|
---|
240 | ExtraProj.m_xw = ScrOffs.m_x;
|
---|
241 | ExtraProj.m_yw = ScrOffs.m_y;
|
---|
242 | ExtraProj.m_zw = ScrOffs.m_z;
|
---|
243 | ChildCam.m_Proj = ExtraProj * ChildCam.m_Proj;
|
---|
244 | BaseProj = ExtraProj * BaseProj;
|
---|
245 | }
|
---|
246 |
|
---|
247 | int gtp_ShaderOpticsUtil::RoundToPowerOfTwoDown(int x)
|
---|
248 | {
|
---|
249 | x |= (x >> 16);
|
---|
250 | x |= (x >> 8);
|
---|
251 | x |= (x >> 4);
|
---|
252 | x |= (x >> 2);
|
---|
253 | x |= (x >> 1);
|
---|
254 | return x & ~(x >> 1);
|
---|
255 | }
|
---|
256 |
|
---|
257 | int gtp_ShaderOpticsUtil::RoundToPowerOfTwoUp(int x)
|
---|
258 | {
|
---|
259 | return RoundToPowerOfTwoDown((x << 1) - 1);
|
---|
260 | }
|
---|
261 |
|
---|
262 | ///////////////////////////////////////////////////////////////////////////////
|
---|