本文整理汇总了C++中MultiGrid::poisson_zone_is_refined方法的典型用法代码示例。如果您正苦于以下问题:C++ MultiGrid::poisson_zone_is_refined方法的具体用法?C++ MultiGrid::poisson_zone_is_refined怎么用?C++ MultiGrid::poisson_zone_is_refined使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类MultiGrid
的用法示例。
在下文中一共展示了MultiGrid::poisson_zone_is_refined方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: get_phi_at
Real MultiGrid::get_phi_at(Real x, Real y, Real z) {
Real p, tmp;
MultiGrid* g;
p = 0.0;
for (int l = 0; l < get_local_node_cnt(); l++) {
g = dynamic_cast<MultiGrid*>(get_local_node(l));
if (MPI_rank() == g->proc()) {
for (int k = 1; k < PNX - 1; k++) {
if (z >= g->MultiGrid::zf(k) && z < g->MultiGrid::zf(k + 1)) {
for (int j = 1; j < PNX - 1; j++) {
if (y >= g->MultiGrid::yf(j) && y < g->MultiGrid::yf(j + 1)) {
for (int i = 1; i < PNX - 1; i++) {
// printf("%e %e %e\n", g->MultiGrid::xf(i), x, g->MultiGrid::xf(i + 1));
if (x >= g->MultiGrid::xf(i) && x < g->MultiGrid::xf(i + 1)) {
if (!g->poisson_zone_is_refined(i, j, k)) {
p = g->phi(i, j, k);
// printf("%e\n", p);
}
}
}
}
}
}
}
}
}
tmp = p;
MPI_Allreduce(&tmp, &p, 1, MPI_DOUBLE_PRECISION, MPI_SUM, MPI_COMM_WORLD );
return p;
}