質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
86.12%
NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

解決済

【python】0除算エラー回避の仕方

mtnow
mtnow

総合スコア13

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

2回答

0グッド

0クリップ

13412閲覧

投稿2019/09/25 01:21

ψ1,ψ2にK1,K2を代入した際に(エラー文)divide by zero encountered in true_divide
がでてしまうのでzの中にあるz=0を除いてzを代入し,K1,K2およびψ1,ψ2を計算したいです.
K1,K2(0<z)

また,z_depth=5のとき,z=0を除けば,zの要素数len(z)=50,yの要素数len(y)=51,となり,y,zの要素数が変わってしまうのですが,ψに代入できますか?

python

1import math 2import numpy as np 3 4########## 5π = math.pi 6Pn = 2200 7z_depth = 5 8c = 0.5 9z = np.arange(0,(z_depth+0.1)*10)/10 10y_c = np.arange((-z_depth*5),(z_depth+0.1)*5)/10 11y = y_c*c 12########## 13yv, zv = np.meshgrid(y, z) 14 15K1 = ((c+yv)**2)+zv**2 #(z<0)の範囲でK1,K2を決定 16K2 = ((c-yv)**2)+zv**2 #z=0,y=-c,+cのとき,K1,K2が0となるため 17 18#ψ1 = π*(1-((K2/K1)**(0.5)))/K1*((K2/K1)**(0.5))*((2*(K2/K1)**(0.5))+(K1+K2-4*(c**2))/K1)**(0.5) 19 20#ψ2 = π*(1+((K2/K1)**(0.5)))/K1*((K2/K1)**(0.5))*((2*(K2/K1)**(0.5))+(K1+K2-4*(c**2))/K1)**(0.5) 21 22input() 23

以下のような質問にはグッドを送りましょう

  • 質問内容が明確
  • 自分も答えを知りたい
  • 質問者以外のユーザにも役立つ

グッドが多くついた質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

気になる質問をクリップする

クリップした質問は、後からいつでもマイページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

下記のような質問は推奨されていません。

  • 間違っている
  • 質問になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

適切な質問に修正を依頼しましょう。

回答2

0

ベストアンサー

numpy の場合 0除算のデフォルトの挙動はErrorではなくて単なるWarningだった気がしますが、本当にエラー(例外)が発生していますか?

とりあえず、0除算時の挙動を

Python

1print(np.geterr()['divide'])

にて確認してみてください。
この結果が 'worm'(ワーニング)や**'print'(ワーニングを標準出力)、'log'(ワーニングをログ出力)'ignore'(無視)の場合はエラーにはならないはずです。
もし、結果が、
''raise**(例外)などになっている場合は

Python

1# ワーニングを表示するが処理は続行する 2np.seterr(divide='warn') 3```Python``` 4または 5```Python 6# ワーニングも出さない 7np.seterr(divide='ignore')

を行い、0除算時の挙動を変更すると良いかと思います。

https://docs.scipy.org/doc/numpy/reference/generated/numpy.seterr.html
https://docs.scipy.org/doc/numpy/reference/generated/numpy.geterr.html

で、上記の設定を行った際の0除算時の結果ですが、

  • 正の値 / 0 の結果は np.inf
  • 負の値 / 0 の結果は -np.inf
  • 0 / 0 の結果は np.nan

として扱われます。

投稿2019/09/25 02:26

magichan

総合スコア15884

良いと思った回答にはグッドを送りましょう。
グッドが多くついた回答ほどページの上位に表示されるので、他の人が素晴らしい回答を見つけやすくなります。

下記のような回答は推奨されていません。

  • 間違っている回答
  • 質問の回答になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

このような回答には修正を依頼しましょう。

0

  • そもそも0除算にならないようにコードを考え直す
  • 除算の前に0かどうかを判断して処理をする
  • 0除算の例外を補足するようにして後でどうにかする

と、こんなもんでしょうか

投稿2019/09/25 01:26

y_waiwai

総合スコア86013

良いと思った回答にはグッドを送りましょう。
グッドが多くついた回答ほどページの上位に表示されるので、他の人が素晴らしい回答を見つけやすくなります。

下記のような回答は推奨されていません。

  • 間違っている回答
  • 質問の回答になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

このような回答には修正を依頼しましょう。

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
86.12%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問

同じタグがついた質問を見る

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。