setting BC position before or after solver affects iteration in solver

[YangHaibin1102]

In UWGeo,

if I set temperatureBCs before Model.solver.set_inner_method("mumps"), it works well

Model.set_temperatureBCs(
top = 273.15 * u.degK,
bottom = 1723.15 * u.degK,
)

Model.solver.set_inner_method("mumps")

but if I reverse the position, putting temperatureBCs after Model.solver.set_inner_method("mumps"), it fails.

Shear3D_nonlinear_UWG.py