再问一个代码优化的问题

shiter
人工智能领域优质创作者
博客专家认证
2016-01-22 03:42:44 




int qx_tree_filter::filter(double*cost,double*cost_backup,int nr_plane)

{

	memcpy(cost_backup,cost,sizeof(double)*m_h*m_w*nr_plane);

	int*node_id=m_node_id;

	int*node_idt=&(node_id[m_nr_pixel-1]);

	for(int i=0;i<m_nr_pixel;i++)//for each pixel in image

	{

		int id=*node_idt--;

		int id_=id*nr_plane;

		int nr_child=m_mst_nr_child[id];

		if(nr_child>0)

		{

			double*value_sum=&(cost_backup[id_]);

			for(int j=0;j<nr_child;j++)

			{

				int id_child=m_mst_children[id][j];

				int id_child_=id_child*nr_plane;

				double weight=m_table[m_mst_weight[id_child]];

				//value_sum+=m_mst_value_sum_aggregated_from_child_to_parent[id_child]*weight;

				double*value_child=&(cost_backup[id_child_]);

				for(int k=0;k<nr_plane;k++)

				{

					value_sum[k]+=(*value_child++)*weight;

				}

			}

		}

		//else

		//{

		//	memcpy(&(cost_backup[id_]),&(cost[id_]),sizeof(double)*nr_plane);

		//}

		//printf("[id-value-weight]: [%d - %3.3f - %3.3f]\n",id,m_mst_[id].value_sum_aggregated_from_child_to_parent,m_mst_[id].weight_sum_aggregated_from_child_to_parent);

	}

	int*node_id0=node_id;

	int tree_parent=*node_id0++;

	int tree_parent_=tree_parent*nr_plane;

	memcpy(&(cost[tree_parent_]),&(cost_backup[tree_parent_]),sizeof(double)*nr_plane);

	for(int i=1;i<m_nr_pixel;i++)//K_00=f(0,00)[K_0-f(0,00)J_00]+J_00, K_00: new value, J_00: old value, K_0: new value of K_00's parent

	{

		int id=*node_id0++;

		int id_=id*nr_plane;

		int parent=m_mst_parent[id];

		int parent_=parent*nr_plane;



		double*value_parent=&(cost[parent_]);//K_0

		double*value_current=&(cost_backup[id_]);//J_00

		double*value_final=&(cost[id_]);

		double weight=m_table[m_mst_weight[id]];//f(0,00)



		for(int k=0;k<nr_plane;k++) 

		{

			double vc=*value_current++;

			*value_final++=weight*((*value_parent++)-weight*vc)+vc;///这句怎么搞？

		}

		//printf("Final: [id-value-weight]: [%d - %3.3f - %3.3f]\n",id,m_mst_[id].value_sum_aggregated_from_parent_to_child,m_mst_[id].weight_sum_aggregated_from_parent_to_child);

	}

	return(0);

}











void qx_nonlocal_cost_aggregation::matching_cost_from_color_and_gradient(unsigned char ***left,unsigned char ***right)

{

	image_zero(m_cost_vol,m_h,m_w,m_nr_plane);

	compute_gradient(m_gradient_left,left);

	compute_gradient(m_gradient_right,right);



	for(int i=0;i<m_nr_plane;i++)

	{

		for(int y=0;y<m_h;y++) //shift the right image by i pixels

		{

			image_copy((&m_image_shifted[y][i]),right[y],m_w-i,3);

			memcpy(&(m_gradient_shifted[y][i]),m_gradient_right[y],sizeof(float)*(m_w-i));

			for(int x=0;x<i;x++) 

			{

				qx_memcpy_u3(m_image_shifted[y][x],right[y][0]);//m_cost_on_border_occlusion;

				m_gradient_shifted[y][x]=m_gradient_right[y][0];//m_cost_on_border_occlusion;

			}

		}

		//for(int y=0;y<m_h;y++) for(int x=0;x<m_w-1;x++) 

		for(int y=0;y<m_h;y++) for(int x=0;x<(m_w);x++) 

		{

			double cost=0;

			for(int c=0;c<3;c++) cost+=abs(left[y][x][c]-m_image_shifted[y][x][c]);

			cost=min(cost/3,m_max_color_difference);

			double cost_gradient=min((double)abs(m_gradient_left[y][x]-m_gradient_shifted[y][x]),m_max_gradient_difference);

			m_cost_vol[y][x][i]=m_weight_on_color*cost+m_weight_on_color_inv*cost_gradient;

		}

	}

}

