From af056b4d82d17da49b5a44d8060d097bd57ecc99 Mon Sep 17 00:00:00 2001 From: Hannes Signer Date: Mon, 21 Aug 2023 11:05:37 +0200 Subject: [PATCH] delete openmp curly bracket --- src/BTCSv2.cpp | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/src/BTCSv2.cpp b/src/BTCSv2.cpp index 775520d..50db8db 100644 --- a/src/BTCSv2.cpp +++ b/src/BTCSv2.cpp @@ -333,7 +333,6 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double ×tep) { MatrixXd concentrations = grid.getConcentrations(); #pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1) - { for (int i = 0; i < rowMax; i++) { @@ -347,14 +346,14 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double ×tep) { concentrations_t1.row(i) = row_t1; } - } + concentrations_t1.transposeInPlace(); concentrations.transposeInPlace(); alphaX.transposeInPlace(); alphaY.transposeInPlace(); #pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1) - { + for (int i = 0; i < colMax; i++) { // swap alphas, boundary conditions and sx/sy for column-wise calculation A = createCoeffMatrix(alphaY, bcLeft, bcRight, rowMax, i, sy); @@ -365,7 +364,7 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double ×tep) { concentrations.row(i) = row_t1; } - } + concentrations.transposeInPlace(); grid.setConcentrations(concentrations);