Hi Wolfgang & Markus, It will be really helpful if one of one of you show me how to implement this curl function. Sorry to ask without trying first. I am not there yet.
Arun. On Wed, May 12, 2010 at 9:49 PM, Wolfgang Bangerth <[email protected]>wrote: > > > 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
