init.compute_betas_from_shifts initialize the vector of optimal values at nodes and
tips, with the value at the root.
init.compute_betas_from_shifts(phy, optimal.value, ...)compute_betas_from_shifts_from_shifts and is designed to
furnish function update.compute_betas_from_shifts with the right structure of data.