SyFi
0.3
|
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