ft: hybrid implementation for mg solve v0

2025-03-30 22:15:38 +02:00
parent 3408d39545
commit e714bdd2c8
2 changed files with 77 additions and 9 deletions
--- a/EnhancedSolver/2D-mpi/src/discretization.c
+++ b/EnhancedSolver/2D-mpi/src/discretization.c
@@ -114,7 +114,9 @@ void initDiscretiztion(Discretization* d, Parameter* params)
    d->f      = allocate(64, size * sizeof(double));
    d->g      = allocate(64, size * sizeof(double));
    d->grid.s = allocate(64, size * sizeof(double));
-
+    #ifdef _OPENMP
+    #pragma omp parallel for  
+    #endif
    for (int i = 0; i < size; i++) {
        d->u[i]      = params->u_init;
        d->v[i]      = params->v_init;