0

Given a gsl_matrix * A object, which is an M-by-N matrix, what is the easiest way to find the kernel of A? I tried using singular value decomposition (specifically the gsl_linalg_SV_decomp method), but that doesn't work when M<N, which is sometimes the case for my matrices.

rhashimoto
  • 15,650
  • 2
  • 52
  • 80
Izzhov
  • 163
  • 4

1 Answers1

1

The easiest way (if you are restricted to that library) would be to pad your matrix with rows of zeros to N x N, which you can then pass to your function. The padded matrix will have the same null space.

rhashimoto
  • 15,650
  • 2
  • 52
  • 80