`numpy.linalg._linalg:_multi_dot_matrix_chain_order`

Signature

def _multi_dot_matrix_chain_order ( arrays , return_costs = False )

Return a np.array that encodes the optimal order of multiplications.

The optimal order array is then used by _multi_dot() to do the multiplication.

Also return the cost matrix if return_costs is True

The implementation CLOSELY follows Cormen, "Introduction to Algorithms", Chapter 15.2, p. 370-378. Note that Cormen uses 1-based indices.

cost[i, j] = min([
cost[prefix] + cost[suffix] + cost_mult(prefix, suffix) for k in range(i, j)])