Implementation of different computational variations of biconjugate residual methods