Ask Your Question

Using sage to check if a vector is in a matrix's null space

asked 2013-11-26 06:04:34 -0500

ensaba gravatar image

I want to check if a given vector is inside a matrix's null space.

Here is what I have done:

#Define a matrix

#Free variables to be used in the given vector
var('s t')

#The null space of the matrix A


#Create a vector

#If ans is in the null space of A, then A*ans is the 0 vector
print A*ans

#If A*ans is the 0 vector, then the following should be true too
print ans in ns

I had expected the last print statement to be true since A*ans results in a zero vector, hence ans is in the null space of A.

Why isn't this the case?

edit retag flag offensive close merge delete

2 answers

Sort by ยป oldest newest most voted

answered 2013-11-26 06:17:44 -0500

updated 2013-11-26 06:18:45 -0500

Your vector ans is not in A.kernel() because of how A.kernel() is defined. From its documentation:

Returns the left kernel of this matrix, as a vector space or free
module. This is the set of vectors "x" such that "x*self = 0".

Note: For the right kernel, use "right_kernel()".  The method
"kernel()" is exactly equal to "left_kernel()".

So since ans * A is not zero, ans is not in A.kernel().

sage: ans in A.kernel()
sage: ans in A.right_kernel()

A.right_kernel() is what you want.

edit flag offensive delete link more

answered 2013-11-26 06:14:35 -0500

tmonteil gravatar image

Actually, the .kernel() method corresponds to the left kernel. What you are looking for is the right kernel:

sage: ns = A.right_kernel()
sage: ans in ns
edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2013-11-26 06:04:34 -0500

Seen: 577 times

Last updated: Nov 26 '13