...全文
302 23 打赏 收藏 转发到动态 举报
写回复
用AI写文章
23 条回复
切换为时间正序
请发表友善的回复…
发表回复
赵4老师 2016-01-27
  • 打赏
  • 举报
 回复
班门弄斧一下哈： 
register double*value_parent=&(cost[parent_]);//K_0
register double*value_current=&(cost_backup[id_]);//J_00
register double*value_final=&(cost[id_]);
register double weight=m_table[m_mst_weight[id]];//f(0,00)
……
register double vc=*value_current++;
            *value_final++=weight*((*value_parent++)-weight*vc)+vc;///这句怎么搞？
赵4老师 2016-01-27
  • 打赏
  • 举报
 回复
引用 21 楼 wangyaninglm 的回复:
[quote=引用 20 楼 zhao4zhong1 的回复:] [quote=引用 15 楼 wangyaninglm 的回复:] [quote=引用 12 楼 zhao4zhong1 的回复:] 请楼主每次提到“快”、“慢”等字眼时，都将我三楼回复的那句话默念至少三遍！
赵老师，我真的没次都是release而且有profile 啊。。。。[/quote] 那我问你你知道耗时最多的那几条语句在第多少行到第多少行吗？[/quote] 你看赵老师你都没认真看截图，我截图里面有呢。。。[/quote] 还真有！
shiter 2016-01-27
  • 打赏
  • 举报
 回复
引用 20 楼 zhao4zhong1 的回复:
[quote=引用 15 楼 wangyaninglm 的回复:] [quote=引用 12 楼 zhao4zhong1 的回复:] 请楼主每次提到“快”、“慢”等字眼时，都将我三楼回复的那句话默念至少三遍！
赵老师，我真的没次都是release而且有profile 啊。。。。[/quote] 那我问你你知道耗时最多的那几条语句在第多少行到第多少行吗？[/quote] 你看赵老师你都没认真看截图，我截图里面有呢。。。
赵4老师 2016-01-27
  • 打赏
  • 举报
 回复
引用 15 楼 wangyaninglm 的回复:
[quote=引用 12 楼 zhao4zhong1 的回复:] 请楼主每次提到“快”、“慢”等字眼时，都将我三楼回复的那句话默念至少三遍！
赵老师，我真的没次都是release而且有profile 啊。。。。[/quote] 那我问你你知道耗时最多的那几条语句在第多少行到第多少行吗？
日立奔腾浪潮微软松下联想 2016-01-26
  • 打赏
  • 举报
 回复
引用 11 楼 wangyaninglm 的回复:
[quote=引用 10 楼 DelphiGuy 的回复:] 可以考虑用intel C++开avx优化试一试，VC++也可以，但是优化水平比intel C++要落后一代，VC++ 2015的优化水平大致相当于intel C++ XE 2013的层次。
高手优化代码这方面的书给我推荐一个吧，就c++代码优化，你说的vs2015开avx是要重写代码还是编译器有相关选项呢？ 就是我还发现一个问题，vs2015的代码releas跑起来居然比vs2010 的release慢。。。同一台机器。。。[/quote] [quote=引用 11 楼 wangyaninglm 的回复:] 代码优化的书不少，但是大部分还是谈通用的优化原则、算法，象高效程序的奥秘这类书，针对具体架构优化的方法，还得看处理器厂家的资料，Intel 64 and IA-32 Architectures Optimization Reference Manual这种。 不用重写代码，就是/arch:AVX、/arch:AVX2这些编译选项，但是VC的优化水平还是比intel c++差不少，经常是你开了特定指令集优化，但是它并没有用到。某些代码编译后比2010编译的代码慢的问题在2013上就有，看反汇编却没有明显的差异。
日立奔腾浪潮微软松下联想 2016-01-26
  • 打赏
  • 举报
 回复
