ceres-solver / ceres-solver / 8e3b7d89e4a995064b5d1a105f21fc6ea0ba9ef9 / . / docs / source / modeling_faqs.rst

.. _chapter-modeling_faqs: | |

.. default-domain:: cpp | |

.. cpp:namespace:: ceres | |

======== | |

Modeling | |

======== | |

#. Use analytical/automatic derivatives. | |

This is the single most important piece of advice we can give to | |

you. It is tempting to take the easy way out and use numeric | |

differentiation. This is a bad idea. Numeric differentiation is | |

slow, ill-behaved, hard to get right, and results in poor | |

convergence behaviour. | |

Ceres allows the user to define templated functors which will | |

be automatically differentiated. For most situations this is enough | |

and we recommend using this facility. In some cases the derivatives | |

are simple enough or the performance considerations are such that | |

the overhead of automatic differentiation is too much. In such | |

cases, analytic derivatives are recommended. | |

The use of numerical derivatives should be a measure of last | |

resort, where it is simply not possible to write a templated | |

implementation of the cost function. | |

In many cases it is not possible to do analytic or automatic | |

differentiation of the entire cost function, but it is generally | |

the case that it is possible to decompose the cost function into | |

parts that need to be numerically differentiated and parts that can | |

be automatically or analytically differentiated. | |

To this end, Ceres has extensive support for mixing analytic, | |

automatic and numeric differentiation. See | |

:class:`CostFunctionToFunctor`. | |

#. When using Quaternions, consider using :class:`QuaternionManifold`. | |

`Quaternions <https://en.wikipedia.org/wiki/Quaternion>`_ are a | |

four dimensional parameterization of the space of three dimensional | |

rotations :math:`SO(3)`. However, the :math:`SO(3)` is a three | |

dimensional set, and so is the tangent space of a | |

Quaternion. Therefore, it is sometimes (not always) beneficial to | |

associate a local parameterization with parameter blocks | |

representing a Quaternion. Assuming that the order of entries in | |

your parameter block is :math:`w,x,y,z`, you can use | |

:class:`QuaternionManifold`. | |

.. NOTE:: | |

If you are using `Eigen's Quaternion | |

<http://eigen.tuxfamily.org/dox/classEigen_1_1Quaternion.html>`_ | |

object, whose layout is :math:`x,y,z,w`, then you should use | |

:class:`EigenQuaternionManifold`. | |

#. How do I solve problems with general linear & non-linear | |

**inequality** constraints with Ceres Solver? | |

Currently, Ceres Solver only supports upper and lower bounds | |

constraints on the parameter blocks. | |

A crude way of dealing with inequality constraints is have one or | |

more of your cost functions check if the inequalities you are | |

interested in are satisfied, and if not return false instead of | |

true. This will prevent the solver from ever stepping into an | |

infeasible region. | |

This requires that the starting point for the optimization be a | |

feasible point. You also risk pre-mature convergence using this | |

method. | |

#. How do I solve problems with general linear & non-linear **equality** | |

constraints with Ceres Solver? | |

There is no built in support in ceres for solving problems with | |

equality constraints. Currently, Ceres Solver only supports upper | |

and lower bounds constraints on the parameter blocks. | |

The trick described above for dealing with inequality | |

constraints will **not** work for equality constraints. | |

#. How do I set one or more components of a parameter block constant? | |

Using :class:`SubsetManifold`. |