diff options
| -rw-r--r-- | src/refinement/coupler2D.hh | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/refinement/coupler2D.hh b/src/refinement/coupler2D.hh index 373a6de..c2b27db 100644 --- a/src/refinement/coupler2D.hh +++ b/src/refinement/coupler2D.hh @@ -109,6 +109,7 @@ void FineCoupler2D<T,DESCRIPTOR>::store() { auto& coarseLattice = this->_coarse.getSuperLattice(); + #pragma omp parallel for for (int y=0; y < this->_coarseSize; ++y) { const auto pos = this->getCoarseLatticeR(y); T rho{}; @@ -160,6 +161,7 @@ void FineCoupler2D<T,DESCRIPTOR>::interpolate() { auto& coarseLattice = this->_coarse.getSuperLattice(); + #pragma omp parallel for for (int y=0; y < this->_coarseSize; ++y) { Cell<T,DESCRIPTOR> coarseCell; coarseLattice.get(this->getCoarseLatticeR(y), coarseCell); @@ -184,6 +186,7 @@ void FineCoupler2D<T,DESCRIPTOR>::couple() const auto& coarseLattice = this->_coarse.getSuperLattice(); auto& fineLattice = this->_fine.getSuperLattice(); + #pragma omp parallel for for (int y=0; y < this->_coarseSize; ++y) { const auto& coarsePos = this->getCoarseLatticeR(y); const auto& finePos = this->getFineLatticeR(2*y); @@ -201,6 +204,7 @@ void FineCoupler2D<T,DESCRIPTOR>::couple() fineLattice.set(finePos, cell); } + #pragma omp parallel for for (int y=1; y < this->_coarseSize-2; ++y) { const auto rho = order4interpolation(_c2f_rho, y); const auto u = order4interpolation(_c2f_u, y); @@ -305,6 +309,7 @@ void CoarseCoupler2D<T,DESCRIPTOR>::couple() const auto& fineLattice = this->_fine.getSuperLattice(); auto& coarseLattice = this->_coarse.getSuperLattice(); + #pragma omp parallel for for (int y=0; y < this->_coarseSize; ++y) { const auto& finePos = this->getFineLatticeR(2*y); const auto& coarsePos = this->getCoarseLatticeR(y); |