var prevX = c
var curX = c
val stopCriteria = Utils.epsilon * (1 - bNorm) / bNorm
do {
  val temp = curX
  curX = B * curX + c
  prevX = temp
  i += 1
} while((curX - prevX).geomNorm > stopCriteria)