python - 错误的 NumPy 数组计算结果

Question

我一直试图弄清楚为什么它不起作用。必须说我是python新手。尝试计算何时 a^d + b^d = c^d; （费马定理）使用 np.arrays 以获得最佳性能。问题：数组 B 中 d > 2 的结果错误。请指出正确的方向？谢谢！

Kubuntu 17.10 Python 3.6.4 |Anaconda 自定义（64 位）| （默认，2018 年 1 月 16 日，18:10:19）键入“copyright”、“credits”或“license”以获取更多信息。IPython 6.4.0——增强的交互式 Python。

这是问题显示时的代码示例：

import numpy as np

d = np.float64(4.)
A = np.float64([[99001457116.0, 99536254152.0],
                [99002341176.0, 99914925909.0],
                [99029463918.0, 99750866948.0]])

B = np.float64((np.sum(A**2, axis=1)**(1/d)).reshape((A.shape[0], 1)))
B

标签： pythonarraysnumpypython-3.6calculus