MULTIMODECODE: an efficient numerical solver for multifield inflation