Adds output param  for GaussEliminate(), to return pivot positions (row, pivot_columns[row]) of coefficient matrix.