fe_ex2.cpp

Go to the documentation of this file.
00001 #include <SyFi.h>
00002 
00003 int main() {
00004   Triangle T(lst(0,0), lst(1,0), lst(0,1), "t"); 
00005   int order = 2; 
00006 
00007   std::map<std::pair<int,int>, ex> A;
00008   std::pair<int,int> index;
00009   LagrangeFE fe; 
00010   fe.set(order); 
00011   fe.set(T); 
00012   fe.compute_basis_functions(); 
00013   for (int i=1; i<= fe.nbf() ; i++) {
00014     index.first = i; 
00015     for (int j=1; j<= fe.nbf() ; j++) {
00016       index.second = j; 
00017       ex nabla = inner(grad(fe.N(i)), grad(fe.N(j))); 
00018       ex Aij = T.integrate(nabla);   
00019       A[index] = Aij;  
00020     }
00021   }
00022 
00023   map<std::pair<int,int>,ex>::iterator iter; 
00024   for (iter = A.begin(); iter != A.end() ; iter++) {
00025       ex intf = (*iter).second; 
00026       cout <<"A["<<(*iter).first.first<<","<<(*iter).first.second<<"]="<<(*iter).second<<endl; 
00027   }
00028 
00029 
00030 
00031 }
00032 
00033 
00034 

Generated on Tue Nov 22 11:10:22 2005 for SyFi by  doxygen 1.4.4