performance/performanceNumcpp.cpp · 淘金小子/NumCpp

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

扫描微信二维码支付

取消

支付完成

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

performanceNumcpp.cpp

performanceNumcpp.cpp 1.60 KB

#include <numcpp/core.h>
#include <numcpp/base.h>
#include <numcpp/expressions.h>

using namespace numcpp;

size_t performance01Numcpp(size_t N, size_t M)
{

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;
 }

return toc(t, false);
}

size_t performance02Numcpp(size_t N, size_t M)
{
 Array<float> x = randn<float>(N);
 Array<float> A = randn<float>(N,N);
 Array<float> y(N);

auto t = tic();

//EnableCaching();

for(auto i=0; i<M; i++)
 {
 y += dot(A,x)*i + dot(A,x)*(1.0/(i+1)); // + dot(A,ones(N))*i;
 }

//DisableCaching();

return toc(t, false);
}

size_t performance03Numcpp(size_t N, size_t M)
{
 Array<float> B = randn<float>(N,N);
 Array<float> A = randn<float>(N,N);
 Array<float> C(N,N);

auto t = tic();

//EnableCaching();

for(auto i=0; i<M; i++)
 {
 C += dot(A,B); // + dot(A,ones(N))*i;
 }

//DisableCaching();

return toc(t, false);
}

int main(int argc, char* argv[])
{
 int l = 6;
 int code = 0;
 if(argc > 1)
 l = atoi(argv[1]);
 if(argc > 2)
 code = atoi(argv[2]);

auto N = pow(2,l);
 auto M = N;
 size_t time;

switch (code) {
 case 1:
 time = performance01Numcpp(N,10);
 break;
 case 2:
 time = performance02Numcpp(N,M);
 break;
 case 3:
 time = performance03Numcpp(N,M);
 break;
 }

std::cout << time << std::endl;

return 0;
}

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

Tobias Knopp 提交于 2013年10月09日 13:55 +08:00 . more template instantiations

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 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103

#include <numcpp/core.h>
#include <numcpp/base.h>
#include <numcpp/expressions.h>

using namespace numcpp;

size_t performance01Numcpp(size_t N, size_t M)
{

 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;
 }

 return toc(t, false);
}

size_t performance02Numcpp(size_t N, size_t M)
{
 Array<float> x = randn<float>(N);
 Array<float> A = randn<float>(N,N);
 Array<float> y(N);


 auto t = tic();

 //EnableCaching();

 for(auto i=0; i<M; i++)
 {
 y += dot(A,x)*i + dot(A,x)*(1.0/(i+1)); // + dot(A,ones(N))*i;
 }

 //DisableCaching();

 return toc(t, false);
}


size_t performance03Numcpp(size_t N, size_t M)
{
 Array<float> B = randn<float>(N,N);
 Array<float> A = randn<float>(N,N);
 Array<float> C(N,N);


 auto t = tic();

 //EnableCaching();

 for(auto i=0; i<M; i++)
 {
 C += dot(A,B); // + dot(A,ones(N))*i;
 }

 //DisableCaching();

 return toc(t, false);
}

int main(int argc, char* argv[])
{
 int l = 6;
 int code = 0;
 if(argc > 1)
 l = atoi(argv[1]);
 if(argc > 2)
 code = atoi(argv[2]);

 auto N = pow(2,l);
 auto M = N;
 size_t time;


 switch (code) {
 case 1:
 time = performance01Numcpp(N,10);
 break;
 case 2:
 time = performance02Numcpp(N,M);
 break;
 case 3:
 time = performance03Numcpp(N,M);
 break;
 }

 std::cout << time << std::endl;

 return 0;
}