0のマイナスべき乗を避けるための例外処理

Question

前提・実現したいこと

プログラム内の関数s1からの戻り値の値を(-0.3)乗しようとした時に
s1からの戻り値が0の時にエラーが出るので，新しく関数avoid_erroを作り，関数avoid_erroの中で条件分岐をして0のマイナスのべき乗を避けるためにプログラムを組んでいます．

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

関数s1からの戻り値が-0.001なので，関数avoid_erroからの戻り値は
-0.001^(-0.3) = -7.943282347242818になるはずなのですが，
avoid_erro関数からの戻り値がnanになってしまいます．

PS C:\Users\□□\Desktop\○○> python error-sample.py
関数s1からの戻り値 -0.001
error-sample.py:36: RuntimeWarning: invalid value encountered in double_scalars
  return erro**(η-1)
関数avoid_erroの戻り値 nan
-7.943282347242818

該当のソースコード

import numpy as np


c1 = 1
c2 = 0.05

η = 0.7

def s1(x1,x2,t):
    return c1*(x1-0.02*np.sin(t))+c2*(x2-0.02*np.cos(t)) 

def avoid_erro(x1, x2, t):
    erro = s1(x1, x2, t)
    print("関数s1からの戻り値", erro)
    if(erro == 0):
        return 0
    else:
        #print(erro**(η-1))
        return erro**(η-1)

print("関数avoid_erroの戻り値",avoid_erro(0, 0, 0))
print(-0.001**(η-1))

試したこと

printを使ってどこで値が変わっているのかを確かめましたがわかりませんでした

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Comments

[yamap55]

あまり詳しくないのでこちらで。違いとしてはs1からの返り値は「numpy.float64」で、-0.001は「float」ですね。
また、浮動小数点数は正しく？値を持てないのでその辺りが原因になるのかと思います。
（erro == -0.001とするとTrueが返ってくるのも嫌な感じですね。これはnpがいい感じに比較してくれてるのだと思いますが。）

※本文の数値が「0.01」となっていますが、正しくは「0.001」ですね。

[meg_]

エラー処理では駄目なのでしょうか？

[NR4200]

yamap_55さんへ
間違いの指摘ありがとうございます．修正しておきます．

meg_さんへ
エラー処理とは，try-exceptの例外処理のことでしょうか？

[tmp]

-0.001^(-0.3)= -7.943282347242818ですが
(-0.001)^(-0.3)は、nan

Answer 1

エラーの原因は計算結果が複素数になるからではないでしょうか。
erro をnumpy.complex型として計算するとエラーが発生しなくなります

import numpy as np

c1 = 1
c2 = 0.05
η = 0.7

def s1(x1,x2,t):
    return c1*(x1-0.02*np.sin(t))+c2*(x2-0.02*np.cos(t))

def avoid_erro(x1, x2, t):
    erro = s1(x1, x2, t)
    print("関数s1からの戻り値", erro)
    if(erro == 0):
        return 0
    else:
        #print(erro**(η-1))
        return np.complex(erro) ** (η-1)

result = avoid_erro(0, 0, 0)
print("関数avoid_erroの戻り値", result)
# 関数avoid_erroの戻り値 (4.668944218504467-6.426250410037962j)
print("関数avoid_erroの戻り値の絶対値", abs(result))
# 関数avoid_erroの戻り値の絶対値 7.943282347242818
print(-0.001**(η-1))
# -7.943282347242818

Answer 2

プログラムを見直すと（）が抜けていたため
(-0.01)^0.7がをしたいのに0.01^0.7になっていました

tmpさんありがとうございました