"where M is a square n1??n2??n3 by n1??n2??n3 matrix and r is a known vector of length n1??n2??n3." . . . .