matlab中的expotest,提高Fortran矩阵指数性能（Expokit比Matlab、Python慢）
我正在进行一个模拟，其中的瓶颈是执行大量复杂的双精度矩阵指数运算，我发现Fortran(Expokit)对于小矩阵很好，但对于较大的矩阵，它的性能比Matlab或Python差。在
我在下面包含了一个显示类似行为的模型程序，尽管它需要更大的矩阵来显示性能差异。看看profiler和source code似乎Expokit大部分时间都在调用zgemm()，所以我唯一的想法是BLAS安装有问题。否则我不明白为什么Fortran的性能会比Matlab或Python差。如果您能深入了解如何改进Fortran矩阵指数代码的性能，我将不胜感激。在Results for 10000 matrices (4x4, 8x8, 30x30, 60x60, 80x80):
Matlab (s): 0.91, 0.97, 2.36, 5.45, 8.69
Python (s): 2.59, 2.89, 9.70, 35.4, 72.7
Fortran, Expokit (s): 0.037, 0.12, 4.14, 30.6, 74.9
Fortran, Expokit, OpenMP with 8 cores (s): 0.0039, 0.016, 0.52, 3.87,
9.53
Fortran代码：subroutine expokit_test()
use omp_lib
use iso_fortran_env
implicit none
integer, parameter :: wp = selected_real_kind(15, 307), size=80
complex(wp), parameter :: i = (0, 1._wp)
integer :: count, a, b
real(wp) :: wtime
complex(wp) :: mat_exp(size, size), mat(size, size), val
val = 1E-8_wp
mat = 0._wp
do a = 1, size
do b = 1, size
mat(a, b) = a * b
end do
end do
call omp_set_num_threads(8)
wtime = omp_get_wtime()
!$omp parallel do default(private) &
!$omp& shared(mat, val)
do count = 1, int(1E4)
mat_exp = expm_complex(-i * mat * val)
end do
!$omp end parallel do
wtime = omp_get_wtime () - wtime
write(6, *) 'expm_complex', sngl(wtime)
end subroutine expokit_test
function expm_complex(A) result(B)
! Calculate matrix exponential of complex matrix A using Expokit
use iso_fortran_env
implicit none
integer, parameter :: wp = selected_real_kind(15, 307)
complex(wp), dimension(:, :), intent(in) :: A
complex(wp), dimension(size(A, 1), size(A, 2)) :: B
integer, parameter :: ideg = 2 ! Pade approximation, 6 is reccomended but 2 appears to be stable
complex(wp) :: t = 1._wp
complex(wp), dimension(4 * size(A, 1) * size(A, 2) + ideg + 1) :: wsp
integer, dimension(size(A, 1)) :: iwsp
integer :: iexp, ns, iflag, n
n = size(A, 1)
call ZGPADM(ideg, n, t, A, n, wsp, size(wsp, 1), iwsp, iexp, ns, iflag)
B = reshape(wsp(iexp : iexp + n * n - 1), [n, n])
end function expm_complex
Matlab代码：
^{pr2}$
Python代码：size = 80
mat = np.ones((size, size))
for a in range(0, size):
for b in range(0, size):
mat[a, b] = ((a+1) * (b+1))
mat = mat + 1E-9
start = time.time()
for loop in range(0, int(1E4)):
test = la.expm(-1j * mat * 1E-8)
end = time.time() - start
print('time taken', end)
相关文章：
Win7下使用Putty代替超级终端通过COM串口连接开发板方法
1、如果电脑（笔记本）没有串口接口，则需要使用一个 USB-Serial 转换线，这里使用 prolific usb-serial USB--串口转换线，首先需要在win7上安装对应的 USB--串口转换线 驱动程序，PL2303_Prolific 驱动程序下载…
java内部类 缺点_Java中的内部类及其优势
Java中提供了定义内部类的选择，这一机制使得代码的书写更为方便和优雅(功能上相关的代码被紧密的组织在了一起)。需要注意的是，内部类和传统的组合(即将一个类的实体定义为另一个类的成员)是完全不同的，其主要特性表现在以下的一些方面&#…
php webapi验签,Asp.netCore3.0 WebApi从0到1手摸手教你写【5】增加接口参数签名验证...
通过前几个教程的学习，对webapi的编写基本上就可以入门了，可以做项目了，今天我们再给接口加个参数签名认证，之前的接口相当于赤果果的暴露在了网络上，只要知道接口地址、接口调用方式和传参就可以畅所欲为的调用接口了…
点击TableView任一行跳转详情页面会跳转两次的解决办法
在做TableView跳转的时候，发现实际上生成了两个detail view。我 navigate back 的时候，也是先看到一次detail view，然后才回到tableView的。 这是因为：performSegue(withIdentifier: , sender: ) 和 prepare(for segue: , sender:…
dos环境下mysql的访问_MYSQL dos环境下使用
有很多朋友虽然安装好了mysql但却不知如何使用它。在这篇文章中我们就从连接MYSQL、修改密码、增加用户等方面来学习一些MYSQL的常用命令。一、连接MYSQL。格式： mysql -h主机地址 -u用户名 －p用户密码1、例1：连接到本机上的MYSQL。首先在打开…
jemeter监听器的使用
打开jemeter，新建线程组，添加http请求，在请求下添加监听器： 一、添加一个jpgc - PerfMon Metrics Collector监听器： 服务器性能监测控件，包括CPU,memory（内存）,networkＩ/…
java 抛异常给上级_java异常处理机制(示例代码)
Exception 类的层次java中所有的异常类是从 java.lang.Exception 类继承的子类。而Exception 类是 Throwable (可抛出的)类的子类。除了Exception类外，Throwable还有一个子类Error 。Java 程序通常不捕获错误。错误一般发生在严重故障时，它们在Java程序处…
前端开发之JavaScript基础篇一
主要内容： 1、JavaScript介绍 2、JavaScript的引入方法和输出及注释 3、javaScript变量和命名规则 4、五种基本数据类型 5、运算符 6、字符串处理 7、数据类型转换 一、JavaScript介绍 1、JavaScript是什么？ javaScript是一种web前端的描述语言&…
v9php 碎片信息,phpcms v9碎片管理及调用技巧分享
今天在这里将分享下Phpcms V9碎片管理及调用技巧。这是关于模板数据自定义、方便客户在后台管理数据调用的一个技巧。在给客户定制模板的时候，往往会涉及到的一个问题就是：有些数据(图片文字，比如LOGO、侧栏的联系方式、首页的幻灯片切换Bann…
iterm2 主题_【超级实用】Iterm2 + ohmyzsh 打造强大的终端编辑器
作者：AndrewHR 地址：http://mrw.so/4D1n7B最终的效果图如下所示：使用iterm2配合oh-my-zsh的命令行，拥有语法高亮、命令自动补全、自动提示符、显示git仓库状态等功能整个配置流程1、安装iterm2首先我们下载的 iterm2 这个软件&…
main方法为什么是静态的
main函数其实也是所在类的一个方法，就比如一个类是test，那么该类的main其实就是test.main(String[] args)，众所周知如果一个方法不是静态的，则要先实例化该类，比如要这样 test tnew test(); 然后才能调用 test.main();…
java类为什么要建两个class_ClassLoader的几个概念、类和对象的解释
首先，转载一篇文章，个人认为是看到过了讲得最清楚的 XD当JVM(Java虚拟机)启动时，会形成由三个类加载器组成的初始类加载器层次结构：bootstrap classloader|extension classloader|system classloaderbootstrap classloader &#…
第四章 python的turtle库的运用
我们可以尝试用python的自带turtle库绘制一条蟒蛇 首先我们设计一下蟒蛇的基本形状 我们先把这段蟒蛇绘制的实例代码贴出来，各位可以在自己的本地运行一下看看效果，然后我们再继续分析代码： 1 #PythonDraw.py2 import turtle3 turtle.setup(6…
oracle导入索引b报错,impdp导入索引很慢
impdp用NETWORK_LINK从远程库导入到本地库，导入表的速度还正常，导入索引的速度特别慢。2个小时才导1300个索引。使用imp的格式：impdp vebackup/abc DIRECTORYDATABAK_DIR NETWORK_LINKprimarydb_link SCHEMASveasms TABLE_EXISTS_ACTIONREPLA…
新的mysql如何使用_如何使用新的MySQL更新日志
使用新的MySQL更新日志的方法未必人人都会，下面就教您如何使用新的MySQL更新日志的方法，希望对您能够有所帮助。如果你只使用一使用新的MySQL更新日志的方法未必人人都会，下面就教您如何使用新的MySQL更新日志的方法，希望对您能够…
find_in_set
在查询表中数据用 ‘,’隔开的记录时 如图所示 用like 查询auth 字段中含A5是查不出来的 如图所示 用find_in_set 可以查询出 auth字段含 ‘A5’的记录 如图所示 用法介绍： find_in_set 可以查询出字段内容以英文逗号隔开的记录 find_in_set(匹配值,字段名)转载于:h…
一堆棋子java代码编程_网易2018校招内推编程题-堆棋子-C++实现
0 1 3 10解法暴力枚举所有可能的点。如图所示，黑点为输入点。所需遍历的点为红线的交点，红圈表示。当时自己写的是遍历了外围红线所构成的封闭矩形里面所有的点了，只有60%的AC率，原因超时。看了康学长的代码，才知道有些…
linux free命令详解和使用实例(查看内存使用率)
free命令可以显示Linux系统中空闲的、已用的物理内存及swap内存,及被内核使用的buffer。在Linux系统监控的工具中，free命令是最经常使用的命令之一1．命令格式： free [参数] 2．命令功能： free 命令显示系统使用和空闲的…
awk linux 获取端口号_Linux提权后获取敏感信息命令
如果不能执行的可能是不同类型的linux。系统版本?cat /etc/issuecat /etc/*-releasecat /etc/lsb-releasecat /etc/redhat-release内核版本？cat /proc/versionuname -auname -mrsrpm -q kerneldmesg | grep Linuxls /boot | grep vmlinuz环境变量？cat /…
lemp-------3多站点访问，，访问控制，，虚拟目录
基于ip vi /etc/nginx/nginx.conf server { listen 192.168.1.142:80; server_name localhost; access_log logs/host.access.log main; location / { root /web2; index index.html index.htm index.php; } error_page 500 502 503 504 /50x.html; location /50x.html { root…
oracle统计id出现次数,oracle 统计sql
oracle 统计月平均交易次数 ：select n_tsc_src_usr_id , floor(count(c_tsc_no)/trunc(months_between(max(d_tsc_req_time),min(d_tsc_req_time))))from tbl_tsc_basegroup by n_tsc_src_usr_idhaving months_between(max(d_tsc_req_time),min(d_tsc_req_time)) &g…
java避免使用orderby_java – Spring安全配置@Order不是唯一的例外
我试图在我的Spring Security配置中注册多个过滤器,但是我总是得到相同的异常：04-Nov-2015 14:35:23.792 WARNING [RMI TCP Connection(3)-127.0.0.1]org.springframework.web.context.support.AnnotationConfigWebApplicationContext.refreshException encountered…
es 启动问题
max file descriptors [4096] for elasticsearch process is too low, increase to at least [65536] vim /etc/security/limits.conf 文件末尾添加： mst hard nofile 65536 mst soft nofile 65536 mst是es启动用户 max virtual memory areas vm.max_map_count [65…
imrot matlab,Matlabtuxiangpipei
文件名大小更新时间Matlab图像匹配的都可以用的到做三维重建\matlabcode\examples-code\11.jpg.......................................\..........\.............\addons\rectifData.mat.......................................\..........\.............\......\shelves.jp…
mysql 5.7 full_MySQL5.7默认打开ONLY_FULL_GROUP_BY 解决方案
MySQL5.7后将sql_mode的ONLY_FULL_GROUP_BY模式默认设置为打开状态，这样一来，很多之前的sql语句可能会出现错误，错误信息如下：Error Code: 1055. Expression #3 of SELECT list is not in GROUP BY clause and contains nonaggreg…
通过web sql实现增删查改
<!DOCTYPE html><html><head lang"en"> <meta charset"UTF-8"> <title></title></head><body><h3>***添加学生***</h3>学号：<input type"text" id"id&qu…
java发送苹果消息慢_Spring-boot JMS 发送消息慢的解决方法
Spring-boot JMS 发送消息慢的问题解决Servicepublic class Producer {Autowiredprivate JmsMessagingTemplate jmsTemplate;public void sendMessage(Destination destination, final String message){jmsTemplate.convertAndSend(destination, message);}}经使用JMeter进行压…
快速幂 + 矩阵快速幂
快速幂 1 #include<iostream>2 #include<algorithm>3 #include<cstring>4 #define LL long long5 using namespace std;6 7 LL Pow(LL a, LL b)8 {9 LL ans 1; 10 while(b ! 0){ 11 if(b&1) 12 ans * a; 13 a …
【Ctsc2011】幸福路径
题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id2306 给定一张有向图，每个点有权值，蚂蚁从某个节点出发，初始体力值为$1$，每走一条边$体力值*p$，每经过一个点会获得幸福值为$点权*体力值$&…
mysql 去重con_python 爬虫 实现增量去重和定时爬取实例
前言： 在爬虫过程中，我们可能需要重复的爬取同一个网站，为了避免重复的数据存入我们的数据库中 通过实现增量去重 去解决这一问题 本文还针对了那些需要实时更新的网站 增加了一个定时爬取的功能；本文作者同开源中国(殊途同归_)&a…