Hello Yuliang,
If cell is a DoFHandler::active_cell_iterator, will
cell->at_boundary(0) and cell->vertex(0) always mean the left vertex
of the cell? I didn't find any information from the documentation.
It depends on the orientation of the cell, but if you did not create the
triangulation by hand, the answer is yes.
I'm programing on a 1D complex valued helmholtz equation with mixed
boundary conditions and I can't use quadrature points or FEFaceValues
to evaluate the boundary terms. Instead, I use cell->at_boundary(0) to
identify the left boundary and FiniteElement->shape_value(i,
cell->vertex(0)) to compute the shape value at the bounary. The
program runs but the solution is not correct compared with the exact
solution. I don't know whether my mathematical formulation or the code
is wrong.
Since, you did not give us the mathematical formulation, we cannot say
anything about that point.
Best Regards,
Markus
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
