> After going through the tutorials and reports I am slowly getting
> familiar with deal. Very good library to learn for sure because of its
> flexibility. I am wondering how to implement maxwell's equation whose
> weak form will have (curl U, curl V). How can I implement this ? I found
> a report on Nedelec elements by Anna which talks about this issue but I
> haven't completely followed it yet in terms of its implementation in
> deal. Can someone provide some hints ?
You would simply have to choose the FE_Nedelec class. Assembling the
bilinear form then proceeds as as discussed in the documentation module on
vector-valued problems.
It might be worthwhile to implement a curl function so that you can do
FEValuesExtractors::Vector E(0);
local_matrix(i,j)
+= fe_values[E].curl(i,q) * ...
Right now you would have to compute the curl from the components of the
gradient matrix but it would be simpler if a function
FEValuesViews::Vector::curl
existed. Let me know if you want help implementing such a function.
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
