max/tug
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/tug.git synced 2025-12-15 18:38:23 +01:00
Code Issues Packages Projects Releases Wiki Activity

delete openmp curly bracket

Browse Source
This commit is contained in:
Hannes Signer 2023-08-21 11:05:37 +02:00
parent c0e46bd82f
commit af056b4d82
1 changed files with 3 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/BTCSv2.cpp
View File

@ -333,7 +333,6 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double &timestep) {
MatrixXd concentrations = grid.getConcentrations(); MatrixXd concentrations = grid.getConcentrations();
#pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1) #pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1)
{
for (int i = 0; i < rowMax; i++) { for (int i = 0; i < rowMax; i++) {
@ -347,14 +346,14 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double &timestep) {
concentrations_t1.row(i) = row_t1; concentrations_t1.row(i) = row_t1;
} }
}
concentrations_t1.transposeInPlace(); concentrations_t1.transposeInPlace();
concentrations.transposeInPlace(); concentrations.transposeInPlace();
alphaX.transposeInPlace(); alphaX.transposeInPlace();
alphaY.transposeInPlace(); alphaY.transposeInPlace();
#pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1) #pragma omp parallel for num_threads(NUM_THREADS_BTCS) private(A, b, row_t1)
{
for (int i = 0; i < colMax; i++) { for (int i = 0; i < colMax; i++) {
// swap alphas, boundary conditions and sx/sy for column-wise calculation // swap alphas, boundary conditions and sx/sy for column-wise calculation
A = createCoeffMatrix(alphaY, bcLeft, bcRight, rowMax, i, sy); A = createCoeffMatrix(alphaY, bcLeft, bcRight, rowMax, i, sy);
@ -365,7 +364,7 @@ static void BTCS_2D(Grid &grid, Boundary &bc, double &timestep) {
concentrations.row(i) = row_t1; concentrations.row(i) = row_t1;
} }
}
concentrations.transposeInPlace(); concentrations.transposeInPlace();
grid.setConcentrations(concentrations); grid.setConcentrations(concentrations);