How to stack same rows of prior distributions inside the with model effectively

@ckrapu Hi Chris, sorry to bother you again. The code I have below uses for loop to dot product of every two columns with beta vector, totally 8 columns in the matrix. I wonder whether this is better way to do this. This is inside the with part. Thanks a lot.

   # Define utility for fixed effects
    utility_fixed = tt.dot(X_fixed[:,[0,1]], beta)
    utility_fixed = tt.reshape(utility_fixed, (69,1))
           
    for i in range(3):
        part = X_fixed[:, [2 + i, 3 + i]]
        dot_prod = tt.dot(part, beta)
        dot_prod = tt.reshape(dot_prod, (69,1))
        utility_fixed = tt.concatenate([utility_fixed, dot_prod], axis = 1)