@@ -61,7 +61,6 @@ class gmd_2012 : public libmpdataxx::solvers::mpdata<ct_params_t>
6161 }
6262 }
6363
64- // this->mem->barrier();
6564 this ->xchng_sclr (gc_node[0 ], this ->i , this ->j );
6665 this ->xchng_sclr (gc_node[1 ], this ->i , this ->j );
6766
@@ -70,34 +69,13 @@ class gmd_2012 : public libmpdataxx::solvers::mpdata<ct_params_t>
7069 {
7170 for (int j = this ->j .first ()-1 ; j <= this ->j .last (); ++j)
7271 {
73- // const auto x = i * this->di;
74- // const auto y = (j+ 0.5) * this->dj - pi / 2;
75- // auto gc_ij_x = gc_exact(t, x + 0.5 * this->di, y);
76- // auto gc_ij_y = gc_exact(t, x, y + 0.5 * this->dj);
77- // this->mem->GC[0](i+h, j) = gc_ij_x.first;
78- // this->mem->GC[1](i, j+h) = gc_ij_y.second;
79-
8072 this ->mem ->GC [0 ](i+h, j) = 0.5 * (gc_node[0 ](i, j) + gc_node[0 ](i+1 , j));
8173 this ->mem ->GC [1 ](i, j+h) = 0.5 * (gc_node[1 ](i, j) + gc_node[1 ](i, j+1 ));
8274 }
8375 }
8476
8577 auto ex = this ->halo - 1 ;
86- // this->xchng_vctr_alng(this->mem->GC);
8778 this ->xchng_vctr_nrml (this ->mem ->GC , this ->i ^ex, this ->j ^ex);
88-
89- // auto div = max(abs(
90- // ( this->mem->GC[0](this->i+h, this->j) - this->mem->GC[0](this->i-h, this->j)
91- // + this->mem->GC[1](this->i, this->j+h) - this->mem->GC[1](this->i, this->j-h)
92- // ) / (*this->mem->G)(this->i, this->j)
93- // ));
94-
95- // if (this->time == 0)
96- // {
97- // std::stringstream ss;
98- // ss << "rank: " << this->rank << " div: " << div << std::endl;
99- // std::cout << ss.str() << std::endl;
100- // }
10179
10280 return true ;
10381 }
0 commit comments