hat_matrix_cpp function

Computing single H matrix in AICc-function using the Mahalanobis distance