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