A variational toolbox for quantum multi-parameter estimation