An answer to this question claims that scipy.odr should compute an orthogonal regression and that it will be equal to the total least squares solution.
You can use scipy.odr it will compute orthogonal regression which should be equal to tls.
However Wikipedia claims that both demming regression and orthogonal regression are special cases of total least squares.
It is a generalization of Deming regression and also of orthogonal regression, and can be applied to both linear and non-linear models.
The SciPy documentation cites Orthogonal distance regression by Boggs and Donaldson 1989, which does not appear to mention total least squares.
Is using scipy.odr actually equivalent to this general case of total least squares?