Add substitute function for molecules.
Fix unittest.
Merge branch 'master' into mol_edit
methods to get bond length. Auto setting of bond lengths for substitutions.
Better check for convergence.