performance/performance01.cpp · 一张弓/NumCpp

【升级通知】Gitee 将于 7 月 11 日晚 23 点进行全站升级维护

代码拉取完成,页面将自动刷新

richgiteeai

Watch

不关注关注所有动态仅关注版本发行动态关注但不提醒动态

1 Star 0 Fork 16

一张弓/NumCpp

forked from Gitee 极速下载/NumCpp

代码 Wiki 统计流水线

服务

加入 Gitee

与超过 1400万开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)

免费加入

已有帐号? 立即登录

文件

master

分支 (2)

管理

master

gh-pages

克隆/下载

HTTPS SSH SVN SVN+SSH 下载ZIP

提示

下载代码请复制以下命令到终端执行

为确保你提交的代码身份被 Gitee 正确识别,请执行以下命令完成配置

git config --global user.name userName 
git config --global user.email userEmail

初次使用 SSH 协议进行代码克隆、推送等操作时,需按下述提示完成 SSH 配置

1 生成 RSA 密钥

2 获取 RSA 公钥内容,并配置到 SSH公钥中

在 Gitee 上使用 SVN,请访问使用指南

使用 HTTPS 协议时,命令行会出现如下账号密码验证步骤。基于安全考虑,Gitee 建议配置并使用私人令牌替代登录密码进行克隆、推送等操作

Username for 'https://gitee.com': userName

Password for 'https://userName@gitee.com': # 私人令牌

分支 2

标签 0

numcpp

performance

performance01.cpp

performance01.cpp 843 Bytes

一键复制编辑原始数据按行查看历史

Tobias Knopp 提交于 2013年08月04日 01:18 +08:00 . rewrite of numcpp

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56

#include <numcpp.h>
#include <Eigen/Core>

using namespace numcpp;

int main()
{
 auto N = 1000000;
 auto M = 10;

 // NumCpp

 auto t = tic();

 for(int l=0;l<10;l++)
 {
 Array<cdouble> x = zeros(N);
 for(auto i=0; i<M; i++)
 x += ones(N);

 for(auto i=0; i<N; i++)
 x(i) += i;

 for(auto i=0; i<M; i++)
 x += i*exp(2*pi*I*i*linspace(0,1,N));

 for(auto i=0; i<M; i++)
 x += i;
 }

 toc(t);

 // Eigen

 t = tic();

 for(int l=0;l<10;l++)
 {
 Eigen::ArrayXcd x = Eigen::ArrayXcd::Zero(N);
 for(auto i=0; i<M; i++)
 x += Eigen::ArrayXcd::Ones(N);

 for(auto i=0; i<N; i++)
 x(i) += i;

 for(auto i=0; i<M; i++)
 x += x += i*exp(std::complex<double>(0,1)*Eigen::ArrayXd::LinSpaced(N,0,1)*i*2*pi);

 for(auto i=0; i<M; i++)
 x += i;
 }

 toc(t);

 return 0;
}