- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在使用英特尔 Fortran 编译器和英特尔 mkl 进行性能检查。我通过像这样的调用将一些数组部分传递给 Fortran 77 接口(interface)
call dgemm( transa,transb,sz_s,P,P,&
a, Ts_tilde,&
sz_s,R_alpha,P,b,tr(:sz_s,:),sz_s)
很明显,tr(:sz_s,:) 在内存中不是连续的,Fortran 77 接口(interface)需要一个连续的 block 并为此创建一个临时 block 。
我想知道的是,如果我在 tr 的代码中显式创建临时数组并在操作之前和之后从该临时数组中来回复制信息,是否会有所不同,或者与编译器本身相同从性能的角度来看创建临时的?我想编译器总是会更高效。
当然,欢迎提出更多消除这些临时文件的建议。
还有一点,如果我显然使用该库的 Fortran 95 接口(interface),并对更简单的测试问题进行类似的调用,则不会针对创建临时对象发出警告。然后我在 mkl 手册中读到 Fortran 95 接口(interface)使用假定的形状数组,这解释了为什么不创建临时数组。
但是在这一点上,我似乎无法使用一些支持功能,例如定时例程。也就是说,intel mkl 有一些计时支持函数,但如果我将它们与下面的 mkl_service 例程一起使用,那么我会收到 dsecnd 的“此名称没有类型,并且必须具有显式类型”错误。也欢迎对此问题的任何想法。一个简单的例子如下
program dgemm95_test
! some modules for Fortran 95 interface
use mkl_service
use mkl95_precision
use mkl95_blas
!
implicit none
!
double precision, dimension(4,3) :: a
double precision, dimension(6,4) :: b
double precision, dimension(5,5) :: r ! result array
double precision, dimension(3,2) :: dummy_b
!
character(len=1) :: transa
character(len=1) :: transb
!
double precision :: alpha, beta, t1, t2, t
integer :: sz1, sz2
! initialize some variables
alpha = 1.0
beta = 0.0
a = 2.3
b = 4.5
r = 0.0
transa = 'n'
transb = 'n'
dummy_b = 0.0
! Fortran 95 interface
t1 = dsecnd()
call gemm( a, b(4:6,1:3:2), r(2:5,3:4),&
transa, transb, alpha, beta )
t2 = dsecnd()
!
write(*,*) r
dummy_b = r(2:4,4:5)
!
end program dgemm95_test
最佳答案
将数组部分传递给旧例程使用的假定大小的数组虚拟参数时,临时参数是绝对必要的,因为数组部分在内存中不是连续的。
您当然可以创建自己的临时数组。是否会更快取决于很多因素。其中重要的是临时数据是分配在堆栈上还是堆上。英特尔 Fortran 编译器能够同时执行这两种操作,有编译器开关来控制行为 (-heap-arrays n
),并且它可以取决于数组大小。堆栈分配速度要快得多,并且通常是默认分配。默认情况下,您可以将其用于自己的临时数组的自动数组也分配在堆栈上。小心堆栈上的大数组,您很容易溢出并导致崩溃。
我建议您进行性能测试并使用更简单的变体(如果它不是太慢的话)。可能是 Fortran 95 界面,但您确实应该测量时间。
至于时间,second()
/dsecnd()
的 MKL 手册页声明您必须包含mkl_lapack.fi
,但没有包含“更不用说任何 Fortran95 接口(interface)了。您也可以将其声明为外部 double ,但我会使用 include。或者使用 system_clock()
作为可移植标准 Fortran 95。
关于fortran - 在 Fortran 中传递不连续的数组部分,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/24474371/
如何检查一个元素是否立即隐藏。即如何通知元素的可见性。 在我的例子中,该元素是通过 slideUp 函数隐藏的。我应该立即收到有关该元素的可见性的通知。 我想到了使用bind()方法。但它没有类似 o
if (srcbloc == NULL) { fprintf(stderr, "warning!: memrip source is null!\n"); exit(1); } if
当我在数据库的旧 View 中清理一些问题时,我遇到了这个“奇怪”的连接条件: from tblEmails [e] join tblPersonEmails [pe]
如何水平对齐多张图像,一张一张地?它们不必适合宽度屏幕:相反,我希望它们超过后者的宽度,如果这有任何意义的话。 我已经检查了很多类似问题的答案,但找不到任何可以解决我的问题的答案。 HTML:
我知道 Cassandra 中的列有 TTL。但是也可以在一行上设置 TTL 吗?在每列上设置 TTL 并不能解决我的问题,如下面的用例所示: 在某些时候,一个进程想要删除一个带有 TTL 的完整行(
我有一个 NSTextField 和 Label,其值绑定(bind)到 View Controller 中的相同 NSString 这里的问题是标签只有在我按 Tab 时才会更新。 如何使其连续,以
例如。 1."abc"; ===>abc 2."ab c"; ===>ab_c 3."ab c"; ===>ab_c 4."ab c" ===>ab_c 对于多个连续空格也是如此。 我怎样
大家好,我想获取前一天或最后一天的信息,只有当我按下按钮时,它才会显示最后一天(星期六)的所有信息,如果我再次单击按钮,它将显示最后一天的信息(星期五)如果我再次点击它(星期四)谢谢你们帮助我 编辑:
我需要从实时音频流中提取ICY元数据,并正在使用mplayer进行此操作,因为它在播放音频流时会输出元数据。我欢迎其他方式执行此操作,目标是将更新的元数据(歌曲信息)保存到文本文件中,只要歌曲(或数据
语音识别有没有解决方案 只有几个字(2 个就够了,10 个就不错了。100 个就很棒了。不需要更多) 也在移动浏览器上运行(是否可以为此使用 flash(而不是 java)?) 可以安装在您自己的服务
我有一个单词列表, list1 = ['hello', 'how', 'are', 'you?', 'i', 'am', 'fine', 'thanks.', 'great!'] 我想加入, list
我正在开发一个程序,但我不断收到“对‘dosell’的 undefined reference ”,我不太明白发生了什么。这是函数的声明: void dosell(int *cash, int *nu
我无法提出执行我要做的事情所需的查询。 我有三个这样的表: client_files ----------------------- client_id file_id ---------
我一直在寻找一个插件/脚本,当到达底部时,它会从头开始继续滚动网站,就像一个连续的循环。 示例:http://unfold.no/和 http://www.aquiesdonde.com.ar/ 我尝
这个问题在这里已经有了答案: How to prevent scanf causing a buffer overflow in C? (6 个答案) 关闭 6 年前。 我一直在使用一个非常简单的程
给定一个整数数组,找到具有相同数量的 x 和 y 的连续子序列的总数。例如 x=1 和 y=2 的数组 [1,2,1] ans = 2 表示它的两个子数组 [1,2] 和 [2,1]。检查每个连续的子
所以,我有一个所有正自然数的数组。我得到了一个阈值。我必须找出总和小于给定阈值的数字(连续)的最大计数。 For example, IP: arr = {3,1,2,1} Threshold = 5
我制作了像内置相机一样的相机应用。 我想实现像内置相机一样的连续对焦功能。(此功能我不触摸屏幕,但相机会尝试自行对焦。) 因此,将其设置为 surfaceCreated : Camera.Pa
我有这样的数据: f x A 1.1 A 2.2 A 3.3 B 3.5 B 3.7 B 3.9 B 4.1 B 4.5 A 5.1 A 5.2 C 5.4 C 5.5 C 6.1 B 6.2 B
假设我有一个包含一组数据点的表,每个数据点由一个时间戳和一个值组成。如果至少有 N 个连续记录(按时间戳排序)高于给定值 X,我将如何编写返回 true (1) 的查询,否则返回 false (0)?
我是一名优秀的程序员,十分优秀!