* colnew.solve:
  - Allow passing initial mesh + mesh values as the initial guess
* mus:
  - Do we really want to include this?
