performance/performanceBlitz.cpp · Job/NumCpp

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

扫描微信二维码支付

取消

支付完成

richgiteeai

Watch

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

1 Star 0 Fork 16

Job/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

performanceBlitz.cpp

performanceBlitz.cpp 1.70 KB

#include <blitz/array.h>
#include <random/normal.h>
#include <numcpp/core/utils.h>

using namespace blitz;

size_t performance01Blitz(size_t N, size_t M)
{

auto t = numcpp::tic();

Array<double,1> linspace(N);

linspace = 56;

for(int l=0;l<10;l++)
 {
 Array<std::complex<double>,1> x(N);
 x = 0;
 for(auto i=0; i<M; i++)
 x += 1;

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

for(auto i=0; i<M; i++)
 x += ((double)i)*exp(2*M_PI*i*linspace * std::complex<double>(0,1));

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

return numcpp::toc(t, false);
}

size_t performance02Blitz(size_t N, size_t M)
{
 Array<float,1> x(N);
 Array<float,2> A(N,N);
 Array<float,1> y(N);

x =1; A=2;

firstIndex k;
 secondIndex j;

auto t = numcpp::tic();

for(auto i=0; i<M; i++)
 {
 y += sum(A(j,k)*x(k),k) + sum(A(j,k)*x(k),k)*(1.0/(i+1));
 }

return numcpp::toc(t, false);
}

size_t performance03Blitz(size_t N, size_t M)
{
 Array<float,2> B(N,N);
 Array<float,2> A(N,N);
 Array<float,2> C(N,N);

B=1; A=2;

firstIndex k;
 secondIndex j;
 thirdIndex l;

auto t = numcpp::tic();

for(auto i=0; i<M; i++)
 {
 C += sum(A(j,k)*B(k,j),k)*i;
 }

return numcpp::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 = performance01Blitz(N,10);
 break;
 case 2:
 time = performance02Blitz(N,M);
 break;
 case 3:
 time = performance03Blitz(N,M);
 break;
 }

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

return 0;
}

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

Tobias Knopp 提交于 2013年04月15日 14:50 +08:00 . blitz tests

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 104 105 106 107 108 109 110

#include <blitz/array.h>
#include <random/normal.h>
#include <numcpp/core/utils.h>

using namespace blitz;

size_t performance01Blitz(size_t N, size_t M)
{

 auto t = numcpp::tic();

 Array<double,1> linspace(N);

 linspace = 56;


 for(int l=0;l<10;l++)
 {
 Array<std::complex<double>,1> x(N);
 x = 0;
 for(auto i=0; i<M; i++)
 x += 1;

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

 for(auto i=0; i<M; i++)
 x += ((double)i)*exp(2*M_PI*i*linspace * std::complex<double>(0,1));

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

 return numcpp::toc(t, false);
}

size_t performance02Blitz(size_t N, size_t M)
{
 Array<float,1> x(N);
 Array<float,2> A(N,N);
 Array<float,1> y(N);

 x =1; A=2;

 firstIndex k;
 secondIndex j;

 auto t = numcpp::tic();

 for(auto i=0; i<M; i++)
 {
 y += sum(A(j,k)*x(k),k) + sum(A(j,k)*x(k),k)*(1.0/(i+1));
 }

 return numcpp::toc(t, false);
}


size_t performance03Blitz(size_t N, size_t M)
{
 Array<float,2> B(N,N);
 Array<float,2> A(N,N);
 Array<float,2> C(N,N);

 B=1; A=2;

 firstIndex k;
 secondIndex j;
 thirdIndex l;

 auto t = numcpp::tic();

 for(auto i=0; i<M; i++)
 {
 C += sum(A(j,k)*B(k,j),k)*i;
 }

 return numcpp::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 = performance01Blitz(N,10);
 break;
 case 2:
 time = performance02Blitz(N,M);
 break;
 case 3:
 time = performance03Blitz(N,M);
 break;
 }

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

 return 0;
}