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

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

ただいまの
回答率

89.07%

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

解決済

回答 2

投稿

  • 評価
  • クリップ 0
  • VIEW 3,089

mtnow

score 13

ψ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の要素数が変わってしまうのですが,ψに代入できますか?

import math
import numpy as np

##########
π = math.pi
Pn = 2200
z_depth = 5
c = 0.5
z = np.arange(0,(z_depth+0.1)*10)/10
y_c = np.arange((-z_depth*5),(z_depth+0.1)*5)/10
y = y_c*c
##########
yv, zv = np.meshgrid(y, z)

K1 = ((c+yv)**2)+zv**2 #(z<0)の範囲でK1,K2を決定
K2 = ((c-yv)**2)+zv**2 #z=0,y=-c,+cのとき,K1,K2が0となるため

#ψ1 = π*(1-((K2/K1)**(0.5)))/K1*((K2/K1)**(0.5))*((2*(K2/K1)**(0.5))+(K1+K2-4*(c**2))/K1)**(0.5)

#ψ2 = π*(1+((K2/K1)**(0.5)))/K1*((K2/K1)**(0.5))*((2*(K2/K1)**(0.5))+(K1+K2-4*(c**2))/K1)**(0.5)

input()
  • 気になる質問をクリップする

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

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

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

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 2

checkベストアンサー

0

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

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

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

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

# ワーニングを表示するが処理は続行する
np.seterr(divide='warn')

Python``` または

# ワーニングも出さない  
np.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

として扱われます。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

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

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

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

0

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

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

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

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

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

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

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

  • ただいまの回答率 89.07%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

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