if theta0 is None:
            theta0 = theta00