gl.grm function

Calculates an identity by descent matrix