and nonlinear observation operator
n = ub.shape[0]