SyFi  0.3
time_element_test.py
Go to the documentation of this file.
00001 
00002 from swiginac import *
00003 from SyFi import *
00004 
00005 t = ReferenceLine()
00006 triangle = ReferenceTriangle()
00007 fe = Lagrange(triangle, 2) 
00008 space_time_domain = SpaceTimeDomain(t, triangle)
00009 
00010 space_time_fe = SpaceTimeElement(t, 3, fe) 
00011 for i in range(0, space_time_fe.nbf()): 
00012     print "fe.N(%d)= %s " % (i, space_time_fe.N(i)) 
00013     print "fe.dof(%d)= %s " %  (i, space_time_fe.dof(i)) 
00014 
00015 
00016 for i in range(0, space_time_fe.nbf()): 
00017     for j in range(0, space_time_fe.nbf()): 
00018         integrand = space_time_fe.N(i)*space_time_fe.N(j)
00019         Aij = space_time_domain.integrate(integrand)
00020         print "A[%d, %d] = %s " %  (i,j, Aij)
00021 
00022         
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Defines