Jupyterでソースコードを動かしたい

Question

前提・実現したいこと

発生している問題・エラーメッセージ

Jupyterでプログラムを動かそうとしたところ、In[*]となってしまい、プログラムを動作、描画することができませんでした。初学者ですのでわからないことが多いですが、どうぞよろしくお願いします。

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D


LIMIT = 1000 # 反復回数の上限
N = 101 # x軸方向の分割数
M = 101 # y軸方向の分割数


def iteration(u):
    """1回分の反復計算"""
    u_next = [[0 for i in range(N)] for j in range(M)] # 次ステップのuij
    for i in range(1, N - 1):
        for j in range(1, M - 1):
            u_next[i][j] = (u[i][j - 1] + u[i -1][j] + u[i + 1][j] + u[i][j + 1]) / 4
            for i in range(1, N - 1):
                for j in range(1, M - 1):
                    u[i][j] = u_next[i][j]


u = [[0 for i in range(N)] for j in range(M)] # uijの初期化
for i in range(M):
    u[i][0] = 100
    

for i in range(LIMIT):
    iteration(u)
print(u) # 結果の出力


x = np.arange(0, N)
y = np.arange(0, M)
X, Y = np.meshgrid(x, y)
fig = plt.figure()
ax = Axes3D(fig)
U = np.array(u)
#ax.plot_wireframe(X, Y, U) # wireframe形式
surf = ax.plot_surface(X, Y, U, cmap = "plasma", alpha = 0.8) # surface形式

fig.colorbar(surf)
plt.show()

Python

Answer 1

こんにちは。
関連した質問がありましたので、もしかしたら参考になるかもしれません。

https://teratail.com/questions/107521
https://ja.stackoverflow.com/questions/64014/pythonのセルが実行出来ません

どちらも Kernel -> Restart で解消されているようです。

添付いただいたソースコードはインデントが省略されてしまっているので、コード面からは現時点でコメントができないのですが、まずは上記を参考にしてみて、その後エラーが違うものになったら、スクリーンショットやメッセージを添えていただければと思います。

Comments

[aa-]

ご回答ありがとうございます。そちらの質問は拝見させていただきましたが、同様の操作でもうまくいきませんでした。

[suama]

ありがとうございます。お試しいただいた上ですね。失礼いたしました。
ソースコードも再添付してくださってありがとうございます！

[suama]

単純にデータ量が多くて集計に時間がかかっているのかもしれませんね。

LIMIT = 1000 とありますが、まずはここを1とか5くらいに減らすとどうでしょう？コード上のエラーか、処理の速度の問題かの切り分けのため。

[aa-]

ありがとうございます。試してみましたが、計算結果は表示されませんでした。プログラムが実行されずに次の行にいってしまいます。悲しい????

[aa-]

何か使えないツールが自分の環境であるのかな～とも思い、調べてみてはいるのですが、それが何かはわからないので困ったものです。

[suama]

Jupyter あまり詳しくなくてもうしわけないですが、処理の数と環境でしょうか。
試しにご提示いただいたソースコードを、 https://https://jupyter.org/try で実行してみました。

以下のように定数を減らすと、10秒以内にデータの配列が出力されて、グラフも描画されました。
-----
LIMIT = 1000 # 反復回数の上限
N = 10 # x軸方向の分割数
M = 10 # y軸方向の分割数
-----
(値は適当に小さくしています。申し訳ありません)
ひとまず　plt.show() まではエラーなく進んだのは確認できています。

定数を多くしていくと、画面したの Busy の表示がなかなか消えません。
処理が終わると idle に変わります。
（ただし、もともとの数だと、20分待っても終わらないです）

参考になりましたら...。

[aa-]

ありがとうございます。提示していただいた通り定数を減らしたところうまくいきました。単純に処理時間の問題だったようです。