Hello Arun,

it is in the fe_values.cc/.h files. But you do not have to search for it. You have just to download deal.II using svn and recompile it. Then it should be available in your copy of deal.II, too.

Best Regards,
Markus



Am 13.05.10 18:17, schrieb arun jaganathan:
Hello Markus,

I am searching for the curl function that you implemented within svn link and not sure what to look out for. Could you tell me the name/location of the function within this link ?

http://www.dealii.org/svn/dealii/trunk/deal.II/

Too bad I cannot even locate it :(

Arun.

On Thu, May 13, 2010 at 1:28 AM, Markus Bürg <[email protected] <mailto:[email protected]>> wrote:

    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]
        <mailto:[email protected]>
                                         www:
        http://www.math.tamu.edu/~bangerth/
        <http://www.math.tamu.edu/%7Ebangerth/>

        _______________________________________________
        dealii mailing list
        http://poisson.dealii.org/mailman/listinfo/dealii


    _______________________________________________
    dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to