#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define NUM_LOOPS 50
int
main(int argc, char *argv[])
{
double **A;
double *u, *v;
long rows = 5000, columns = 5000;
time_t begin_t, end_t;
long i, j, k;
// Inicjujemy wektory i macierz.
u
= (double*)malloc(columns
* sizeof(double));
v
= (double*)malloc(rows
* sizeof(double));
A
= (double**)malloc(rows
* sizeof(double*));
for (i = 0; i < rows; i++)
A
[i
] = (double*)malloc(columns
* sizeof(double));
// Wypełniamy wektor u oraz macierz A dowolnymi wartościami.
for (i = 0; i < columns; i++) {
u[i] = (double)(i / 1000.0f);
for ( j = 0; j < rows; j++ )
A[j][i] = (double)(i * j / 1000.0f);
}
// Pętla zewnętrzna.
for (k = 0; k < NUM_LOOPS; k++) {
// Obliczamy v = A*u z wykorzystaniem dyrektyw OpenMP.
#pragma omp parallel for shared(A, v, u, rows, columns) private(i, j)
for (i = 0; i < rows; i++) {
v[i] = 0.0f;
for (j = 0; j < columns; j++)
v[i] = v[i] + A[i][j] * u[j];
}
}
// Zwalniamy pamięć.
for (i = 0; i < rows; i++)
return 0;
}
Replies to Untitled
Title |
Name |
Language |
UNIX |
When |
Re: Untitled |
Sharp Cassowary |
c |
1423309640 |
9 Years ago. |
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}