Hello,
I have following problem: I have transport equation in 2D, which has 1D
boundary part, on which Neumann boundary condition is put.
I compute the boundary integral in the following way
for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
if (cell->face(face)->at_boundary()
&&
(((cell->face(face)->boundary_indicator()) == 3)))
{
.....
I am computing over quadrature points
for(unsigned int q=0; q<n_face_q_points; ++q)
where n_face_q_points and quadrature on the face is defined:
const unsigned int n_face_q_points = face_quadrature_formula.size();
QGauss<dim-1> face_quadrature_formula(degree+2);
const Quadrature<dim-1> *quadrature_face = &fe_face_values.get_quadrature();
But now I need for computations to know what are real coordinates of each
quadrature point (in 2D)
I have solved it by following way, but it works bad (bad results of
computations, it could be compiled and runned)
Point<dim> point = *new
Point<dim>(1.5,quadrature_face->point(q)(0)); //I know, that x
coordinate of boundary should be 1.5
globalvector=mapping->transform_unit_to_real_cell(cell,
quadrature_face->point(q));
ikska[pocitadlo_q]=globalvector(1);
//desired coordinate y of quadrature point q
Global 2D mapping is defined
const Mapping<dim> *mapping = &fe_values.get_mapping();
Please help me!
Marek Capek
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
