> 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
Yes, wrong approach :-) The way to do what you want to do is to say fe_face_values.quadrature_point(q) or fe_face_values.get_quadrature_points() to get at the coordinates of the quadrature points *in real space* (as opposed to asking the quadrature object which gives you the coordinates on the reference cell and then transforming them yourself). You need to add update_quadrature_points to the update flags of the FEFaceValues object constructor call, of course. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
