[OpenMP] Solution 5 : Dyadic Product

This material is part of the course Practical Introduction to Parallel Programming and meant for educational purposes only.

Task

Write a parallel program that computes the dyadic product

\begin{equation*} \mathbf{x}\mathbf{y}^\top= \begin{bmatrix} x_1y_1 & x_1b_2 & \cdots & x_1y_N\\ x_2y_1 & x_2b_2 & \cdots & x_2y_N\\ \vdots & \vdots & \ddots & \vdots\\ x_Ny_1 & x_Nb_2 & \cdots & x_Ny_N \end{bmatrix} \end{equation*}

for two vectors \(\mathbf{x},\mathbf{y}\in\mathbb{R}^N\).

Test your program for different scheduling strategies and chunk sizes, and different values for \(N\).

Cole este comando no seu terminal:

A senha para conectar é

Questão 1: Dyadic Product

#include <stdio.h>

#ifdef _OPENMP
#include <omp.h>
#else
#define omp_get_num_threads() 1
#define omp_get_max_threads() 1
#define omp_get_thread_num()  0
#endif

int main( int argc, char** argv )
{
    const int n=1000;
    int x[n],y[n];
    int d[n][n];

    double t0 = omp_get_wtime();

#pragma omp parallel default(none) shared(x,y,d)
    {
        /* Initialization */
#pragma omp for schedule(static)
        for (int i=0; i<n; i++)
        {
            x[i] = i; y[i]= i ;
        }

        /* Parallel dyadic product */
#pragma omp for schedule(static) collapse(2)
        for (int i=0; i<n; i++)
            for (int j=0; j<n; j++)
            {
                d[i][j] = x[i]*y[j];
            }
    }

    double t1 = omp_get_wtime();

    if (n<10)
    {
        printf("Dyadic product (in %g seconds):\n", t1-t0);
        for (int i=0; i<n; i++)
        {
            for (int j=0; j<n; j++)
                printf("%d ", d[i][j]);
            printf("\n");
        }
    }
    else
        printf("Dyadic product (in %g seconds):\nSet n<10 to produce output\n", t1-t0);

    return 0;
}

Questão 2: OpenMP environment variables

Specify the OpenMP environment variables that should be used for running your program (export OMP_NUM_THREADS=1 ...). Use semicolons to separate multiple variables

Questão 3: Command Line Argument

Specify the command line arguments that should be passed to your program when it is run. Arguments must be given in quotes and separated by commas. Leave this filed empty if you do not want to specify command line arguments.

Example: "1","int","arg=2" is interpreted as three arguments 1, int, and arg=2.

Autores	Erik-Jan Lingen, Matthias Möller
Prazo de entrega	Sem prazo
Limite de submissão	No limitation

Informações

Entrar