MHJacNL is renamed MHBilinear, and allows to handle multiple (multi-harmonic) arguments in the function. It can also assemble contributions either in the system matrix instead of the Jacobian.