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