質問編集履歴
5
```にしました
title
CHANGED
File without changes
|
body
CHANGED
@@ -8,7 +8,7 @@
|
|
8
8
|

|
9
9
|
|
10
10
|
### 該当のソースコード
|
11
|
-
|
11
|
+
```python
|
12
12
|
Ex = np.zeros([Lx, Ly])
|
13
13
|
Ey = np.zeros([Lx, Ly])
|
14
14
|
|
@@ -26,13 +26,13 @@
|
|
26
26
|
|
27
27
|
plt.xlim([0, Lx]) # 描くXの範囲
|
28
28
|
plt.ylim([0, Ly]) # 描くYの範囲
|
29
|
-
|
29
|
+
```
|
30
30
|
### 試したこと
|
31
31
|
|
32
32
|
delta_Lxとdelta_Lyを逆にした.
|
33
33
|
|
34
34
|
### 補足情報(FW/ツールのバージョンなど)
|
35
|
-
|
35
|
+
```python
|
36
36
|
delta_Lx = 0.01 # グリッドx幅
|
37
37
|
delta_Ly = 0.01 # グリッドy幅
|
38
38
|
LLx = 1
|
@@ -47,10 +47,10 @@
|
|
47
47
|
phi = np.zeros([Lx+1, Ly+1])
|
48
48
|
phi_Bound = np.zeros([Lx+1, Ly+1])
|
49
49
|
phi_Bound[0, :] = V # 全ての行の0列目を取得
|
50
|
-
|
50
|
+
```
|
51
51
|
|
52
52
|
# 電荷密度の設定
|
53
|
-
|
53
|
+
```python
|
54
54
|
eps0 = 1
|
55
55
|
charge = np.zeros([Lx+1, Ly+1])
|
56
56
|
c = 1500
|
@@ -67,11 +67,9 @@
|
|
67
67
|
aa_recta = 0.5 * (np.cos(np.pi / Lx) + np.cos(np.pi / Ly))
|
68
68
|
omega_SOR_recta = 2 / (1 + np.sqrt(1 - aa_recta ** 2)) # SOR法における加速パラメータ
|
69
69
|
print("omega_SOR_rect=", omega_SOR_recta)
|
70
|
-
|
70
|
+
```
|
71
|
-
|
72
|
-
|
73
71
|
# メイン
|
74
|
-
|
72
|
+
```python
|
75
73
|
delta = 1.0
|
76
74
|
n_iter = 0
|
77
75
|
conv_check = []
|
@@ -97,8 +95,7 @@
|
|
97
95
|
|
98
96
|
im = plt.imshow(phi, cmap = 'jet')
|
99
97
|
anim.append([im])
|
100
|
-
|
98
|
+
```
|
101
|
-
|
102
99
|
前進差分法により離散化をしたときの電場のグラフです.
|
103
100
|

|
104
101
|
このように中心差分法の時も描かれないといけませんが, 上記に貼った通り描かれていないのが現状です,
|
4
intendしました
title
CHANGED
File without changes
|
body
CHANGED
@@ -63,6 +63,7 @@
|
|
63
63
|
charge[i, j] = c
|
64
64
|
|
65
65
|
# for SOR method
|
66
|
+
’’’Python
|
66
67
|
aa_recta = 0.5 * (np.cos(np.pi / Lx) + np.cos(np.pi / Ly))
|
67
68
|
omega_SOR_recta = 2 / (1 + np.sqrt(1 - aa_recta ** 2)) # SOR法における加速パラメータ
|
68
69
|
print("omega_SOR_rect=", omega_SOR_recta)
|
3
’’’としました
title
CHANGED
File without changes
|
body
CHANGED
@@ -8,7 +8,7 @@
|
|
8
8
|

|
9
9
|
|
10
10
|
### 該当のソースコード
|
11
|
-
|
11
|
+
’’’Python
|
12
12
|
Ex = np.zeros([Lx, Ly])
|
13
13
|
Ey = np.zeros([Lx, Ly])
|
14
14
|
|
@@ -26,13 +26,13 @@
|
|
26
26
|
|
27
27
|
plt.xlim([0, Lx]) # 描くXの範囲
|
28
28
|
plt.ylim([0, Ly]) # 描くYの範囲
|
29
|
-
|
29
|
+
’’’
|
30
30
|
### 試したこと
|
31
31
|
|
32
32
|
delta_Lxとdelta_Lyを逆にした.
|
33
33
|
|
34
34
|
### 補足情報(FW/ツールのバージョンなど)
|
35
|
-
|
35
|
+
’’’Python
|
36
36
|
delta_Lx = 0.01 # グリッドx幅
|
37
37
|
delta_Ly = 0.01 # グリッドy幅
|
38
38
|
LLx = 1
|
@@ -47,10 +47,10 @@
|
|
47
47
|
phi = np.zeros([Lx+1, Ly+1])
|
48
48
|
phi_Bound = np.zeros([Lx+1, Ly+1])
|
49
49
|
phi_Bound[0, :] = V # 全ての行の0列目を取得
|
50
|
-
|
50
|
+
’’’
|
51
51
|
|
52
52
|
# 電荷密度の設定
|
53
|
-
|
53
|
+
’’’Python
|
54
54
|
eps0 = 1
|
55
55
|
charge = np.zeros([Lx+1, Ly+1])
|
56
56
|
c = 1500
|
@@ -66,10 +66,11 @@
|
|
66
66
|
aa_recta = 0.5 * (np.cos(np.pi / Lx) + np.cos(np.pi / Ly))
|
67
67
|
omega_SOR_recta = 2 / (1 + np.sqrt(1 - aa_recta ** 2)) # SOR法における加速パラメータ
|
68
68
|
print("omega_SOR_rect=", omega_SOR_recta)
|
69
|
-
|
69
|
+
’’’
|
70
70
|
|
71
|
+
|
71
72
|
# メイン
|
72
|
-
|
73
|
+
’’’Python
|
73
74
|
delta = 1.0
|
74
75
|
n_iter = 0
|
75
76
|
conv_check = []
|
@@ -95,7 +96,7 @@
|
|
95
96
|
|
96
97
|
im = plt.imshow(phi, cmap = 'jet')
|
97
98
|
anim.append([im])
|
98
|
-
|
99
|
+
’’’
|
99
100
|
|
100
101
|
前進差分法により離散化をしたときの電場のグラフです.
|
101
102
|

|
2
前進差分法の図を追加した.
title
CHANGED
File without changes
|
body
CHANGED
@@ -95,4 +95,8 @@
|
|
95
95
|
|
96
96
|
im = plt.imshow(phi, cmap = 'jet')
|
97
97
|
anim.append([im])
|
98
|
-
''''''''''''''''''''''''''''''''''''''''''''''''''
|
98
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''
|
99
|
+
|
100
|
+
前進差分法により離散化をしたときの電場のグラフです.
|
101
|
+

|
102
|
+
このように中心差分法の時も描かれないといけませんが, 上記に貼った通り描かれていないのが現状です,
|
1
Pythonで書いたコードの部分の上に''''''Python, 下に'''''''と追加をしました.
title
CHANGED
File without changes
|
body
CHANGED
@@ -8,7 +8,7 @@
|
|
8
8
|

|
9
9
|
|
10
10
|
### 該当のソースコード
|
11
|
-
|
11
|
+
'''''''''''''''''''''''''''''''''''''''''''Python
|
12
12
|
Ex = np.zeros([Lx, Ly])
|
13
13
|
Ey = np.zeros([Lx, Ly])
|
14
14
|
|
@@ -26,11 +26,13 @@
|
|
26
26
|
|
27
27
|
plt.xlim([0, Lx]) # 描くXの範囲
|
28
28
|
plt.ylim([0, Ly]) # 描くYの範囲
|
29
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''
|
29
30
|
### 試したこと
|
30
31
|
|
31
32
|
delta_Lxとdelta_Lyを逆にした.
|
32
33
|
|
33
34
|
### 補足情報(FW/ツールのバージョンなど)
|
35
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''Python
|
34
36
|
delta_Lx = 0.01 # グリッドx幅
|
35
37
|
delta_Ly = 0.01 # グリッドy幅
|
36
38
|
LLx = 1
|
@@ -45,8 +47,10 @@
|
|
45
47
|
phi = np.zeros([Lx+1, Ly+1])
|
46
48
|
phi_Bound = np.zeros([Lx+1, Ly+1])
|
47
49
|
phi_Bound[0, :] = V # 全ての行の0列目を取得
|
50
|
+
'''''''''''''''''''''''''''''''''''''''''''''''''''
|
48
51
|
|
49
52
|
# 電荷密度の設定
|
53
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''Python
|
50
54
|
eps0 = 1
|
51
55
|
charge = np.zeros([Lx+1, Ly+1])
|
52
56
|
c = 1500
|
@@ -62,10 +66,10 @@
|
|
62
66
|
aa_recta = 0.5 * (np.cos(np.pi / Lx) + np.cos(np.pi / Ly))
|
63
67
|
omega_SOR_recta = 2 / (1 + np.sqrt(1 - aa_recta ** 2)) # SOR法における加速パラメータ
|
64
68
|
print("omega_SOR_rect=", omega_SOR_recta)
|
69
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''''
|
65
70
|
|
66
|
-
|
67
|
-
|
68
71
|
# メイン
|
72
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''Python
|
69
73
|
delta = 1.0
|
70
74
|
n_iter = 0
|
71
75
|
conv_check = []
|
@@ -90,4 +94,5 @@
|
|
90
94
|
n_iter += 1
|
91
95
|
|
92
96
|
im = plt.imshow(phi, cmap = 'jet')
|
93
|
-
anim.append([im])
|
97
|
+
anim.append([im])
|
98
|
+
''''''''''''''''''''''''''''''''''''''''''''''''''
|