 constant vector, x is the solution vector.
  #
  # ps is the pivot, a vector which indicates the permutation of rows 