nVIDIA SDK White Paper
Vertex Texture Fetch Water
作者:
Jeremy Zelsnack
本文版权归原作者所有,仅供个人学习使用,请勿转载,勿用于任何商业用途。
由于本人水平有限,难免出错,不清楚的地方请大家以原著为准。也欢迎大家和我多多交流。
Blog:http://blog.csdn.net/soilwork
水面模拟
水面模拟的目标是在合理的计算量下,创造视觉上真是可信的水面。其中可能的方法之一是基于流体力学中的Navier-Stokes方程求解。这种方法可以创建极度真是的水面,不幸的是,计算量太过沉重。一种较为简单,同时资源消耗较少的方法是在一张均匀网格上解2D波方程。牺牲一点点真实性换取实现的简洁和速度。本文描述了2D波方程技术。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq1.jpg)
上面的就是2D波方程。这个方程表示点将在Y方向上上下起伏。C表示波传播的速度。直观的说,这个方程表示点上下起伏的加速度,和表面起伏改变的快慢是成正比的。
看到这里,你可能会问,“如何来求解这样一个方程呢?”GPU并不知道偏导数是什么;只知道空间位置而已。如果我们把方程左边对t进行两次积分,就得到了y值。如果继续把计算出的y值与点x和z的值组合到一起,就得到了3D坐标,这就是GPU所能理解的了。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq2.jpg)
接下来的难题就是如何对方程的右边进行积分。有许多种积分的技术可以使用,每种都有其优点和缺点。对我的程序来说,需要的是稳定性和速度。稳定性是相当重要的,因为公正的说,这个方程太过于呆板。方程右边的C
2将使左边的值迅速增大。因此,不能使用类似欧拉的简单积分方法――它们不太稳定。Runge-Kutta方法也是不行的,因为它需要计算中间值,同时,还对储存速率有很高要求。这对于我们的程序来说太复杂,太耗费资源。Verlet积分计算量很轻,同时也很稳定,对储存速率也没有要求。不幸的是我忘记了Verlet方程(-_-# 汗~~),并且在30秒内都没有找到任何关于Verlet方程的参考资料(=。=)。因此,我重新发明了一种类似于Vertet的方法。
如果我们知道了点在t=0时刻的位置P(t0),以及t=1时刻的位置P(t1),加速度a(t1),那么使用方程3,只需要一些高中的物理知识,就可以求解t=2时刻的位置P(t2)。这里的推导假设采样步长是常数,因此,t1-t0 = t2 – t1 = t(n+1) – tn。当然,也可以步长值也可以为非常量,但这就需要更多的数学计算。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq3.jpg)
方程3和真正的Verlet积分(方程4)相当类似。唯一的不同,就是在加速度系数上有区别。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq4.jpg)
好了,方程3描述了一种计算积分的方法。但是有一个问题,它并不像看起来那么稳定,而是抖动的(jittery)。我们需要做一些额外的处理,来消耗由于不精确的积分而带来的额外能量。可以添加一个常量因子来减小通过推测而得来的速率。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq5.jpg)
好了,现在我们有了积分的技术,几乎可以对2D波方程的右边进行积分了。现在来看看如何计算方程右边的偏导数。方法之一,就是找一个和我们所要模拟的顶点相称的表面。有了这个表面,我们就能分析推测偏导数的计算公式。
那么什么样的表面才和要模拟的网格类似呢?我们需要一个快速,简单同时能在一定程度上保持精确的模型。3次样条曲线集合正好能满足我们的要求。可以对问题进行简化,变为对两个2D三次样条曲线的问题。这样就推导出了公式6。对一个三次方程来说,需要指定4个常数。其中,3个显而易见的常数就是当前以及邻近顶点的坐标。第四个常量则是当前顶点位置下的导数。图一描述了这种情况。假设切线和邻近的顶点是平行的,同时,网格也是均匀的,则可以得到一个2阶偏导数方程(公式7)。使用偏导数,右边的2D波公式就变为了方程8。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq6.jpg)
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/eq7.jpg)
好了,现在我们简单分析了计算方程右边的算法,以及积分方法,可以实际计算网格中顶点的y值了。由于我们假设网格中的所有顶点之间的距离都为1个单位(在世界坐标中不一定合适),因此需要调整波的速度,让他看起来正常。
好了,现在可以使用简单的迭带来求解波形方程,我们让GPU来完成这个任务吧。我们把均匀网格的高度保存为一张纹理。同时,还需要追踪上一张高度图的值。这样就可以计算下一张高度图的值了。下面是使用HLSL求解波形方程的代码:
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/HLSL.jpg)
顶点纹理将保存为16bit的浮点纹理。使用8bit的整数是不可行的,16bit的整数虽然可以满足精度的要求,但无论如何浮点数都是最好以及最方便的选择。示例的程序中,我们把高度和法线混合起来,保存为一张D3DFMT_A16B16G16R16的纹理。这样做,可以减少vertex shader中纹理拾取的次数。
需要注意的是,GeForce 6系列的硬件在vertex shader中并不支持D3DFMT_A16B16G16R16的格式。这就需要把fp16的波形方程结果转换为fp32格式。你可能会问,问什么开始不直接把波形方程的结果保存为fp32格式呢?答案是fp2的纹理渲染起来很慢,此外,fp32的渲染目标不支持混合或者纹理过滤。由于fp32缺乏混合机制,将导致水面控制更加复杂。而缺乏过滤,则会降低渲染质量。
![](http://blog.csdn.net/images/blog_csdn_net/soilwork/T2G.jpg)
好了,有了模拟水面的方法,接下来看看如何控制它。幸运的是,使用Verlet积分对此进行控制。Verlet积分本身就表示了高度偏移的速率。在GeForce 6系列的硬件上,高度偏移值将直接渲染为fp16格式的高度纹理。这是最自然也是最高效的方法。但是在Geforce 6系列之前的硬件上,这却是不可行的,应为缺乏对fp16 alpha混合的支持。没有混合,就不得不把高度偏移值渲染为额外的非fp16渲染目标中,之后,再把这张偏移值纹理渲染到高度图之上。这是很麻烦的方法,同时精度也不高。
再进一步来看看这个方法的性能如何,水面模拟到底有多快。虽然看起来有很多操作。幸运的是,这对GPU来说并算不了什么。你基本上只是在一个简单的pixel shader中渲染了一张小渲染目标(比如128x128)。凭感觉来说,这对现代GPU算不了什么。作为一个稍微“科学”的基准来看,在GeForce 6800GT的显卡上,关闭了水面模拟之后,帧速率从263fps变为了268fps。这表示我们的模拟只需要0.07毫秒,这并不算长。与实际的模型相比,顶点纹理拾取的性能才是瓶颈。
~~~~~~~~~~~~~~~~~~~·未完待续·~~~~~~~~~~~~~~~~~~~~~~~~
这篇文章是SDK里Vertex Texture Fetch Water的白皮书,今天头晕晕的,好多地方翻的不恰当,要看的同志还是对照原文一起来-_-#。 又快天亮了,偶睡觉去鸟~~~~
分类: GPU程序设计2006-05-19 06:34 5367人阅读
看了那么多文章,终于动手自己写水面渲染的程序了。选择了使用顶点纹理来实现(FFT。。。看了2遍信号与系统,还是没有明白怎么变换求解-_-#)。弄了3天,总算把水面的初步动画实现了,画面效果还没进行优化,纹理也没贴,不过效果还不错,CPU占用率出乎意料的低,网格大小为256x256,并且没有任何裁减的情况下,6600GT可以跑到70帧左右。本来打算用MDX2.0的,结果写了一点点,就出了一大堆问题,没有文档实在是麻烦的事情,还是回到了1.0。对顶点纹理进行采样的sampler参数设置差点让我崩溃了,查了无数资料,都没有关于如何在MDX里进行设置的例子。最后,模仿nVIDIA SDK里一个C++程序的例子写了一个,居然成功了~~。VS2003里不能调试shader实在是痛苦啊,但是拼写错误就害死人啊T_T。
记下一些容易犯的低级错误 ,以免再错-_-:
1. 顶点索引的顺序是很重要的, 切记需要按顺时针顺序来定义;
2. 矩阵相乘的顺序是有讲究的;
3. 拜托,写demo的时候先画个草图,设计好物体,观察点的位置,不要弄了半天发现物体根本不在视线内~_~;
4. 顶点纹理一定要是32bit的float point格式,否则就算程序可以运行,也不会有结果 +_+;
5. 注意tex2Dlod这个函数接收的纹理坐标参数为float4,而不是普通的float2 =.=;
6. 顶点纹理sampler的参数和pixar shader是不一样的,需要用特殊的标识符 @_@;
7. 拼写~~!!一定要注意拼写的正确性啊~~ 0_0;
![](https://p-blog.csdn.net/images/p_blog_csdn_net/soilwork/water.jpg)