summaryrefslogtreecommitdiff
path: root/src/refinement
diff options
context:
space:
mode:
Diffstat (limited to 'src/refinement')
-rw-r--r--src/refinement/coupler2D.hh5
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);