Tugas Arsitektur Komputer Paralel4

TUGAS ARSITEKTUR KOMPUTER
PARALEL
EDGE DETECTION IN OPENMP USING C
LANGUAGE
Nama
: Grasiadi Hersanto
NIM
: 13212124
New.pgm
Original.pgm
New.pgm
Original.pgm
Berikut besar alokasi memory beserta lama eksekusi menggunakan

OpenMP:
Berikut besar alokasi memory beserta lama eksekusi menggunakan

OpenMP:
Proses
OpenM
P
Waktu Eksekusi
0.257 s
Alokasi Memori
992k
Serial
1.59s
892k
/* EDGE DETECTION PROGRAM WITH OPENMP

*/
/
****************************************************************************
**/
/
****************************************************************************
**/
/* Name : Grasiadi Hersanto
*/
/* NIM : 13212124
*/
/* EXCERSICE : 4
*/
/* FILE NAME : edge_detection_openmp.c
*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <omp.h>
#define IM_SIZE 256 // maximum image size
#define FILTER_SIZE 3 // filter size
double t1,t2,result;
int tid, nthreads, chunk=10;
/* Declare all function */
void read_pict(int pict[][IM_SIZE], int *r, int *c);
void image_filter(int pict[][IM_SIZE], int r, int c, double filter[]
[FILTER_SIZE],
int new_pict[][IM_SIZE]);
void write_pict(int pict[][IM_SIZE], int r, int c);
/* Main Function */
int main( int argc, char *argv[] )
{
int width, height; // actual image size
int pict[IM_SIZE][IM_SIZE];
int new_pict[IM_SIZE][IM_SIZE];
double flt[FILTER_SIZE][FILTER_SIZE]={{-1.25, 0, -1.25},
{0, 10, 0},
{-1.25, 0, -1.25}};
read_pict(pict, &height, &width);
image_filter(pict, height, width, flt, new_pict);
write_pict(new_pict, height, width);
return(0);
/* End Main Function */
}
/* Begin read_pict function

*/
/
****************************************************************************
**/
/* Purpose : This function reads the image from original.pgm
*/
/
****************************************************************************
**/
/* Variable Definitions
*/
/* Variable Name
Type
Description
*/
/* pict[][]
int
array address
*/
/* r
int *
# of rows pointer
*/
/* c
int *
# of column pointer
*/
/* max
int
maximum pixel value
*/
/* i
int
loop counter
*/
/* j
int
loop counter
*/
/* line[200]
char
string
*/
/* in
FILE *
input file pointer
*/
/
****************************************************************************
**/
void read_pict(int pict[][IM_SIZE], int *r, int *c)
{
int i, j, max;
FILE *in;
char line[200];
in=fopen("original.pgm", "r"); //get image
if(in == NULL)
{
printf("Error reading original.pgm\n"); //Error reading image
exit(1);
}
fgets(line, 199, in); //get PGM Type
line[strlen(line)-1]='\0';
if(strcmp(line, "P2") != 0)
{
printf("Cannot process %s PGM format, only P2 type\n", line);

exit(2);
}
fgets(line, 199, in);
fscanf(in, "%d %d", c, r); //get size width x height
fscanf(in, "%d", &max); //get maximum pixel value
for(i = 0; i < *r; i++)
for(j = 0; j < *c; j++)
fscanf(in, "%d", &pict[i][j]);
fclose(in);
return;
/* End read_pict function */
}
/* Begin image_filter function
*/
/
****************************************************************************
**/
/* Purpose : This function filter the image and create a new image
*/
/
****************************************************************************
**/
*/
/* Variable Name
Type
Description
*/
/* pict[][]
int
array address
*/
/* new_pict[][]
int
array address
*/
/* r
int
# of rows
*/
/* c
int
# of column
*/
/* i
int
loop counter
*/
/* j
int
loop counter
*/
/* m
int
loop counter
*/
/* n
int
loop counter
*/
/* coeff
double
coefficient
*/
/* sum
double
sum
*/
/
****************************************************************************
**/
/* Source Code:
*/
void image_filter(int pict[][IM_SIZE], int r, int c, double filter[]
[FILTER_SIZE],
int new_pict[][IM_SIZE])
{
double coeff = 0;
double sum;
int i, j, m, n;
#pragma omp parallel shared(nthreads,chunk) private(tid,i,j,m,n,sum)
{
tid = omp_get_thread_num();
t1 = omp_get_wtime(); //start time
#pragma omp for schedule (static, chunk)
for(i = 0; i < r; i++)
{
new_pict[i][0] = pict[i][0];
new_pict[i][c-1] = pict[i][c-1];
}
for(j = 0; j < c; j++)
{
new_pict[0][j] = pict[0][j];
new_pict[r-1][j] = pict[r-1][j];
}
for(i = 0; i < FILTER_SIZE; i++)
for(j = 0; j < FILTER_SIZE; j++)
coeff += filter[i][j];
/* Image Filtering */
for(i = 1; i < r-1; i++)
for(j = 1; j < c-1; j++)
{
sum = 0;
for(m = 0; m < FILTER_SIZE; m++)
for(n = 0; n < FILTER_SIZE; n++)
sum += pict[i+(m-1)][j+(n-1)]*filter[m][n];
new_pict[i][j] = (int)sum;
}
if(coeff != 0)
{
for(i = 1; i < r-1; i++)
for(j = 1; j < c-1; j++)

new_pict[i][j] = (int)(new_pict[i][j]/coeff);
}
/* Pixel Checking */
for(i = 1; i < r-1; i++)
for(j = 1; j < c-1; j++)
{
if(new_pict[i][j] < 0)
new_pict[i][j] = 0;
else if(new_pict[i][j] > 255)
new_pict[i][j] = 255;
}
}
return;
/* End image_filter function */
}
/* Begin write_pict function
*/
/
****************************************************************************
**/
/* Purpose : This function write the filtered image to new.pgm
*/
/
****************************************************************************
**/
*/
/* Variable Name
Type
Description
*/
/* pict[][]
int
array address
*/
/* new_pict[][]
int
array address
*/
/* r
int
# of rows
*/
/* c
int
# of column
*/
/* i
int
loop counter
*/
/* j
int
loop counter
*/
/* out
FILE *
output FILE pointer
*/
/
****************************************************************************
**/
/* Source Code:
*/
void write_pict(int pict[][IM_SIZE], int r, int c)
{
int i, j;
FILE *out;
out=fopen("new.pgm", "w");
fprintf(out, "P2\n");
fprintf(out, "# new.pgm\n");
fprintf(out, "%d %d\n", r, c);
fprintf(out, "255\n");
for(i = 0; i < r; i++)
{
for(j = 0; j < c; j++)
fprintf(out, "%5d", pict[i][j]);
fprintf(out, "\n");
}
fclose(out);
return;
/* End write_pict function
*/
}

Tugas Arsitektur Komputer Paralel4

Uploaded by

Copyright:

Available Formats

Tugas Arsitektur Komputer Paralel4

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Tugas Arsitektur Komputer Paralel4

Uploaded by

Copyright:

Available Formats

TUGAS ARSITEKTUR KOMPUTER

Berikut besar alokasi memory beserta lama eksekusi menggunakan

Berikut besar alokasi memory beserta lama eksekusi menggunakan

/* EDGE DETECTION PROGRAM WITH OPENMP

/* Begin read_pict function

printf("Cannot process %s PGM format, only P2 type\n", line);

for(j = 1; j < c-1; j++)

You might also like