引用 11 楼 wangyaninglm 的回复:
[quote=引用 10 楼 DelphiGuy 的回复:] 可以考虑用intel C++开avx优化试一试，VC++也可以，但是优化水平比intel C++要落后一代，VC++ 2015的优化水平大致相当于intel C++ XE 2013的层次。
高手优化代码这方面的书给我推荐一个吧，就c++代码优化，你说的vs2015开avx是要重写代码还是编译器有相关选项呢？ 就是我还发现一个问题，vs2015的代码releas跑起来居然比vs2010 的release慢。。。同一台机器。。。[/quote] 前面的回复怎么乱套了。 代码优化的书不少，但是大部分还是谈通用的优化原则、算法，象高效程序的奥秘这类书，针对具体架构优化的方法，还得看处理器厂家的资料，Intel 64 and IA-32 Architectures Optimization Reference Manual这种。 不用重写代码，就是/arch:AVX、/arch:AVX2这些编译选项，但是VC的优化水平还是比intel c++差不少，经常是你开了特定指令集优化，但是它并没有用到。某些代码编译后比2010编译的代码慢的问题在2013上就有，看反汇编却没有明显的差异。
赵4老师 2016-01-26
  • 打赏
  • 举报
 回复
请楼主每次提到“快”、“慢”等字眼时，都将我三楼回复的那句话默念至少三遍！
yshuise 2016-01-26
  • 打赏
  • 举报
 回复
我以前用的是微软编译器vs2008，其他的没有用
yshuise 2016-01-26
  • 打赏
  • 举报
 回复
胆子要大，对for循环inline是必杀技，不信你可以测试下
shiter 2016-01-26
  • 打赏
  • 举报
 回复
引用 16 楼 yshuise 的回复:
inline即可，早就试过了，很有效。建议楼主试试。
这么长俩函数定义成inline？
yshuise 2016-01-26
  • 打赏
  • 举报
 回复
inline即可，早就试过了，很有效。建议楼主试试。
shiter 2016-01-26
  • 打赏
  • 举报
 回复
引用 12 楼 zhao4zhong1 的回复:
请楼主每次提到“快”、“慢”等字眼时，都将我三楼回复的那句话默念至少三遍！
赵老师，我真的没次都是release而且有profile 啊。。。。
shiter 2016-01-25
  • 打赏
  • 举报
 回复
引用 10 楼 DelphiGuy 的回复:
可以考虑用intel C++开avx优化试一试，VC++也可以，但是优化水平比intel C++要落后一代，VC++ 2015的优化水平大致相当于intel C++ XE 2013的层次。
高手优化代码这方面的书给我推荐一个吧，就c++代码优化，你说的vs2015开avx是要重写代码还是编译器有相关选项呢？ 就是我还发现一个问题，vs2015的代码releas跑起来居然比vs2010 的release慢。。。同一台机器。。。
日立奔腾浪潮微软松下联想 2016-01-25
  • 打赏
  • 举报
 回复
可以考虑用intel C++开avx优化试一试，VC++也可以，但是优化水平比intel C++要落后一代，VC++ 2015的优化水平大致相当于intel C++ XE 2013的层次。
paschen 版主 2016-01-24
  • 打赏
  • 举报
 回复
楼主也是一个效率控
shiter 2016-01-24
  • 打赏
  • 举报
 回复
引用 6 楼 DelphiGuy 的回复:
数据结构未知啊
变量都是double，以及一些double，float数组
shiter 2016-01-24
  • 打赏
  • 举报
 回复
引用 8 楼 paschen 的回复:
楼主也是一个效率控
还请大牛指导，我也是刚接触优化这个方向
赵4老师 2016-01-23
  • 打赏
  • 举报
 回复
引用 4 楼 wangyaninglm 的回复:
[quote=引用 2 楼 zhao4zhong1 的回复:] 至少得生成Release版。 在MSDN里面搜“MMX, SSE, and SSE2 Intrinsics”
引用 2 楼 zhao4zhong1 的回复:
至少得生成Release版。 在MSDN里面搜“MMX, SSE, and SSE2 Intrinsics”
引用 3 楼 zhao4zhong1 的回复:
无profiler不要谈效率！！尤其在这个云计算、虚拟机、模拟器、CUDA、多核 、多级cache、指令流水线、多种存储介质、……满天飞的时代！
赵老师，我这个p'rofile而且是release的哇[/quote] 那就好，我只是提醒你注意。
日立奔腾浪潮微软松下联想 2016-01-23
  • 打赏
  • 举报
 回复
数据结构未知啊
shiter 2016-01-22
  • 打赏
  • 举报
 回复
