Nice catch! I forgot to mention it is bilinear. But thanks for all the
explanation.
在 2018年3月28日星期三 UTC-7下午12:42:27，Wolfgang Bangerth写道：
>
>
> > |
> > // The vector I want to output is the solution for vector-valued
> problem.
> > // I want to output the cell-average value for the 0th componen
Hi Weixiong,
The issue that you’re seeing here is actually related to a discussion that
we’re currently having regarding the behaviour of boundaries and manifolds (and
will hopefully be fixed soon in the developer version). You’re hitting the
issue discussed here
https://github.com/dealii/deali
I changed the strategy that I don't assign manifold first. I merge two
triangulations. Then, I loop over cells to see if cell has a side on the
merged interior boundary. If so, I set manifold id to it.
const Point<2> center (0,0);
const SphericalManifold<2> boundary_description(center);
tri
I think I'd better describe a little bit more.
I want a square domain with a cylinder inside it. What I tried doing was to
get a square with a cylindrical hole merged with a cylinder.
Here are how I did it:
1. generate the cylinder using GridGenerator::general_cell and
merge_triangulation to
Hi,
I just started to try manifold and tria merging functionalities. My
question comes from the problem that I have to merge two curved
triangulations. When refining, how do I make sure the merged boundaries are
still curved?
I currently set manifold before merging triangulation but it seems n
Hi Prof. Bangerth,
I would like to post that snippet of code. To be honest, I verified the
solution in a naive way, that is, I compare the two different number of
processes solutions in paraview. The comparison(including data range, plot
curve over time at a point) shows they are exactly the same.
|
// The vector I want to output is the solution for vector-valued problem.
// I want to output the cell-average value for the 0th component.
Vectormodi (solu);
Quadratureqq(fe.get_unit_support_points());
FEValuesfv(fe,qq,update_values |update_quadrature_points);
std::vectorselected_dofs (dof_ha
Qing,
I meet a MPI problem when I simulate a three point bending test. I try
to compute the load force traction on the top boundary. The results are
correct
with the single process. But When I switched to a multi-process, the
results changed.
Let's say, if we wan to get the load force in y di
