Implemented optimization in AugmentedLagrangianUnitFlow algorithm when
both flow and right hand side are linear operators

Here we do not need to use secant for unit flow iteration. We compute
the flowrate for normalized rhs and then using linearity to find the
proper parameter of rhs