blitzpp / blitz

Blitz++ Multi-Dimensional Array Library for C++
https://github.com/blitzpp/blitz/wiki
Other
404 stars 83 forks source link

[efficient matrix multiplication] with only blitz-inbuilt-features #118

Closed ClmnsRck closed 5 years ago

ClmnsRck commented 5 years ago

The title says it all, but what is the fastest way to do Matrix-multiplication on two blitz::Array's?

A = sum(B(i,j) * C(j,k), j);

This is terribly slow.

citibeth commented 5 years ago

Write a loop? Use Eigen::Map<Eigen::Matrix<>> on it?

On Tue, Apr 23, 2019 at 05:26 ClmnsRck notifications@github.com wrote:

The title says it all, but what is the fastetst way to do Matrix-multiplication on two blitz::Array's?

A = sum(B(i,j) * C(j,k), j);

This is terribly slow.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/blitzpp/blitz/issues/118, or mute the thread https://github.com/notifications/unsubscribe-auth/AAOVY52YPEKRY4VR5CPZAUDPR3I5NANCNFSM4HHWFDVA .

ClmnsRck commented 5 years ago

Well then it would not be using only blitz-built-in features, would it? i'll try the loop ...

citibeth commented 5 years ago

Not using blot built in features doesn’t bother me. Blitz is good at storage, Eigen is good at algorithms.

On Tue, Apr 23, 2019 at 06:28 ClmnsRck notifications@github.com wrote:

Well then it would not be using only blitz-built-in features, would it? i'll try the loop ...

— You are receiving this because you commented.

Reply to this email directly, view it on GitHub https://github.com/blitzpp/blitz/issues/118#issuecomment-485747050, or mute the thread https://github.com/notifications/unsubscribe-auth/AAOVY57DJLBIOYRCBTMVIGTPR3QCZANCNFSM4HHWFDVA .