42 |
t=0 |
t=0 |
43 |
|
|
44 |
# define the location of the point source |
# define the location of the point source |
45 |
L=Locator(domain,xc) |
L=Locator(domain,numarray.array(xc)) |
46 |
# find potential at point source |
# find potential at point source |
47 |
u_pc=L.getValue(u) |
u_pc=L.getValue(u) |
48 |
print "u at point charge=",u_pc |
print "u at point charge=",u_pc |
70 |
t+=h |
t+=h |
71 |
n+=1 |
n+=1 |
72 |
print n,"-th time step t ",t |
print n,"-th time step t ",t |
|
L=Locator(domain,xc) |
|
73 |
u_pc=L.getValue(u) |
u_pc=L.getValue(u) |
74 |
print "u at point charge=",u_pc |
print "u at point charge=",u_pc |
75 |
|
|