The driven cavity benchmark: 2D boundary conditions
space Xh (omega, approx,
"vector");
Xh.block("top"); Xh.block("bottom");
if (omega.dimension() == 3) {
Xh.block("back"); Xh.block("front");
Xh[1].block("left"); Xh[1].block("right");
} else {
Xh.block("left"); Xh.block("right");
}
return Xh;
}
return uh;
}
string valued = (omega.dimension() == 3) ? "vector" : "scalar";
space Ph (omega, approx, valued);
Ph.block("top"); Ph.block("bottom");
if (omega.dimension() == 3) {
Ph.block("back"); Ph.block("front");
} else {
Ph.block("left"); Ph.block("right");
}
return Ph;
}
}
};