為什么MATLAB在矩陣乘法中速度這么快？我正在使用CUDA、C+、C#和Java編寫一些基準(zhǔn)，并使用MATLAB進(jìn)行驗(yàn)證和矩陣生成。但當(dāng)我用MATLAB進(jìn)行乘法時(shí)，2048x2048甚至更大的矩陣幾乎立即被乘以。 1024x1024 2048x2048 4096x4096
--------- --------- ---------
CUDA C (ms) 43.11 391.05 3407.99
C++ (ms) 6137.10 64369.29 551390.93
C# (ms) 10509.00 300684.00 2527250.00
Java (ms) 9149.90 92562.28 838357.94
MATLAB (ms) 75.01 423.10 3133.90只有CUDA是有競(jìng)爭(zhēng)力的，但我認(rèn)為至少C+會(huì)有點(diǎn)接近，而不是60x慢點(diǎn)。所以我的問(wèn)題是-MATLAB是怎么做到這么快的？C+代碼：float temp = 0;timer.start();for(int j = 0; j < rozmer; j++){
for (int k = 0; k < rozmer; k++)
{
temp = 0;
for (int m = 0; m < rozmer; m++)
{
temp = temp + matice1[j][m] * matice2[m][k];
}
matice3[j][k] = temp;
}}timer.stop();編輯：我也不知道如何看待C#結(jié)果。算法與C+和Java完全相同，但是有一個(gè)巨大的飛躍2048從…1024?編輯2：最新的MATLAB和4096x4096結(jié)果
3 回答
梵蒂岡之花
TA貢獻(xiàn)1900條經(jīng)驗(yàn) 獲得超5個(gè)贊
>> A = rand(1024); gA = gpuArray(A);% warm up by executing the operations a couple of times, and then:>> tic, C = A * A; tocElapsed time is 0.075396 seconds.>> tic, gC = gA * gA; tocElapsed time is 0.008621 seconds.
gpuArray
使用R2014a更新timeitgputimeit
>> A = rand(1024); gA = gpuArray(A);>> timeit(@()A*A)ans = 0.0324>> gputimeit(@()gA*gA)ans = 0.0022
使用R2018b更新
>> timeit(@()A*A)ans = 0.0229>> gputimeit(@()gA*gA)ans = 4.8019e-04
小唯快跑啊
TA貢獻(xiàn)1863條經(jīng)驗(yàn) 獲得超2個(gè)贊
歷史：
BLAS：
爭(zhēng)取更好的業(yè)績(jī)：
矩陣乘法的技術(shù)細(xì)節(jié)：
dgemm
matice2[m][k]matice2[0][k]matice2[1][k]matice28*1024*1024
timer.start();float temp = 0;//transpose matice2for (int p = 0; p < rozmer; p++){
for (int q = 0; q < rozmer; q++)
{
tempmat[p][q] = matice2[q][p];
}}for(int j = 0; j < rozmer; j++){
for (int k = 0; k < rozmer; k++)
{
temp = 0;
for (int m = 0; m < rozmer; m++)
{
temp = temp + matice1[j][m] * tempmat[k][m];
}
matice3[j][k] = temp;
}}timer.stop();dgemm
慕的地6264312
TA貢獻(xiàn)1817條經(jīng)驗(yàn) 獲得超6個(gè)贊
C=A*B
- 3 回答
- 0 關(guān)注
- 1978 瀏覽
添加回答
舉報(bào)
0/150
提交
取消