c - Optimising and why openmp is much slower than sequential way? -

-

i newbie in programming openmp. wrote simple c program multiply matrix vector. unfortunately, comparing executing time found openmp slower sequential way.

here code (here matrix n*n int, vector n int, result n long long):

#pragma omp parallel private(i,j) shared(matrix,vector,result,m_size) for(i=0;i<m_size;i++) {     for(j=0;j<m_size;j++)   {       result[i]+=matrix[i][j]*vector[j];   } }

and code sequential way:

for (i=0;i<m_size;i++)         for(j=0;j<m_size;j++)             result[i] += matrix[i][j] * vector[j];

when tried these 2 implementations 999x999 matrix , 999 vector, execution time is:

sequential: 5439 ms parallel: 11120 ms

i cannot understand why openmp slower sequential algo (over 2 times slower!) can solve problem?

because when openmp distributes work among threads there lot of administration/synchronisation going on ensure values in shared matrix , vector not corrupted somehow. though read-only: humans see easily, compiler may not.

things try out pedagogic reasons:

0) happens if matrix , vector not shared?

1) parallelize inner "j-loop" first, keep outer "i-loop" serial. see happens.

2) not collect sum in result[i], in variable temp , assign contents result[i] after inner loop finished avoid repeated index lookups. don't forget init temp 0 before inner loop starts.


Comments

Post a Comment

Popular posts from this blog

c++ - Function signature as a function template parameter -

-
i'd use function signature template argument. works great classes, when try same trick function templates, msvc throws error: error c2768: 'func' : illegal use of explicit template arguments here's code: template <typename signature> void func(); template <typename r, typename a1> void func<r(a1)>(); what should make work? you cannot partially specialize function template, not supported language. can create partially specialized class template static member function, , possibly trampoline function instantiate class template , invoke static function. something this: namespace detail { template<typename signature> struct helper; template<typename r, typename a1> struct helper<r(a1)> { static void call() { // stuff r , a1... } }; } template<typename signature> void func() { detail::helper<signature>::call(); }
Read more

algorithm - What are some ways to combine a number of (potentially incompatible) sorted sub-sets of a total set into a (partial) ordering of the total set? -

-
i have sorting problem i'm pretty sure doesn't have canonical "answer" - has number of different approaches, each pros/cons. i'm interested in hearing different approaches. suppose bunch of people {p_i | i=1,...,m} ranking favorite ice cream flavors. in total there n different flavors, single person p_i ranks n_i << n flavors or familiar with. i'm interested in combining these sub-rankings plausible overall ranking of n flavors. for concrete situation: m , n both 1000 (by coincidence), , each n_i 20 . can assume there's enough overlap in flavors people rank such no flavor "isolated". again, i'm interested in hearing different ways approach this, if there isn't single clear answer. thanks! one way approach problem construct directed cyclic graph representing of constraints on ordering of elements. nodes in graph represent different elements compared , have edge first node second node if had ranked first e...
Read more

How to call a javascript function after the page loads with a chrome extension? -

-
there's function on page want call after loads. not own page, can't insert html it. can please provide specific examples how chrome extension or greasemonkey? (function submitaction() ). both chrome extensions , greasemonkey offer ways specify when script inject if executed. in chrome extension, you set in manifest : "content_scripts": [ { "run_at" : "document_idle", in greasemonkey compatible userscripts, set in metadata // @run-at document-end in both cases, have add script if want ensure resources loaded when function executed : document.addeventlistener('load', submitaction);
Read more