Acceleration of 3D pin-by-pin calculations based on the heterogeneous variational nodal method