Hi all,
I have a mesh (generated with gmsh) that has an internal line inside of a
domain to evaluate an output along the line. As shown in the below mesh
figure, the line is composed of element edges from (-2,-0.5) to (-2.0,0.5).
I set the internal line as a physical line and try to get an output, but it
seems my output functional evaluation has never called. So I switched my
strategy to pass the coordinates of the element edges along the internal
line such as
virtual void boundary_assembly(FEMContext &c) {
AutoPtr<Elem> elem_side = c.elem->build_side(c.side);
const Real centroid_x = elem_side->centroid()(0);
const Real centroid_y = elem_side->centroid()(1);
if( (-2.0 == centroid_x) && ((-0.5 <= centroid_y) || (centroid_y <=
0.5)) ) {
.......
}
but again my output evaluation has never called --- it looks like there is
no element edges satisfying (-2.0 == centroid_x). I'm afraid probably the
way I specify coordinates may not correct. Plz suggest me some ideas to
fix the problem.
------------------------------------------------------------------------------
RSA(R) Conference 2012
Save $700 by Nov 18
Register now
http://p.sf.net/sfu/rsa-sfdev2dev1
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
