aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorJaroCamphuijsen <jjecamphuijsen@gmail.com>2016-06-08 12:25:59 +0200
committerJaroCamphuijsen <jjecamphuijsen@gmail.com>2016-06-08 12:25:59 +0200
commit208a434c1ae6ab2191e516d73bc8724f686a2624 (patch)
tree12b75397afcc54783cf6fcd3940337f25d9c3302
parent9eba1b7f8bccdb1898bce23a242eeba5c1e7d9e0 (diff)
Fixed circular domain
-rw-r--r--Set 3/eigenmodes.py8
1 files changed, 4 insertions, 4 deletions
diff --git a/Set 3/eigenmodes.py b/Set 3/eigenmodes.py
index 48fa613..31eee5c 100644
--- a/Set 3/eigenmodes.py
+++ b/Set 3/eigenmodes.py
@@ -16,7 +16,7 @@ def circDomain(size, M):
for j in range(size):
if sqrt((i-size/2.)**2 + (j-size/2.)**2) > size/2:
M[(i*size + j)] = 0
- # M[:(i*size + j)] = 0
+ M[:,(i*size + j)] = 0
# M = np.delete(M,(i*size + j),0)
return M
@@ -25,7 +25,7 @@ def rectDomain(size, M, frac):
for j in range(size):
if i > size*frac:
M[(i*size + j)] = 0
- # M[:(i*size + j)] = 0
+ M[:,(i*size + j)] = 0
# M = np.delete(M,(i*size + j),0)
return M
@@ -202,10 +202,10 @@ def plot_time():
if __name__=="__main__":
- # plot_eigenmodes()
+ plot_eigenmodes()
# plot_L_dependence()
# plot_dstep_dependence()
- plot_time()
+ # plot_time()