mxpoisson_ex.cpp File Reference

#include <SyFi.h>

Go to the source code of this file.

Functions

int main ()


Function Documentation

int main  ) 
 

Definition at line 5 of file mxpoisson_ex.cpp.

References DiscontinuousLagrangeFE::compute_basis_functions(), P0::compute_basis_functions(), RaviartThomas::compute_basis_functions(), compute_mixed_Poisson_element_matrix(), print(), DiscontinuousLagrangeFE::set(), StandardFE::set(), RaviartThomas::set(), and usage().

00005            {
00006 
00007   ReferenceTriangle domain; 
00008 
00009   // First the lowest order. 
00010   // The lowest order is special, because we must use P0 
00011   RaviartThomas v_fe;
00012   v_fe.set(1);
00013   v_fe.set(domain);
00014   v_fe.compute_basis_functions();
00015 
00016   P0 p_fe; 
00017   p_fe.set(domain);
00018   p_fe.compute_basis_functions();
00019 
00020   usage(v_fe, p_fe); 
00021 
00022   Dof dof; 
00023   std::map<std::pair<int,int>, ex> A; 
00024   compute_mixed_Poisson_element_matrix(v_fe, p_fe, dof, A);  
00025   print(A); 
00026 
00027 
00028   // Then for instance order 3
00029   int order = 3; 
00030   RaviartThomas v_fe2;
00031   v_fe2.set(order);
00032   v_fe2.set(domain);
00033   v_fe2.compute_basis_functions();
00034 
00035   DiscontinuousLagrangeFE p_fe2; 
00036   p_fe2.set(order); 
00037   p_fe2.set(domain);
00038   p_fe2.compute_basis_functions();
00039 
00040   usage(v_fe2, p_fe2); 
00041 
00042   Dof dof2; 
00043   std::map<std::pair<int,int>, ex> A2; 
00044   compute_mixed_Poisson_element_matrix(v_fe2, p_fe2, dof2, A2);  
00045   print(A2); 
00046 
00047 
00048 
00049 }


Generated on Mon Jan 9 18:08:08 2006 for SyFi by  doxygen 1.4.4