From: Lorenzo Botti <bottilorenzo@gm...>  20090529 10:57:54

Hi all,I'm trying to avoid inverse mapping for side integration with the following procedure. 1) search for the elem local nodes indices on the side (they are stored in side_nodes_map). 2) use a method that returns reference space nodes coordinates for the element. 3) map the side quadrature points on the element's side using reference space coordinates. Nodes ordering match map ordering. This works for boundary side integration. For internal side integration (DG) 4) search for the neighbor local nodes indices on the element's side (the local indices ordering will be different from side_nodes_map). This could be done with: loop over side nodes (side is the side of the element) loop over neighbor nodes if side node id is equal to neighbor node id store the neighbor loop index. Maybe you have better ideas... 5) use the same method of point 2) 6) map the side quadrature points on the neighbor's side using reference space coordinates. Nodes ordering match map ordering. This works for conforming meshes and uses the side mapping we have already computed. For non conforming meshes I would like to use embedded matrix to compute neighbor reference space side coordinates. Is it possible to get the side child number from the element child number? I know that sides have no parents but it would be useful. Thanks Lorenzo 