A simple and rapid method for calculating identity-by-descent matrices using multiple markers