SyFi
0.3
|
00001 00002 00003 from swiginac import * 00004 from SyFi import * 00005 00006 t = ReferenceTriangle() 00007 fe = Lagrange(t,5) 00008 00009 u = 0 00010 us = [] 00011 for i in range(0, fe.nbf()): 00012 ui = symbol("u_%d" % i) 00013 us.append(ui) 00014 u += ui*fe.N(i) 00015 00016 00017 J = matrix(fe.nbf(), fe.nbf()) 00018 for i in range(0, fe.nbf()): 00019 Fi = u*fe.N(i) 00020 for j in range(0, fe.nbf()): 00021 uj = us[j] 00022 integrands = diff(Fi, uj) 00023 J[j,i] = t.integrate(integrands) 00024 00025 00026 #print J 00027 00028 00029