Does right_kernel_matrix support assumptions over the symbolic ring?

asked 2016-08-09 09:51:57 -0600

exodus gravatar image

I'm trying to a basis for the kernel of a symbolic matrix. However, exactly quite what the kernel is depends on the assumptions that have been made. It seems that although other matrix functions (e.g. calculating the inverse) respect assumptions the kernel function doesn't. Here is an example:

var('x,y')

M = Matrix(SR, [[x, 0],[0,y]])
print(~M) # Gives expected inverse [[1/x, 0], [0,1/y]]
print(M.right_kernel_matrix()) # Gives [] as expected (if x, y != 0 then M is invertible so has trivial kernel)

assume(x == 0)
# print(~M) # Throws an error (ZeroDivisionError)
print(M.right_kernel_matrix()) # Still gives [], whereas the kernel now has basis [(1, 0)]

Am I missing something here? Is it possible to use this function with assumptions?

edit retag flag offensive close merge delete

Comments

Most likely computing the inverse of a symbolic matrix is done using Maxima, which knows about the assumptions, while the other methods don't, so it doesn't. (More or less only methods which use Maxima will know about them; usually one doesn't have it used for matrix stuff so this is an outlier. Yes, Sage's assumption framework could use improvement.)

kcrisman gravatar imagekcrisman ( 2016-08-09 20:57:35 -0600 )edit