message: Overwrite inputs to `LinearOperatorBlockLowerTriangular` methods instead of modifying them in-place.
diff: @@ -726,7 +726,7 @@ class LinearOperatorBlockLowerTriangular(linear_operator.LinearOperator):
  # `Ax_0 + Bx_1 + Dx_2 = y_0` as `Ax_0 = y_0*`, where
  # `y_0* = y_0 - Bx_1 - Dx_2`.
  for j in reversed(range(index + 1, len(self.operators))):
  - y -= self.operators[j][index].matmul(
  + y = y - self.operators[j][index].matmul(
  solution_list[len(self.operators) - 1 - j],
  adjoint=adjoint)
  # Continuing the example above, solve `Ax_0 = y_0*` for `x_0`.
  @@ -748,7 +748,7 @@ class LinearOperatorBlockLowerTriangular(linear_operator.LinearOperator):
  # `Dx_0 + Ex_1 + Fx_2 = y_2` as `Fx_2 = y_2*`, where
  # `y_2* = y_2 - D_x0 - Ex_1`.
  for i, operator in enumerate(row[:-1]):
  - y -= operator.matmul(solution_list[i], adjoint=adjoint)
  + y = y - operator.matmul(solution_list[i], adjoint=adjoint)
  # Continuing the example above, solve `Fx_2 = y_2*` for `x_2`.
  solution_list.append(row[-1].solve(y, adjoint=adjoint))
  