Case study of a differentiable heterogeneous multiphysics solver for a nuclear fusion application