Computación Paralela (Sistemas SMP)

Computación Paralela (Sistemas SMP)

Introducción teórica
La práctica consiste el realizar un programa para calcular el valor de π, calculando para ello el valor de una integral definidacuyo resultado sabemos que es π/4. El cálculo del valor de la integral se realizara por el medio del algoritmo de los trapecios y usando programación multihebra. La función es la siguiente:

Elmétodo de los trapecios se usa para el cálculo de integrales definidas, la integral de una función entre dos puntos es el área que esta función encierra entre estos puntos.
El método de los trapeciosdivide esta área en múltiples trapecios, calcula el área de cada trapecio y calcula el resultado final como la suma de todas estas áreas.
Cuando mayor sea el número de trapecios mayor precisiónobtenemos en el resultado.



La figura representa uno de los trapecios en que hemos dividido el segmento de los límites de la integral (a,b). La anchura de este segmento será h y viene dada porla siguiente ecuación donde n es el número de trapecios.

La parte sombreada se toma como el área de este trapecio y viene dada por:

Código fuente
#include
#include
#include#include


#define workers 8
#define trapecios 10000000
#define trapeciosperworker trapecios/workers
#define a 0.0
#define b 1.0

#define check(status,string) if (status != 0) { \errno = status; \ printf(stderr, "%s status %d: %s\n", \ string, status, strerror(status)); \
}
pthread_mutex_t resultado_integral= PTHREAD_MUTEX_INITIALIZER; pthread_mutex_tcond_mutex= PTHREAD_MUTEX_INITIALIZER; pthread_cond_t cond_var= PTHREAD_COND_INITIALIZER;
long double total = 0.0; int thread_hold = 1; int count_fins = 1; pthread_t threads[workers];
staticvoid unlock_cond(void *arg) { int status;
status = pthread_mutex_unlock(&cond_mutex); check(status, "Mutex_unlock");
} double f(double x) {
double return_val = 1 / (1 + x*x);...