Fortran OMP 随记 (1)

今天测试了Fortran OMP求$\Sigma_1^n i $ 串行和并行的速度。代码如下：

program mainuse OMP_LIBimplicit noninteger N,M,i, time_begin, time_endreal(kind=8) t, acall OMP_set_num_threads(4)N=2000000t=0.0call SYSTEM_CLOCK(time_begin)do i=1,Na=float(i)t=t+float(i)M = OMP_get_num_threads()enddocall SYSTEM_CLOCK(time_end)write(*, "('t =  ', F20.5, ' running on ', I3, ' threads,','using time =', I7)") t,M,time_end-time_begint=0.0call SYSTEM_CLOCK(time_begin)!$OMP PARALLEL DOdo i = 1, N!$OMP CRITICALa = float(i)t = t+a!$OMP END CRITICALM = OMP_get_num_threads()enddo!$OMP end parallel docall SYSTEM_CLOCK(time_end)write(*, "('t =  ', F20.5, ' running on ', I3, ' threads,','using time =', I7)") t,M,time_end-time_begint=0.0call SYSTEM_CLOCK(time_begin)!$OMP PARALLEL DO REDUCTION(+:t)do i = 1, Nt = t+float(i)M = OMP_get_num_threads()enddo!$OMP end parallel docall SYSTEM_CLOCK(time_end)write(*, "('t =  ', F20.5, ' running on ', I3, ' threads,','using time =', I7)") t,M,time_end-time_begin！The following codes are not correct.! It is right only under the serialt=0.0call SYSTEM_CLOCK(time_begin)!$OMP PARALLEL DOdo i=1,Na=float(i)t=t+aM=OMP_get_num_threads()enddo!$OMP end parallel docall SYSTEM_CLOCK(time_end)write(*, "('t =  ', F20.5, ' running on ', I3, ' threads,','using time =', I7)") t,M,time_end-time_beginend

下面是我的运行结果：

t =   2000001000000.00000 running on   1 threads,using time =     31
t =   2000001000000.00000 running on   4 threads,using time =   4383
t =   2000001000000.00000 running on   4 threads,using time =     47
t =    540364135587.00000 running on   4 threads,using time =     62

可能是因为电脑配置的原因，我双核四线程测试的并行结果并不好。有兴趣的小伙伴可以一起做个测试，把你们的结果放上来，看看怎么样

Fortran OMP 随记 (1)相关推荐

Perl特殊变量（flw等翻译）
NAME perlvar - Perl 预定义变量 DESCRIPTION 预定义名称后面列出的名称对 Perl 来说具有特殊含义. 大多数标点名称都有合理的助记方法或类 ...
粒子群算法Fortran代码(OMP并行)
粒子群算法可用于解决强非线性优化问题,原理较为简单(参加:最优化算法之粒子群算法(PSO)_青萍之末的博客-CSDN博客_粒子群算法),这里给出Fortran代码实现模块( module POS).该 ...
matlab中的expotest,提高Fortran矩阵指数性能（Expokit比Matlab、Python慢）
我正在进行一个模拟,其中的瓶颈是执行大量复杂的双精度矩阵指数运算,我发现Fortran(Expokit)对于小矩阵很好,但对于较大的矩阵,它的性能比Matlab或Python差.在我在下面包含了一个 ...
陆志勤：记交大数学系二三事
作者简介:陆志勤,上海交大数学系本科.硕士,柯朗研究所博士,现任教于加州大学尔湾分校数学系,研究方向为微分几何.陆志勤教授不仅研究成果丰硕,而且有许多高水平的数学普及文章和报告,见他的主页 https ...
fortran 语法学习之二
1.参考offloading_success.c的写法,将C++翻译为fortran的写法如下参考IBM Docs,第一个参数的默认属性是tofrom !$OMP TARGET MAP (from: ...
cuda FORTRAN 统一内存 managed
这里是官方的一篇很好的帖子,里面有实例和讲解 managed memory 的使用方法. http://www.pgroup.com/lit/articles/insider/v6n1a2.htm 可 ...
Fortran学习3：控制流2：循环
一个程序中循环是非常重要的,他会控制变量进行多次相同运算,在数组操作中有着不可替代的地位,今天就走进Fortran的循环,感受一下循环的乐趣. 哦,对了,Fortran程序在stop之前,可以加一个p ...
Fortran简明自学手册
Fortran简明自学手册因为CSDN不支持对Fortran代码的高亮显示,这儿贴出学习手册PDF版百度网盘的下载链接: https://pan.baidu.com/s/1A_6N0TQERg2n ...
Fortran语法汇总（下）（持续更新中）
目录第5章 I/O操作 5.1 格式描述符 5.2 OPEN语句 5.3 CLOSE语句 5.4 INQUIRE语句 5.5 READ语句 5.6 WRITE语句 5.7 文件定位语句 5.8 I/ ...

Fortran OMP 随记 (1)

可能是因为电脑配置的原因，我双核四线程测试的并行结果并不好。有兴趣的小伙伴可以一起做个测试，把你们的结果放上来，看看怎么样

Fortran OMP 随记 (1)相关推荐

最新文章

热门文章

Fortran OMP 随记 (1)

可能是因为电脑配置的原因，我双核四线程测试的并行结果并不好。 有兴趣的小伙伴可以一起做个测试，把你们的结果放上来，看看怎么样

Fortran OMP 随记 (1)相关推荐

最新文章

热门文章

可能是因为电脑配置的原因，我双核四线程测试的并行结果并不好。有兴趣的小伙伴可以一起做个测试，把你们的结果放上来，看看怎么样