Main Page | Namespace List | Class Hierarchy | Class List | Directories | File List | Namespace Members | Class Members | File Members | Search for

# raviartthomas_ex2.cpp

Go to the documentation of this file.

    00001 #include <SyFi.h>
    00002
    00003 using namespace GiNaC;
    00004
    00005 int main () {
    00006   int order = 3;
    00007
    00008   ReferenceTriangle triangle("t");
    00009   RaviartThomas vfe;
    00010   vfe.set(triangle);
    00011   vfe.set(order);
    00012   vfe.compute_basis_functions();
    00013
    00014   DiscontinuousLagrangeFE pfe;
    00015   pfe.set(triangle);
    00016   pfe.set(order);
    00017   pfe.compute_basis_functions();
    00018
    00019   for (int i=0; i< vfe.nbf(); i++)
    00020     cout <<"vfe.N("<<i<<")="<<vfe.N(i)<<endl;
    00021
    00022   for (int i=0; i< pfe.nbf(); i++)
    00023     cout <<"pfe.N("<<i<<")="<<pfe.N(i)<<endl;
    00024
    00025   return 0;
    00026 }
    00027
    00028
    00029
    00030
    00031
    00032
    00033