Matrix multiplication of pyomo variable with constant matrix.

I have the pyomo variable defined in the form of matrix, let’s say 21×201.

I want to add a constraint that matrix multiplication of pyomo variables such as :
indx = 11 ; A is diagonal matrix containing the number 10 at the diagonal:
x[1:21, 2:indx] = A * x[1:21, 1:10]

How can I achieve this in python and using pyomo model.

Thanking in advance.