There was no problem with the multimode Jacobian! Just need to select initial c near threshold to be large enough not to make Newton oscillate unstably.