A novel matrix multiplication framework for modeling genotype-by-environment interaction in genomic prediction