extend MatrixProductOperator for multi-dim fields

9 jobs for matrix_product_operator_mf in 17 minutes and 6 seconds (queued for 2 seconds)