コード例 #1
0
ファイル: test.py プロジェクト: tbetcke/PyPWDG
def basiccube():
    dir = na([[1,0,0]])
    k = 15
    maxvol = 0.02
    g = dg3d.impedance(dir,k)
    gg = dg3d.planeWaves(dir,k)[0]
    s = dg3d.solver(dg3d.cubemesh(maxvol), 2, 12, k).solve(g)
    v = dg3d.Visualiser(s,gg)
コード例 #2
0
ファイル: test.py プロジェクト: tbetcke/PyPWDG
    t.append(time.time())
    B2 = numpy.sum(A2, axis=0)
    t.append(time.time())
    l = len(A2)
    B3 = numpy.reshape(numpy.dot(numpy.ones(l), numpy.reshape(A2, (l, -1))), (n,n))
    t.append(time.time())
    op = lambda a,b:a+b
    B4 = reduce(op, A2)
    t.append(time.time())
     
    print na(t[1:]) - na(t[:-1])

if __name__ == '__main__':
    print "hello"
    dir = na([[1,0,0]])
    k = 20
    maxvol = 0.04
    print k, maxvol
    g = dg3d.impedance(dir,k)
    gg = dg3d.planeWaves(dir,k)[0]
    s = dg3d.solver(dg3d.cubemesh(maxvol), 4, 10, k)
    
    s.solve(g)
    
    v = dg3d.Visualiser(s,gg)
    v.showuaveragereal()
    
    sys.exit(0)