SyFi 0.3
syfi_test.py
Go to the documentation of this file.
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 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Defines