Dear all,
I am trying to solve an elasticity problem on a domain() with two different materials. I have the solution. But I have problem trying to get the solution just on the interface of the two subdomains and plotting the solution.
Take the official example "demo_subdomains-poisson.py" for instance, how should I modify it so that I can achieve the above goal? (In the example the unknown is a scalar function instead of a vector-valued function.) The code of the example can be found at:
http://fenicsproject.org/documentation/dolfin/1.6.0/python/demo/documented/subdomains-poisson/python/documentation.html
Thank you very much for your attention.
