Hello Wolfgang,
thank you for the hint.
This is just what I was looking for. I think I can simply copy the
VectorTools::no_normal_flux_constraints, because this function actually
only deals with the difficulties on finding the correct normal vectors.
What I have to do is to implement an equivalent to the
internals::VectorTools::add_constraints function, where the creation
and addition of the constraints takes place, but this looks straight
forward.
Best Regards,
Markus
Wolfgang Bangerth schrieb:
n \times E_s = n \times E_i
[...]
Did anyone already implement something like this or does anyone know how
to do it?
I don't know of anyone who would already have done that. But you may want
to look at the VectorTools::compute_no_normal_flux_constraints function,
which has quite some documentation on how to compute the constraints for
the boundary condition
n cdot E = 0
I would imagine that you can compute constraints for the tangential
conditions in a similar way.
Let us know if you need help!
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
|
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii