Praveen, the easiest (and widely popular) way to implement weak boundary
conditions is Nitsche's method. While you refer to a Lagrange multiplier
method, this one is an augmented Lagrangian.
LocalIntegrators::Laplace::nitsche_matrix and
LocalIntegrators::Laplace::nitsche_residual can be used to
implement it.
Best,
Guido
On 07/07/2011 08:34 AM, Praveen C wrote:
Hello
Is it possible to enforce dirichlet BC in weak form, the so called
primal mixed approach. An example is given here, see Remark 7.1.4
http://books.google.co.in/books?id=j7WmQ1yFjkAC&lpg=PP1&dq=quarteroni%20numerical%20approximation&pg=PA227#v=onepage&q&f=false
<http://books.google.co.in/books?id=j7WmQ1yFjkAC&lpg=PP1&dq=quarteroni%20numerical%20approximation&pg=PA227#v=onepage&q&f=false>
This requires defining finite elements on the boundary. I dont
remember seeing any examples in deal.ii which do this.
praveen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii