00001 #include <SyFi.h> 00002 00003 using namespace GiNaC; 00004 00005 int main() { 00006 00007 ReferenceTriangle domain; 00008 00009 VectorCrouzeixRaviart v_fe; 00010 v_fe.set_size(2); 00011 v_fe.set(domain); 00012 v_fe.compute_basis_functions(); 00013 00014 P0 p_fe; 00015 p_fe.set(domain); 00016 p_fe.compute_basis_functions(); 00017 00018 usage(v_fe, p_fe); 00019 00020 00021 Dof dof; 00022 std::map<std::pair<int,int>, ex> A; 00023 compute_Stokes_element_matrix(v_fe, p_fe, dof, A); 00024 print(A); 00025 00026 } 00027