C++随机数导致的bug

这周我在写path tracing的时候遇到了一个很奇怪的bug:光线追踪深度设为1的时候结果看起来是正常的,但是深度高于1的时候就会出现很奇怪的渲染结果,如下图所示(为了便于调试,我用了一个cornell-box的场景,mesh的颜色是统一的,并且设置的只是点光源):

我仔细排查了逻辑上的代码,并没有发现几何计算、光照计算等模块有明显的问题(如果有问题,应该在depth为1的时候就测试出来了)。这时我观察到,这个颜色呈现出“一块一块”的表现,我猜测可能是由于在光线做完第一次碰撞后,利用随机数产生出射光线进行第二次碰撞,但是随机数分布并不均匀导致的。
于是我试着把代码中所有的srand全部去掉,果然结果就正常了。也就是说,是srand导致了产生的随机数分布不均匀
那么我的srand是怎么用的呢?

Pixel PathTracingRenderer::sample_pixel(Camera* camera, int x, int y, int width, int height)
{
    srand(time(0));
    float sample_x = 0.0f;
    float sample_y = 0.0f;
    Spectrum r;
    Pixel p;
    for (int i = 0; i < sample_count; i++)
    {
            //compute sample result
    }
    return p;
}

我将srand写到了遍历每个像素的代码里,这明显是不对的。我们可以做个试验来验证一下:
下图是我不使用srand,完全随机生成的一堆二维点,可以看到基本上是随机分布的:

假设我这样写:

for (int i = 0; i < 100; i++)
{
    srand(time(0));
    std::cout << rand() % 101 << "," << rand() % 101 << ";"<<std::endl;
}

生成的结果会是如下所示的(请原谅我看到结果时笑出了声):

之后我又做了很多次试验,发现只有每次生成随机数间隔时间差不多1s,这样写才能生成比较随机的结果:

所以导致这个bug的答案也就呼之欲出了:srand会重置随机生成的结果,输入同一个种子,那么第一次rand的结果是一样的。之所以会导致如图所示的问题,就是因为相邻像素之间重复调用srand,而相隔的时间又很短,不断重置随机数结果,始终生成了相同的随机数,那么自然就会出现上图这种“颜色聚集”的情况了。
于是我就把代码改成了每次只有在渲染前才会srand一次,这样一来结果就正常了(和PBRT相同的场景与光照比对基本一致):

原文地址:https://www.cnblogs.com/wickedpriest/p/12443131.html

时间: 2024-08-30 01:43:04

C++随机数导致的bug的相关文章

对非正确使用浮点型数据而导致项目BUG的问题探讨

乘法分配律 在上小学的时候就已经学习过乘法分配律,乘法分配律的具体内容是:两个数的和与一个数相乘,可以先把他们分别与这个数相乘,再相加,得数不变.乘法分配律的定义还可以用表达式"(a+b)×c = a×c+b×c"的形式给出.乘法分配律的反用"a×c+b×c = (a+b)×c"同样成立.例如"10.2×(3+7) = 10.2×3+10.2×7 = 102"(反用形式为"10.2×3+10.2×7 = 10.2×(3+7) = 102

结构体类型重声明导致的bug一个

bug前提条件 当模块比较多,头文件较多,某个结构体类型会在当前模块中重新声明进而引用其成员,而不直接包含其他模块的头文件.这样的好处是不引入不需要的类型声明到此模块,头文件包含的交叉:坏处是,增加了bug的几率,耦合太大!比如下面一种情况发生而导致bug: 已知两个模块A和B,同一个结构类型struct node在两个模块中分别声明,其中B模块无意或者有意调整了结构类型中的某些域.那么这个时候,若B模块中引用A模块中此类型实例然后访问成员变量,就会引发bug!如下重现示例 bug重新示例代码

在docker以FPM-PHP运行php,慢日志导致的BUG分析

问题描述: 最近将IOS书城容器化,切换流量后.正常的业务测试了一般,都没发现问题.线上的错误监控系统也没有报警,以为迁移工作又告一段落了,暗暗的松了一口气.紧接着,报警邮件来了,查看发现是一个苹果支付相关接口调用的curl错误,错误码为"56",错误描述为:“Failure with receiving network data”接收网络数据失败. 机器 : 192.168.1.1当前URL : /xxx/recharge/apple?xxxxxxxxxxxxx接口URL : htt

高并发压力下导致数据库bug

环境信息:  linux 6.1 + oracle11.2.0.3 RAC 问题现象: 学校晚上6点选课,人数大概有3000,7点时,数据库报错如下(数据库到6点多还是可以连接的),数据库hung住了. Tue Dec 16 18:00:33 2014Dumping diagnostic data in directory=[cdmp_20141216180033], requested by (instance=2, osid=24917 (M001)), summary=[incident=

关于硬盘满了导致的bug

在测试的过程中,我们往往会遇到很多稀奇古怪的bug,而研发有时一时半会也找不到原因,往往是基于我们之前的测试经验,来迅速排查问题. 1.在刚开始进入测试行业不久的时候,遇到一个问题,在服务器端,抓包发现每次抓的包都是无数据,刚开始排查的原因是语句出错,其实不然,发现在home下,挂的盘满了,导致无法存储数据 2.在参加项目是,也遇到一个同样的问题.搜人的项目,查看日志,都发现没有任何的异常,但是处理的图片没有保存,并且在数据库也没有保存图片的路径,依据经验,硬盘空间满了,df -h 看了一下,证

PHP 数组 foreach引用导致的bug

<?php $arr = array('a', 'b', 'c'); foreach($arr as &$v) {} foreach($arr as $v) {} var_dump($arr); /* 结果为: array(3) { [0]=> string(1) "a" [1]=> string(1) "b" [2]=> &string(1) "b" } */ //第一个foreach用引用赋值的方式将数

一次疏忽导致的bug

NBB_PUT_SHORT 这个宏是按char* 类型算指针的实际工作中,没有注意这一点,输入指针类型 强转为了 SHORT* 导致 填充的内容错误 这是第一次在工作中遇到 指针类型相关的问题值得记录一笔 1 // ConsoleApplication1.cpp : 定义控制台应用程序的入口点. 2 // 3 4 #include "stdafx.h" 5 #include <stdio.h> 6 #include <string.h> 7 #include &

input name属性太长导致的BUG

<td ><input id="sttotalResourceEqualsCost" name="sttotalResourceEqualsCost" type="text" class="easyui-numberbox" data-options="min:0,precision:2" style="width:80px" /> 这个一个easyui的numb

Unity3D BakeMesh 因二进制文件读写字符串的方式导致莫名bug

由其他人提供的插件,bake mesh static void writeString(BinaryWriter bw, string str, int cap) { var ba = new byte[cap]; var bs = System.Text.Encoding.ASCII.GetBytes(str); Array.Copy(bs, ba, bs.Length); bw.Write(ba); } 直接导致了当字符串过长的时候,内容读取的缺失和后续连锁导致序列化的错误. 解决方案: 1