引用 2 楼 zhao4zhong1 的回复:
至少得生成Release版。 在MSDN里面搜“MMX, SSE, and SSE2 Intrinsics”
引用 2 楼 zhao4zhong1 的回复:
至少得生成Release版。 在MSDN里面搜“MMX, SSE, and SSE2 Intrinsics”
引用 3 楼 zhao4zhong1 的回复:
无profiler不要谈效率！！尤其在这个云计算、虚拟机、模拟器、CUDA、多核 、多级cache、指令流水线、多种存储介质、……满天飞的时代！
赵老师，我这个p'rofile而且是release的哇
加载更多回复（3）
MySQL性能优化全过程落地方案（工业界真实数据场景实操）
在互联网公司中，MySQL用的非常广泛，通常系统产品实现初期 一般要求实现功能，然后再后期数据量不断增大的情况下，系统的性能就出现了问题，其中MySQL性能影响占据很大一部分，所以我们需要了解MySQL原理，以及优化的专业知识，才能对mysql性能优化得心应手。学习该课程的价值：1增加自己的知识深度，学习完MySQL原理，其他数据库产品类似，达到触类旁通。2、一般企业或者大厂都要求会MySQL性能优化，面试过程也会问到，学会后不仅提高自己的专业水平，而且可以给自己跳槽涨薪机会。本课程基于电商环境测试数据讲解，数据量400w+。课程结尾会基于真实的环境进行MySQL调优讲解，让大家在实战中掌握MySQL性能知识，能够把调优技术真实用到实际生产环境中。本课程包含的技术： IDEA集成开发工具 MySQL 架构存储引擎SQL优化基础知识MySQL日志执行计划MySQL索引Spring Boot Spring CloudMysql性能优化实战等
【Lua进阶系列】之代码优化，别再依赖if..else了
【Lua基础系列】之ifElse代码优化 大家好，我是Lampard~~ 欢迎来到Lua基础系列的博客 前文再续，书接上一回。 今天讲解一下我在工作中遇到的代码优化问题。需求是这样的：我需要在项目的通用资源弹窗中加入一个新的资源道具，但是以前的通用弹窗道具可以换算成元宝，有自己的容量。但是新的道具是“英雄经验”是没有上述的属性的，所以有部分的逻辑是不适用于新增的道具，因此我一开始的处理方式是直接在各个判断的逻...
【代码优化】有了这45个小技巧，再也不怕代码写得烂了！!
比如在项目中不同的类型的业务可能需要上传各种各样的附件，此时就可以定义好不同的一个附件的枚举，来区分不同业务的附件。不要在代码中直接写死，不定义枚举，代码阅读起来非常困难，直接看到数字都是懵逼的。。
Web前端开发学习3：SEO代码优化
一代码优化概述        关于代码优化的知识是纯理论的知识，学习的很枯燥。在学到CSS时，不免遇到CSS+div进行代码优化的知 识，因此在网上看了一些关于这方面的知识，简单的整合一下，梳理自己所了解的代码优化问题。         所谓代码优化是指对程序代码进行等价（指不改变程序的运行结果）变换。程序代码可以是中间代码，也可以是 目标代码。等价的含义是使得变换后的代码运行结果与变换前
KEIL里的优化选项解读
优化大家都知道，-O0,-O1,-O2,-O3,等级越高优化程度越高，代码的执行效率也就越高。3、Optimize for time的3种状态是什么意思？1、低优化等级代码能跑，开启高优化等级后代码就跑不起来了。2、高优化等级代码能跑，开启低优化等级后代码就跑不起来了。第3个问题还是在群里问大神才知道。问题1和2的解决思路看看。
C++ 语言

64,662

社区成员

250,488

社区内容

发帖
与我相关
我的任务
社区描述
C++ 语言相关问题讨论，技术干货分享，前沿动态等
c++ 技术论坛（原bbs）
社区管理员
  • C++ 语言社区
  • encoderlee
  • paschen
加入社区
  • 近7日
  • 近30日
  • 至今

加载中

查看更多榜单
社区公告
  1. 请不要发布与C++技术无关的贴子
  2. 请不要发布与技术无关的招聘、广告的帖子
  3. 请尽可能的描述清楚你的问题，如果涉及到代码请尽可能的格式化一下

试试用AI创作助手写篇文章吧

+ 用AI写文章