init_fp_diffmatrix Subroutine