teratail header banner
teratail header banner
質問するログイン新規登録

質問編集履歴

5

```にしました

2021/01/16 10:01

投稿

nagi-kimura
nagi-kimura

スコア3

title CHANGED
File without changes
body CHANGED
@@ -8,7 +8,7 @@
8
8
  ![イメージ説明](781075d5452af12c1fdbd3d97bfe49a1.png)
9
9
 
10
10
  ### 該当のソースコード
11
- ’’’Python
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
- ’’’Python
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
- ’’’Python
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
- ’’’Python
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
  ![イメージ説明](156319d1c8da9bb316adae636ae5ced1.png)
104
101
  このように中心差分法の時も描かれないといけませんが, 上記に貼った通り描かれていないのが現状です,

4

intendしました

2021/01/16 10:01

投稿

nagi-kimura
nagi-kimura

スコア3

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

’’’としました

2021/01/16 09:42

投稿

nagi-kimura
nagi-kimura

スコア3

title CHANGED
File without changes
body CHANGED
@@ -8,7 +8,7 @@
8
8
  ![イメージ説明](781075d5452af12c1fdbd3d97bfe49a1.png)
9
9
 
10
10
  ### 該当のソースコード
11
- '''''''''''''''''''''''''''''''''''''''''''Python
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
- ''''''''''''''''''''''''''''''''''''''''''''''''''Python
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
- ''''''''''''''''''''''''''''''''''''''''''''''''''Python
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
- ''''''''''''''''''''''''''''''''''''''''''''''''''Python
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
  ![イメージ説明](156319d1c8da9bb316adae636ae5ced1.png)

2

前進差分法の図を追加した.

2021/01/16 09:41

投稿

nagi-kimura
nagi-kimura

スコア3

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
+ ![イメージ説明](156319d1c8da9bb316adae636ae5ced1.png)
102
+ このように中心差分法の時も描かれないといけませんが, 上記に貼った通り描かれていないのが現状です,

1

Pythonで書いたコードの部分の上に''''''Python, 下に'''''''と追加をしました.

2021/01/16 09:27

投稿

nagi-kimura
nagi-kimura

スコア3

title CHANGED
File without changes
body CHANGED
@@ -8,7 +8,7 @@
8
8
  ![イメージ説明](781075d5452af12c1fdbd3d97bfe49a1.png)
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
+ ''''''''''''''''''''''''''''''''''''''''''''''''''