solving for nonzero x in Ax=0

[mkrems]

is there a specific lapack or scalapack routine or a general way of doing this?

thanks!

[buttari]

is there a specific lapack or scalapack routine or a general way of doing this?
thanks!

Hello,
a singular value decomposition will return you an orthogonal basis for the null space of the matrix. You can do this with the DGESVD subroutine in LAPACK (or the PDGESVD in ScaLAPACK).

Regards

Alfredo