i.e., linear model M and linear observation operator H
n = ub.shape[0]
```

```