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
00027
00028
00029