\n \n \n \n \n\n
\n \n \n \n
/*<\/span><\/div> * =====================================================================================<\/span><\/div> *       Filename:  hello.c<\/span><\/div> *         Author:  Raimon Grau Cuscó (rg), raimonster@gmail.com<\/span><\/div> * =====================================================================================<\/span><\/div> */<\/span><\/div> <\/div>#include <omp.h><\/span><\/div>#include <stdlib.h><\/span><\/div> <\/div>int<\/span> main<\/span> ()<\/span><\/div>{<\/span><\/div>        double<\/span> start<\/span>,<\/span>end<\/span>,<\/span> r<\/span>,<\/span> sh<\/span>;<\/span><\/div>        srand<\/span>(<\/span>8<\/span>);<\/span><\/div>        //srand(time(0));<\/span><\/div>        sh<\/span>=<\/span>0<\/span>;<\/span><\/div>        start<\/span> =<\/span> omp_get_wtime<\/span>();<\/span><\/div>        printf<\/span>(<\/span>"Total number of procs: %d<\/span>\\n<\/span>"<\/span>,<\/span> omp_get_num_procs<\/span>());<\/span><\/div>        printf<\/span>(<\/span>"Total number of max available threads: %d<\/span>\\n<\/span>"<\/span>,<\/span> omp_get_max_threads<\/span>());<\/span><\/div>#pragma omp parallel shared(sh)<\/span><\/div>{<\/span><\/div>        printf<\/span>(<\/span>"Hello world! %d<\/span>\\n<\/span>"<\/span>,<\/span> omp_get_thread_num<\/span>());<\/span><\/div>        sh<\/span>+=<\/span>rand<\/span>();<\/span><\/div>        printf<\/span>(<\/span>"%f<\/span>\\n<\/span>"<\/span>,<\/span> sh<\/span>);<\/span><\/div>}<\/span><\/div>        end<\/span> =<\/span> omp_get_wtime<\/span>();<\/span><\/div>        printf<\/span>(<\/span>"Elapsed time: %f<\/span>\\n<\/span>"<\/span>,<\/span> end<\/span> -<\/span> start<\/span>);<\/span><\/div>        printf<\/span>(<\/span>"acc rand: %f<\/span>\\n<\/span>"<\/span>,<\/span>sh<\/span>);<\/span><\/div> <\/div>        return<\/span> 0<\/span>;<\/span><\/div>}<\/span><\/div> <\/div><\/pre><\/div>\n          \n        <\/div>\n\n        \n          view raw<\/a>\n          hello.c<\/a>\n          This Gist<\/a> brought to you by GitHub<\/a>.\n        <\/div>\n      <\/div>\n    \n            \n\n      \n        \n          \n          \n          \n            #include <stdio.h><\/span><\/div>#include <stdlib.h><\/span><\/div>#include <unistd.h><\/span><\/div>#include <math.h><\/span><\/div>#include <omp.h><\/span><\/div> <\/div>#define MATSIZE 700<\/span><\/div> <\/div>#define DEBUG<\/span><\/div> <\/div>void<\/span> init_mat<\/span> (<\/span>double<\/span> M<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>],<\/span> double<\/span> v<\/span>)<\/span><\/div>{<\/span><\/div>   int<\/span> i<\/span>,<\/span> j<\/span>;<\/span><\/div> <\/div>   #pragma omp for collapse(2)<\/span><\/div>   for<\/span> (<\/span>i<\/span>=<\/span>0<\/span>;<\/span> i<\/span><<\/span>MATSIZE<\/span>;<\/span> i<\/span>++<\/span>)<\/span> {<\/span><\/div>      for<\/span> (<\/span>j<\/span>=<\/span>0<\/span>;<\/span> j<\/span><<\/span>MATSIZE<\/span>;<\/span> j<\/span>++<\/span>)<\/span> {<\/span><\/div>         M<\/span>[<\/span>i<\/span>][<\/span>j<\/span>]<\/span> =<\/span> sin<\/span>((<\/span>double<\/span>)<\/span> i<\/span>);<\/span><\/div>      }<\/span><\/div>   }<\/span><\/div>}<\/span><\/div> <\/div>void<\/span> matmul<\/span> (<\/span>double<\/span> C<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>],<\/span><\/div>              double<\/span> A<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>],<\/span> double<\/span> B<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>])<\/span><\/div>{<\/span><\/div>   int<\/span> i<\/span>,<\/span> j<\/span>,<\/span> k<\/span>;<\/span><\/div>   double<\/span> aux<\/span>;<\/span><\/div> <\/div>   #pragma omp for private(j,k,aux)<\/span><\/div>   for<\/span> (<\/span>i<\/span>=<\/span>0<\/span>;<\/span> i<\/span><<\/span>MATSIZE<\/span>;<\/span> i<\/span>++<\/span>)<\/span> {<\/span><\/div>      for<\/span> (<\/span>j<\/span>=<\/span>0<\/span>;<\/span> j<\/span><<\/span>MATSIZE<\/span>;<\/span> j<\/span>++<\/span>)<\/span> {<\/span><\/div>         aux<\/span> =<\/span> 0<\/span>;<\/span><\/div>         for<\/span> (<\/span>k<\/span>=<\/span>0<\/span>;<\/span> k<\/span><<\/span>MATSIZE<\/span>;<\/span> k<\/span>++<\/span>)<\/span> {<\/span><\/div>            aux<\/span> +=<\/span> A<\/span>[<\/span>i<\/span>][<\/span>k<\/span>]<\/span> *<\/span>B<\/span>[<\/span>k<\/span>][<\/span>j<\/span>];<\/span><\/div>         }<\/span><\/div>        C<\/span>[<\/span>i<\/span>][<\/span>j<\/span>]<\/span> +=<\/span> aux<\/span>;<\/span><\/div>      }<\/span><\/div>   }<\/span><\/div>}<\/span><\/div> <\/div>int<\/span> comp<\/span> (<\/span>double<\/span> A<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>],<\/span> double<\/span> B<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>])<\/span><\/div>{<\/span><\/div>   int<\/span> i<\/span>,<\/span> j<\/span>,<\/span> count<\/span> =<\/span> 0<\/span>;<\/span><\/div> <\/div>   for<\/span> (<\/span>i<\/span>=<\/span>0<\/span>;<\/span> i<\/span><<\/span>MATSIZE<\/span>;<\/span> i<\/span>++<\/span>)<\/span> {<\/span><\/div>      for<\/span> (<\/span>j<\/span>=<\/span>0<\/span>;<\/span> j<\/span><<\/span>MATSIZE<\/span>;<\/span> j<\/span>++<\/span>)<\/span> {<\/span><\/div>         if<\/span> (<\/span>abs<\/span> (<\/span>A<\/span>[<\/span>i<\/span>][<\/span>j<\/span>]<\/span>-<\/span> B<\/span>[<\/span>i<\/span>][<\/span>j<\/span>])<\/span> ><\/span> 1e-10<\/span>)<\/span> count<\/span>++<\/span>;<\/span><\/div>      }<\/span><\/div>   }<\/span><\/div>   return<\/span> count<\/span>;<\/span><\/div>}<\/span><\/div> <\/div>int<\/span> main<\/span> (<\/span>int<\/span> argc<\/span>,<\/span> char<\/span> *<\/span> argv<\/span>[])<\/span><\/div>{<\/span><\/div>   double<\/span> start<\/span>,<\/span> end<\/span>;<\/span><\/div>   double<\/span> A<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>];<\/span><\/div>   double<\/span> B<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>];<\/span><\/div>   double<\/span> C<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>];<\/span><\/div> <\/div>   init_mat<\/span> (<\/span>A<\/span>,<\/span> 0.1<\/span>);<\/span><\/div>   init_mat<\/span> (<\/span>B<\/span>,<\/span> 0.1<\/span>);<\/span><\/div>   init_mat<\/span> (<\/span>C<\/span>,<\/span> 0.0<\/span>);<\/span><\/div> <\/div>   start<\/span> =<\/span> omp_get_wtime<\/span>();<\/span><\/div>   matmul<\/span> (<\/span>C<\/span>,<\/span> A<\/span>,<\/span> B<\/span>);<\/span><\/div>   end<\/span> =<\/span> omp_get_wtime<\/span>();<\/span><\/div>   printf<\/span> (<\/span>"matmul1 time: %f seconds<\/span>\\n<\/span>"<\/span>,<\/span> end<\/span>-<\/span>start<\/span>);<\/span><\/div> <\/div>#ifdef DEBUG<\/span><\/div>   double<\/span> D<\/span>[<\/span>MATSIZE<\/span>][<\/span>MATSIZE<\/span>];<\/span><\/div>   init_mat<\/span> (<\/span>D<\/span>,<\/span> 0.0<\/span>);<\/span><\/div>   omp_set_num_threads<\/span>(<\/span>1<\/span>);<\/span><\/div>   matmul<\/span> (<\/span>D<\/span>,<\/span>A<\/span>,<\/span>B<\/span>);<\/span><\/div>   printf<\/span> (<\/span>"Comparing...<\/span>\\n<\/span>"<\/span>);<\/span><\/div>   if<\/span> (<\/span>comp<\/span> (<\/span>C<\/span>,<\/span> D<\/span>))<\/span> printf<\/span> (<\/span>"matmul unsuccessful<\/span>\\n<\/span>"<\/span>);<\/span><\/div>   else<\/span>  printf<\/span> (<\/span>"matmul ok<\/span>\\n<\/span>"<\/span>);<\/span><\/div>#endif<\/span><\/div> <\/div>}<\/span><\/div><\/pre><\/div>\n          \n        <\/div>\n\n        \n          view raw<\/a>\n          matmul.c<\/a>\n          This Gist<\/a> brought to you by GitHub<\/a>.\n        <\/div>\n      <\/div>\n    \n  \n<\/div>\n')