Hello Wolfgang,
I already implemented this function some time ago, but planned to send
it all together as one package, when the higher order Nédélec elements
are ready to release. I will send the curl function to you today such
that it can be included.
Best Regards,
Markus
Am 13.05.10 04:49, schrieb Wolfgang Bangerth:
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
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
