fn accumulate_product( target: &mut DMatrix<Complex>, left: MatRef<'_, Complex>, right: MatRef<'_, Complex>, )
target += left * right.
target += left * right