文号
t18436
阅读
26163
回复
12
收藏
5
大爱的CUDA：图像亮度直方图统计
小俊
2009/08/07
软件综合
IP:广东
（发些代码赚些零用KCB来花花）

CPU实现大法：计算每个点的亮度值（0~255），然后用一个256大小的数组统计每个亮度的点的数量。



对1280*1024大小的图像，在QX6600（2.4GHz）上单线程运行时间为17.5ms。

hist_cpu.jpg
+777
科创币
破93
2009/08/07
70马要赞扬
来自：
计算机科学/软件综合
12
全部只看作者
小俊作者
17年5个月前IP:未同步
137898
1楼
CUDA实现大法：和CPU法原理差不多，但实现起来要分两步。第一步把整个图像分成很多个部分，每个部分用一个线程块（含很多个线程）独立统计；第二步则把各个部分独立统计的结果再累加起来得到总的统计结果。之所以用这个方法是从GPU内部架构考虑的，因为每个部分并行独立统计的结果可以放在对应的片上shared memory中，效率很高。代码比较复杂：



同样大小的图像，在Geforce GTX 285上计算时间为0.56ms，比QX6600上单线程提速约30倍。

hist_gpu.jpg
引用
评论
novakon
17年5个月前IP:未同步
137938
2楼
小俊那里是否能便宜拿到显卡?geforce8XXX系列我想弄一块,越便宜越好,最好100元以内解决
引用
评论
phpskycn
17年5个月前IP:未同步
138010
3楼
我能300搞到块9600GT-512MB
等我买新机器装上
引用
评论
小俊作者
17年5个月前IP:未同步
138140
4楼
引用第3楼novakon于2009-08-08 10:05发表的  :
小俊那里是否能便宜拿到显卡?geforce8XXX系列我想弄一块,越便宜越好,最好100元以内解决


弄不到。我现在台式机用的显卡和主板是在NV美国总部的垃圾堆捡回来的。
引用
评论
darkstorm
17年5个月前IP:未同步
138264
5楼
对了，经研究发现最后一步可以用对global atomic实现
第二个kernel可以省掉
以前理解有误
引用
评论
小俊作者
17年5个月前IP:未同步
138281
6楼
引用第6楼darkstorm于2009-08-08 23:06发表的  :
对了，经研究发现最后一步可以用对global atomic实现
第二个kernel可以省掉
以前理解有误


是不是在__syncthreads之后，直接把各线程块的shared memory的内容通过AtomicAdd累加到global memory中呢？
引用
评论
小俊作者
17年5个月前IP:未同步
138284
7楼
这样？

atomicAdd((int *)&out[tid], smem[tid]);
atomicAdd((int *)&out[tid + 128], smem[tid + 128]);
引用
评论
acmilan
17年5个月前IP:未同步
138337
8楼
好像在科创搞计算机的都一夜暴富了呵。。。
我还是一穷光蛋。。。
引用
评论
darkstorm
17年5个月前IP:未同步
138391
9楼
引用第8楼小俊于2009-08-08 23:40发表的  :
这样？

atomicAdd((int *)&out[tid], smem[tid]);
atomicAdd((int *)&out[tid + 128], smem[tid + 128]);


然
引用
评论
小俊作者
17年5个月前IP:未同步
138828
10楼
引用第10楼darkstorm于2009-08-09 11:39发表的  :


然


今天试了一下，改成对global memory原子加之后，计算时间从0.56ms减少到0.55ms，还是有一定作用的，多谢指教！
引用
评论
novakon
17年5个月前IP:未同步
141557
11楼
引用第5楼小俊于09-08-08 18:30发表的  :


弄不到。我现在台式机用的显卡和主板是在NV美国总部的垃圾堆捡回来的。

啊(突然头很痛)。
引用
评论
ki帛
5个月29天前IP:重庆
949914
12楼

话说咱KC的大佬们都那么超前吗1f605那时估计CUDA才刚出来没多久，知名度不是很高吧

引用
评论
全部只看作者
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的