用cupy（python）计算矩阵的行列式

Posted 2023-02-16

【中文标题】用cupy（python）计算矩阵的行列式

【英文标题】：calculating the determinant of a matrix with cupy (python)

【发布时间】：2020-05-14 14:14:58

【问题描述】：

我目前正在编辑 Python 中用于 GPU 加速的 Cupy 代码。大多数编辑都需要从 numpy 直接交换到 cupy（例如 cp.array 而不是 np.array），但有一些小例外。但是，我在这段代码中遇到了问题，它可以与 numpy 一起使用，我不确定为什么我突然收到这个错误。我有以下代码

ab = (b - a)
ac = (c - a)
ap1 = (p1 - a)
ap2 = (p2 - a)

其中 a,b,c,p1 和 p2 都是一维 Cupy 数组

从这里我想使用以下来计算行列式。

testP1 = 1 if cp.linalg.det([ab,ac,ap1]) >0 else 0
testP2 = 1 if cp.linalg.det([ab,ac,ap2]) >0 else 0

这在 numpy 中工作正常（输入是 numpy 而不是 cupy 数组），但是，我现在收到以下错误

AttributeError: 'list' 对象没有属性 'ndim'

【问题讨论】：

不要使用列表作为det 的输入。使用一个或 numpy 或 cupy 数组

您是否尝试将[ab,ac,ap1] 显式转换为cp.array 对象？也许 numpy 隐含地进行转换，而 cupy 没有。

嗨@g2i 就是这样，在尝试最简单的选项之前，我愚蠢地尝试了所有其他组合，包括在数组类型之间传输。谢谢

嗨，MattC1990，在写下一个问题之前，请先熟悉asking questions！编码愉快！

【参考方案1】：

正如 g2i 所评论的，cupy 不会隐式转换为数组，而 numpy 会，因此必须显式转换为 cupy 数组