diff --git a/src/config.c b/src/config.c
index 7d1c1b7..0c590a4 100644
--- a/src/config.c
+++ b/src/config.c
@@ -3,12 +3,12 @@
 
 typedef struct configuration {
   int n;
-  double north;
-  double east;
-  double south;
-  double west;
-  double init_value;
-  double threshold;
+  float north;
+  float east;
+  float south;
+  float west;
+  float init_value;
+  float threshold;
 } configuration;
 
 int load_config(configuration *config) {
@@ -38,17 +38,17 @@ int load_config(configuration *config) {
     if (strcmp(property, "N") == 0) {
       sscanf(value, "%d", &(config->n));
     } else if (strcmp(property, "NORTH") == 0) {
-      sscanf(value, "%lf", &(config->north));
+      sscanf(value, "%f", &(config->north));
     } else if (strcmp(property, "EAST") == 0) {
-      sscanf(value, "%lf", &(config->east));
+      sscanf(value, "%f", &(config->east));
     } else if (strcmp(property, "SOUTH") == 0) {
-      sscanf(value, "%lf", &(config->south));
+      sscanf(value, "%f", &(config->south));
     } else if (strcmp(property, "WEST") == 0) {
-      sscanf(value, "%lf", &(config->west));
+      sscanf(value, "%f", &(config->west));
     } else if (strcmp(property, "INIT_VALUE") == 0) {
-      sscanf(value, "%lf", &(config->init_value));
+      sscanf(value, "%f", &(config->init_value));
     } else if (strcmp(property, "THRESHOLD") == 0) {
-      sscanf(value, "%lf", &(config->threshold));
+      sscanf(value, "%f", &(config->threshold));
     } else {
       printf("Unknown property %s\n", property);
     }
diff --git a/src/config.h b/src/config.h
index 32635dd..a141c59 100644
--- a/src/config.h
+++ b/src/config.h
@@ -1,11 +1,11 @@
 typedef struct configuration {
   int n;
-  double north;
-  double east;
-  double south;
-  double west;
-  double init_value;
-  double threshold;
+  float north;
+  float east;
+  float south;
+  float west;
+  float init_value;
+  float threshold;
 } configuration;
 
 int load_config(configuration *config);
diff --git a/src/impl/mpi_line.c b/src/impl/mpi_line.c
index 853b43a..d4f05cb 100644
--- a/src/impl/mpi_line.c
+++ b/src/impl/mpi_line.c
@@ -12,12 +12,12 @@
 #define TAG_BORDER 0
 #define TAG_MATRIX 1
 
-double *compute_jacobi(int rank, int numprocs, int n, double init_value, double threshold, borders b, int *iterations) {
-  double *complete_x;
-  double *x;
-  double *new_x;
-  double *tmp_x;
-  double max_diff, global_max_diff, new_value;
+float *compute_jacobi(int rank, int numprocs, int n, float init_value, float threshold, borders b, int *iterations) {
+  float *complete_x;
+  float *x;
+  float *new_x;
+  float *tmp_x;
+  float max_diff, global_max_diff, new_value;
   int i, j;
   int nb = n + 2; // n plus the border
   int rows, rows_to_transmit;
@@ -67,7 +67,7 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
     for (i = 1; i <= rows; i++) {
       for (j = 1; j <= n; j++) {
         new_value = 0.25 * (x[IDX(nb, i - 1, j)] + x[IDX(nb, i, j + 1)] + x[IDX(nb, i + 1, j)] + x[IDX(nb, i, j - 1)]);
-        max_diff = (double) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
         new_x[IDX(nb, i, j)] = new_value;
       }
     }
@@ -77,40 +77,40 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
     if (rank % 2 == 0) {
       if (rank != numprocs - 1) {
         // Send and receive south border
-        MPI_Send(&x[IDX(nb, rows, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD);
-        MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+        MPI_Send(&x[IDX(nb, rows, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD);
+        MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
       }
       if (rank != 0) {
         // Send and receive north border
-        MPI_Send(&x[IDX(nb, 1, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD);
-        MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+        MPI_Send(&x[IDX(nb, 1, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD);
+        MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
       }
     } else {
       // Receive and send north border
-      MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
-      MPI_Send(&x[IDX(nb, 1, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD);
+      MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+      MPI_Send(&x[IDX(nb, 1, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD);
       if (rank != numprocs - 1) {
         // Receive and send south border
-        MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
-        MPI_Send(&x[IDX(nb, rows, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD);
+        MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+        MPI_Send(&x[IDX(nb, rows, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD);
       }
     }
     /* LOG(printf("[Process %d/%d] max_diff: %f\n", rank, numprocs, max_diff)); */
-    MPI_Allreduce(&max_diff, &global_max_diff, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD);
+    MPI_Allreduce(&max_diff, &global_max_diff, 1, MPI_FLOAT, MPI_MAX, MPI_COMM_WORLD);
     /* LOG(printf("[Process %d/%d] global_max_diff: %f\n", rank, numprocs, global_max_diff)); */
     (*iterations)++;
   } while (global_max_diff > threshold);
 
   if (rank == 0) {
     complete_x = create_sa_matrix(nb, nb);
-    memcpy(complete_x, x, (rows + ((rank == numprocs - 1) ? 2 : 1)) * (nb) * sizeof(double));
+    memcpy(complete_x, x, (rows + ((rank == numprocs - 1) ? 2 : 1)) * (nb) * sizeof(float));
     rows_to_transmit = n / numprocs;
     receive_pos = rows + 1;
     for (i = 1; i < numprocs; i++) {
       if (i == numprocs - 1) {
         rows_to_transmit++;
       }
-      MPI_Recv(&complete_x[IDX(nb, receive_pos, 0)], rows_to_transmit * (nb), MPI_DOUBLE, i, TAG_MATRIX, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+      MPI_Recv(&complete_x[IDX(nb, receive_pos, 0)], rows_to_transmit * (nb), MPI_FLOAT, i, TAG_MATRIX, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
       receive_pos += n / numprocs;
     }
   } else {
@@ -119,7 +119,7 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
     if (rank == numprocs - 1) {
       rows_to_transmit++;
     }
-    MPI_Send(&x[IDX(nb, 1, 0)], rows_to_transmit * (nb), MPI_DOUBLE, 0, TAG_MATRIX, MPI_COMM_WORLD);
+    MPI_Send(&x[IDX(nb, 1, 0)], rows_to_transmit * (nb), MPI_FLOAT, 0, TAG_MATRIX, MPI_COMM_WORLD);
   }
 
   return complete_x;
diff --git a/src/impl/mpi_line_async.c b/src/impl/mpi_line_async.c
index 849c3fe..a7e21f3 100644
--- a/src/impl/mpi_line_async.c
+++ b/src/impl/mpi_line_async.c
@@ -12,12 +12,12 @@
 #define TAG_BORDER 0
 #define TAG_MATRIX 1
 
-double *compute_jacobi(int rank, int numprocs, int n, double init_value, double threshold, borders b, int *iterations) {
-  double *complete_x;
-  double *x;
-  double *new_x;
-  double *tmp_x;
-  double max_diff, global_max_diff, new_value;
+float *compute_jacobi(int rank, int numprocs, int n, float init_value, float threshold, borders b, int *iterations) {
+  float *complete_x;
+  float *x;
+  float *new_x;
+  float *tmp_x;
+  float max_diff, global_max_diff, new_value;
   int i, j;
   int nb = n + 2; // n plus the border
   int rows, rows_to_transmit;
@@ -66,18 +66,18 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
   do {
     if (rank != numprocs - 1) {
       // Send south border
-      MPI_Isend(&x[IDX(nb, rows, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD, &request_south);
+      MPI_Isend(&x[IDX(nb, rows, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD, &request_south);
     }
     if (rank != 0) {
       // Send north border
-      MPI_Isend(&x[IDX(nb, 1, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD, &request_north);
+      MPI_Isend(&x[IDX(nb, 1, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD, &request_north);
     }
     max_diff = 0;
     global_max_diff = 0;
     for (i = 1; i <= rows; i++) {
       for (j = 1; j <= n; j++) {
         new_value = 0.25 * (x[IDX(nb, i - 1, j)] + x[IDX(nb, i, j + 1)] + x[IDX(nb, i + 1, j)] + x[IDX(nb, i, j - 1)]);
-        max_diff = (double) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
         new_x[IDX(nb, i, j)] = new_value;
       }
     }
@@ -86,28 +86,28 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
     x = tmp_x;
     if (rank != numprocs - 1) {
       // Receive south border
-      MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_DOUBLE, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+      MPI_Recv(&x[IDX(nb, rows + 1, 0)], nb, MPI_FLOAT, rank + 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
     }
     if (rank != 0) {
       // Receive north border
-      MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_DOUBLE, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+      MPI_Recv(&x[IDX(nb, 0, 0)], nb, MPI_FLOAT, rank - 1, TAG_BORDER, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
     }
     LOG(printf("[Process %d/%d] max_diff: %f\n", rank, numprocs, max_diff));
-    MPI_Allreduce(&max_diff, &global_max_diff, 1, MPI_DOUBLE, MPI_MAX, MPI_COMM_WORLD);
+    MPI_Allreduce(&max_diff, &global_max_diff, 1, MPI_FLOAT, MPI_MAX, MPI_COMM_WORLD);
     /* LOG(printf("[Process %d/%d] global_max_diff: %f\n", rank, numprocs, global_max_diff)); */
     (*iterations)++;
   } while (global_max_diff > threshold);
 
   if (rank == 0) {
     complete_x = create_sa_matrix(nb, nb);
-    memcpy(complete_x, x, (rows + ((rank == numprocs - 1) ? 2 : 1)) * (nb) * sizeof(double));
+    memcpy(complete_x, x, (rows + ((rank == numprocs - 1) ? 2 : 1)) * (nb) * sizeof(float));
     rows_to_transmit = n / numprocs;
     receive_pos = rows + 1;
     for (i = 1; i < numprocs; i++) {
       if (i == numprocs - 1) {
         rows_to_transmit++;
       }
-      MPI_Recv(&complete_x[IDX(nb, receive_pos, 0)], rows_to_transmit * (nb), MPI_DOUBLE, i, TAG_MATRIX, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
+      MPI_Recv(&complete_x[IDX(nb, receive_pos, 0)], rows_to_transmit * (nb), MPI_FLOAT, i, TAG_MATRIX, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
       receive_pos += n / numprocs;
     }
   } else {
@@ -116,7 +116,7 @@ double *compute_jacobi(int rank, int numprocs, int n, double init_value, double
     if (rank == numprocs - 1) {
       rows_to_transmit++;
     }
-    MPI_Send(&x[IDX(nb, 1, 0)], rows_to_transmit * (nb), MPI_DOUBLE, 0, TAG_MATRIX, MPI_COMM_WORLD);
+    MPI_Send(&x[IDX(nb, 1, 0)], rows_to_transmit * (nb), MPI_FLOAT, 0, TAG_MATRIX, MPI_COMM_WORLD);
   }
 
   return complete_x;
diff --git a/src/impl/omp.c b/src/impl/omp.c
index 5888122..b88fa14 100644
--- a/src/impl/omp.c
+++ b/src/impl/omp.c
@@ -8,11 +8,11 @@
 #include "../config.h"
 #include "../utils.h"
 
-double *compute_jacobi(int n, double init_value, double threshold, borders b, int *iterations) {
-  double *x;
-  double *new_x;
-  double *tmp_x;
-  double max_diff, new_value;
+float *compute_jacobi(int n, float init_value, float threshold, borders b, int *iterations) {
+  float *x;
+  float *new_x;
+  float *tmp_x;
+  float max_diff, new_value;
   int i, j;
   int nb = n + 2; // n plus the border
 
@@ -45,7 +45,7 @@ double *compute_jacobi(int n, double init_value, double threshold, borders b, in
     for (i = 1; i <= n; i++) {
       for (j = 1; j <= n; j++) {
         new_value = 0.25 * (x[IDX(nb, i - 1, j)] + x[IDX(nb, i, j + 1)] + x[IDX(nb, i + 1, j)] + x[IDX(nb, i, j - 1)]);
-        max_diff = (double) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
         new_x[IDX(nb, i, j)] = new_value;
       }
     }
diff --git a/src/impl/sequential.c b/src/impl/sequential.c
index 2767447..81b8259 100644
--- a/src/impl/sequential.c
+++ b/src/impl/sequential.c
@@ -8,11 +8,11 @@
 #include "../config.h"
 #include "../utils.h"
 
-double *compute_jacobi(int n, double init_value, double threshold, borders b, int *iterations) {
-  double *x;
-  double *new_x;
-  double *tmp_x;
-  double max_diff, new_value;
+float *compute_jacobi(int n, float init_value, float threshold, borders b, int *iterations) {
+  float *x;
+  float *new_x;
+  float *tmp_x;
+  float max_diff, new_value;
   int i, j;
   int nb = n + 2; // n plus the border
 
@@ -42,7 +42,7 @@ double *compute_jacobi(int n, double init_value, double threshold, borders b, in
     for (i = 1; i <= n; i++) {
       for (j = 1; j <= n; j++) {
         new_value = 0.25 * (x[IDX(nb, i - 1, j)] + x[IDX(nb, i, j + 1)] + x[IDX(nb, i + 1, j)] + x[IDX(nb, i, j - 1)]);
-        max_diff = (double) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
         new_x[IDX(nb, i, j)] = new_value;
       }
     }
diff --git a/src/impl/sse.c b/src/impl/sse.c
index b16ca8a..550bdb0 100644
--- a/src/impl/sse.c
+++ b/src/impl/sse.c
@@ -8,14 +8,15 @@
 #include "../config.h"
 #include "../utils.h"
 
-double *compute_jacobi(int n, double init_value, double threshold, borders b, int *iterations) {
-  double *x;
-  double *new_x;
-  double *tmp_x;
-  double max_diff, new_value;
+float *compute_jacobi(int n, float init_value, float threshold, borders b, int *iterations) {
+  float *x;
+  float *new_x;
+  float *tmp_x;
+  float max_diff, new_value;
   int i, j;
   int nb = n + 2; // n plus the border
-  int n_mult = (n % 2 == 0) ? n : n - 1;
+  int n_mult = (n / 4) * 4;
+  __m128 new_value_vec, tmp_vec;
   /* Initialize boundary regions */
   x = create_sa_matrix(nb, nb);
   new_x = create_sa_matrix(nb, nb);
@@ -40,24 +41,25 @@ double *compute_jacobi(int n, double init_value, double threshold, borders b, in
   do {
     max_diff = 0;
     for (i = 1; i <= n; i++) {
-      for (j = 1; j <= n_mult; j += 2) {
-        __m128d new_value_vec, tmp_vec;
-        new_value_vec = _mm_loadu_pd(&x[IDX(nb, i - 1, j)]);
-        tmp_vec = _mm_loadu_pd(&x[IDX(nb, i + 1, j)]);
-        new_value_vec = _mm_add_pd(new_value_vec, tmp_vec);
-        tmp_vec = _mm_loadu_pd(&x[IDX(nb, i, j - 1)]);
-        new_value_vec = _mm_add_pd(new_value_vec, tmp_vec);
-        tmp_vec = _mm_loadu_pd(&x[IDX(nb, i, j + 1)]);
-        new_value_vec = _mm_add_pd(new_value_vec, tmp_vec);
-        tmp_vec = _mm_set1_pd(0.25);
-        new_value_vec = _mm_mul_pd(new_value_vec, tmp_vec);
-        _mm_storeu_pd(&new_x[IDX(nb, i, j)], new_value_vec);
-        max_diff = (double) fmax(max_diff, fabs(new_x[IDX(nb, i, j)] - x[IDX(nb, i, j)]));
-        max_diff = (double) fmax(max_diff, fabs(new_x[IDX(nb, i, j + 1)] - x[IDX(nb, i, j + 1)]));
+      for (j = 1; j <= n_mult; j += 4) {
+        new_value_vec = _mm_loadu_ps(&x[IDX(nb, i - 1, j)]);
+        tmp_vec = _mm_loadu_ps(&x[IDX(nb, i + 1, j)]);
+        new_value_vec = _mm_add_ps(new_value_vec, tmp_vec);
+        tmp_vec = _mm_loadu_ps(&x[IDX(nb, i, j - 1)]);
+        new_value_vec = _mm_add_ps(new_value_vec, tmp_vec);
+        tmp_vec = _mm_loadu_ps(&x[IDX(nb, i, j + 1)]);
+        new_value_vec = _mm_add_ps(new_value_vec, tmp_vec);
+        tmp_vec = _mm_set1_ps(0.25);
+        new_value_vec = _mm_mul_ps(new_value_vec, tmp_vec);
+        _mm_storeu_ps(&new_x[IDX(nb, i, j)], new_value_vec);
+        max_diff = (float) fmax(max_diff, fabs(new_x[IDX(nb, i, j)] - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_x[IDX(nb, i, j + 1)] - x[IDX(nb, i, j + 1)]));
+        max_diff = (float) fmax(max_diff, fabs(new_x[IDX(nb, i, j + 2)] - x[IDX(nb, i, j + 2)]));
+        max_diff = (float) fmax(max_diff, fabs(new_x[IDX(nb, i, j + 3)] - x[IDX(nb, i, j + 3)]));
       }
       for (j = n_mult; j <= n; j++) {
         new_value = 0.25 * (x[IDX(nb, i - 1, j)] + x[IDX(nb, i, j + 1)] + x[IDX(nb, i + 1, j)] + x[IDX(nb, i, j - 1)]);
-        max_diff = (double) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
+        max_diff = (float) fmax(max_diff, fabs(new_value - x[IDX(nb, i, j)]));
         new_x[IDX(nb, i, j)] = new_value;
       }
     }
diff --git a/src/main/main.c b/src/main/main.c
index 289e1b1..8bf645f 100644
--- a/src/main/main.c
+++ b/src/main/main.c
@@ -5,15 +5,15 @@
 #include "../config.h"
 #include "../utils.h"
 
-double *compute_jacobi(int n, double init_value, double threshold, borders b, int *iterations);
+float *compute_jacobi(int n, float init_value, float threshold, borders b, int *iterations);
 
 int main(int argc, char* argv[]) {
   int n;
-  double init_value, threshold;
+  float init_value, threshold;
   borders b;
   int config_loaded;
   configuration config;
-  double *x;
+  float *x;
   int iterations;
   struct timeval start, end;
   long secs_used, micros_used;
@@ -36,7 +36,7 @@ int main(int argc, char* argv[]) {
 
   secs_used = (end.tv_sec - start.tv_sec);
   micros_used = ((secs_used * 1000000) + end.tv_usec) - (start.tv_usec);
-  printf("Wall clock time: %fs\n", (double) micros_used / 1000000);
+  printf("Wall clock time: %fs\n", (float) micros_used / 1000000);
   printf("Iterations: %d\n", iterations);
   if (n < 10) {
     print_sa_matrix(x, n + 2, n + 2);
diff --git a/src/main/main_mpi.c b/src/main/main_mpi.c
index 6d3961b..3013fa9 100644
--- a/src/main/main_mpi.c
+++ b/src/main/main_mpi.c
@@ -5,18 +5,18 @@
 #include "../config.h"
 #include "../utils.h"
 
-double *compute_jacobi(int rank, int numprocs, int n, double init_value, double threshold, borders b, int *iterations);
+float *compute_jacobi(int rank, int numprocs, int n, float init_value, float threshold, borders b, int *iterations);
 
 int main(int argc, char* argv[]) {
   int rank;
   int numprocs;
   int n;
-  double init_value, threshold;
-  double north, south, east, west;
+  float init_value, threshold;
+  float north, south, east, west;
   borders b;
   int config_loaded;
   configuration config;
-  double *x;
+  float *x;
   double startwtime = 0.0, endwtime;
   int iterations;
 
@@ -39,12 +39,12 @@ int main(int argc, char* argv[]) {
     west = config.west;
   }
   MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&init_value, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&threshold, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&north, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&south, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&east, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
-  MPI_Bcast(&west, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&init_value, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&threshold, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&north, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&south, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&east, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
+  MPI_Bcast(&west, 1, MPI_FLOAT, 0, MPI_COMM_WORLD);
 
   b.north = north;
   b.south = south;
diff --git a/src/utils.c b/src/utils.c
index 7140bf8..f2115a7 100644
--- a/src/utils.c
+++ b/src/utils.c
@@ -2,18 +2,18 @@
 #include <stdlib.h>
 #include "utils.h"
 
-double *create_sa_matrix(int rows, int cols) {
-  double *x;
+float *create_sa_matrix(int rows, int cols) {
+  float *x;
 
-  x = (double *) malloc(rows * cols * sizeof(double));
+  x = (float *) malloc(rows * cols * sizeof(float));
   return x;
 }
 
-void destroy_sa_matrix(double *x) {
+void destroy_sa_matrix(float *x) {
   free(x);
 }
 
-void print_sa_matrix(double *x, int rows, int cols) {
+void print_sa_matrix(float *x, int rows, int cols) {
   int i, j;
   for (i = 0; i < rows; i++) {
     for (j = 0; j < cols; j++) {
@@ -24,18 +24,18 @@ void print_sa_matrix(double *x, int rows, int cols) {
   fflush(stdout);
 }
 
-double **create_matrix(int rows, int cols) {
+float **create_matrix(int rows, int cols) {
   int i;
-  double **x;
+  float **x;
   
-  x = (double **) malloc(rows * sizeof(double));
+  x = (float **) malloc(rows * sizeof(float));
   for (i = 0; i < rows; i++) {
-    x[i] = (double *) malloc(cols * sizeof(double));
+    x[i] = (float *) malloc(cols * sizeof(float));
   }
   return x;
 }
 
-void destroy_matrix(double **x, int rows) {
+void destroy_matrix(float **x, int rows) {
   int i;
 
   for (i = 0; i < rows; i++) {
@@ -44,7 +44,7 @@ void destroy_matrix(double **x, int rows) {
   free(x);
 }
 
-void print_matrix(double **x, int rows, int cols) {
+void print_matrix(float **x, int rows, int cols) {
   int i, j;
   for (i = 0; i < rows; i++) {
     for (j = 0; j < cols; j++) {
diff --git a/src/utils.h b/src/utils.h
index eb0b37c..ac3cbdb 100644
--- a/src/utils.h
+++ b/src/utils.h
@@ -6,7 +6,7 @@
 #ifdef ENABLE_LOG
 #  define LOG(x) x
 #else
-#  define LOG(x) (void) 0 
+#  define LOG(x) (void) 0
 #endif
 
 /*
@@ -17,30 +17,30 @@
 #define IDX(cols, r, c) ((r) * (cols) + (c))
 
 typedef struct borders {
-  double north;
-  double east;
-  double south;
-  double west;
+  float north;
+  float east;
+  float south;
+  float west;
 } borders;
 
 
 /*
  * Create a matrix stored in a single array.
  */
-double *create_sa_matrix(int rows, int cols);
+float *create_sa_matrix(int rows, int cols);
 
 /*
  * Destroy a single array matrix.
  */
-void destroy_sa_matrix(double *x);
+void destroy_sa_matrix(float *x);
 
 int sa_index(int cols, int r, int c);
 
 /*
  * Print a single array matrix.
  */
-void print_sa_matrix(double *x, int rows, int cols);
+void print_sa_matrix(float *x, int rows, int cols);
 
-double **create_matrix(int rows, int cols);
-void destroy_matrix(double **x, int rows);
-void print_matrix(double **x, int rows, int cols);
+float **create_matrix(int rows, int cols);
+void destroy_matrix(float **x, int rows);
+void print_matrix(float **x, int rows, int cols);