add the possibility to get the force of contact tasks computed by the solver during optimization process