vectorproduct_block_diagonal function

Vector-Matrix Multiplication for Block Diagonal Matrices