python - sage中矩阵元素的条件定义

Question

在 Sage 中，尝试通过以下方式定义具有单元格条件的矩阵:

matrix([[(if gcd(i, j) == 0: log(radical((i+j)*i*j)) else: -1.0) for j in srange(1, 5)] for i in srange(1, 5)])

我收到语法错误:

...
matrix([[(if gcd(i, j) == _sage_const_0 : log(radical((i+j)*i*j)) else: -_sage_const_1p0 ) for j in srange(_sage_const_1 , _sage_const_5 )] for i in srange(_sage_const_1 , _sage_const_5 )])
           ^
SyntaxError: invalid syntax

这里有什么问题吗？如何解决这个问题？

Answer 1

你的问题实际上是 Python 问题，而不是 Sage 本身。 Python 对列表推导式有一些过滤，但它看起来不是这样的。参见例如this question .

所以let's try it :

matrix([[log(radical((i+j)*i*j)) if gcd(i,j)==0 else -1.0 for j in srange(1,5)] for i in srange(1,5)])

顺便问一下，你真的想要if gcd(i,j)==1吗？你不太可能在这个中得到零的 gcd!