The problem I am currently having is directly related to my previous post [here](https://groups.google.com/forum/#!searchin/dealii/krei/dealii/eq_zP0jrSJU/SyQXA7A-DAAJ).
The physics are exactly the same as described in that thread: I have a two domain system: vacuum and metal, both domains have their own meshes. [Here is a picture](http://i.imgur.com/CejXi1y.png). The meshes are exactly matching on the boundary. I want to solve the Laplace equation for electric fields in the vacuum part (does not depend on anything in the metal part) and then use the electric field as a boundary condition to the metal part. I am currently using VectorTools::point_gradient to evaluate the electric field from the Laplace problem in the quadrature points of the boundary cell faces of the metal mesh. As a result, it takes over an hour to assemble the system (solving takes less than a second). Now, considering that the mesh is exactly matching at the boundary (all the quadrature points etc), how could I efficiently evaluate the electric field at the boundary? And I apologize if I'm mistaken, but there doesn't seem to be any tutorials on these kinds of problems where you have multiple domains with multiple meshes. Or is it usually done by using a single mesh and somehow defining different domains in that mesh? (haven't delved deeply into Step-46, but there they do something like this, right?) -- 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.