new steadystate benchmarks including bicgstab based solver of modified L rho = 0