script_base.c

#include <stdio.h>
#include <stdlib.h>

// Usado para monitorear tiempo de ejecucion
#include <time.h>


// Definir el tamaño de la cuadrícula y el número de pasos de simulación
#define SIZE 3000    // Tamaño de la cuadrícula (3000x3000)
#define STEPS 10000  // Número de iteraciones para la simulación

int main() {
    // Variables para medir tiempo
    clock_t start, end;
    double cpu_time_used;

    // Reservar memoria dinámica para la cuadrícula como un array unidimensional
    double *grid = malloc(SIZE * SIZE * sizeof(double));
    if (grid == NULL) {
        // Verificar si la asignación de memoria falló
        fprintf(stderr, "Error: Unable to allocate memory for grid.\n");
        return EXIT_FAILURE;
    }

    start = clock(); // Comenzar medicion del tiempo

    int i, j, step; // Variables de control para los bucles

    // Inicializar toda la cuadrícula con valores de 0.0
    for (i = 0; i < SIZE * SIZE; i++) {
        grid[i] = 0.0;
    }

    // Establecer un valor inicial alto en el centro de la cuadrícula
    grid[(SIZE / 2) * SIZE + (SIZE / 2)] = 10000.0;

    // Inicio de la simulación de difusión
    for (step = 0; step < STEPS; step++) {
        // Actualizar los valores de cada celda (excepto los bordes)
        for (i = 1; i < SIZE - 1; i++) {
            for (j = 1; j < SIZE - 1; j++) {
                // Actualizar la celda actual como el promedio de sus vecinos
                grid[i * SIZE + j] = 0.25 * (
                    grid[(i - 1) * SIZE + j] +   // Celda de arriba
                    grid[(i + 1) * SIZE + j] +   // Celda de abajo
                    grid[i * SIZE + (j - 1)] +   // Celda de la izquierda
                    grid[i * SIZE + (j + 1)]     // Celda de la derecha
                );
            }
        }
    }

    end = clock(); // Terminar medicion del tiempo
    cpu_time_used = ((double) (end - start)) / CLOCKS_PER_SEC;

    // Imprimir el valor final en el centro de la cuadrícula
    printf("Final grid center value: %f\n", grid[(SIZE / 2) * SIZE + (SIZE / 2)]);
    printf("Time: %f", cpu_time_used);

    // Liberar la memoria asignada para la cuadrícula
    free(grid);

    return 0; // Terminar el programa con éxito
}