Dear users, I am trying to implement my finite element model in dealii. the weak form of my problem is <https://lh3.googleusercontent.com/--o9P1GVqOns/V3Y5WGBQY7I/AAAAAAAAGDY/G76lijaZJbUNY5EmCDMkhqKS-rxwWBQGQCLcB/s1600/CodeCogsEqn%2B%25281%2529.gif> where as <https://lh3.googleusercontent.com/-l5QxInGUqRc/V3Y5deBB5PI/AAAAAAAAGDc/CPEEiORknPQycVeR4OFBRoMcq3qyfUkCwCLcB/s1600/CodeCogsEqn%2B%25282%2529.gif> and Ti is some internal boundary in the domain. My idea was to implement it using the method use in the example 46 and followed as 1) i made two seperate domains in gmash and marked the elements on the seleted edge with material_id =left and material_id=right. 2) I merged them to the single triangulation so the material id information is transfered to the new mesh with all outboundary id = 0 3) i want to loop over the cell with material id=left and check that if cell->face->neighbour-cell ->material_id=right than do the face integration on the face. but i did not found any function that can check neighbour-cell material id. Is there any easy way to do this problem or any alternate suggestion. and second thing I also wish to add some stiffness to some node. can i mark the node for this purpose to add some stiffness to stiffness matrix. thanks in advance asim khan -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.