Hi, I think we're talking about different things here. I'm totally aware of the boundary conditions of my system and found examples how to implement these constraints using deal.II. What I meant was looking for example at the "vector-valued problems" module description, one has the final form (v,u) - (div v,p) - (q, div u) = (q, f) + (nv, p)_\Omega.
In the implementation however the last term in the above form, the boundary term, is neglected. I'm now wondering how would this look like if I don't neglect it. Somehow I then have to be able to evaluate the shape functions and their gradients on the surface of each cell. How can this be done? Best regards, Klaus ___________________________________________________________ Schon gehört? WEB.DE hat einen genialen Phishing-Filter in die Toolbar eingebaut! http://produkte.web.de/go/toolbar _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
