MILP initialization for solving parabolic PDEs with PINNs