连老师，您好。冒昧打扰您，想请教您一个问题，是关于用stata做分块矩阵组内乘法的

我有两个矩阵A和B，行、列分别是ab*b 和ab*1 （a=30，b=150）。A矩阵由a个b*b的子矩阵mi组成，B矩阵由a个b*1的子矩阵ni组成（i=1,2，……，150）。其中，a可以由变量"province"来识别，b可以由变量“industry”来识别。我想让每个mi和ni做组内的矩阵乘法，但不知道应该怎么编程。

原始数据不是“矩阵”格式的，而是“变量”格式的，所以一开始我想试着先用reshape命令再进行矩阵运算。可是因为a和b的取值都比较大，所以不太可行。请教一下连老师，能不能通过编程实现循环的矩阵乘法运算

A: a(1,1)    a(1,2 ) ...       a(1,150)       B:   b(1,1)  

   ...                                          ..... .  

   a(150,1)   a(150,2) ...       a(150,150)         b(150,1)

a(151,1)   a(151,2) ...       a(151,150)        b(151,1)                                     

.......                                                                        .......   

a(300,1)   a(300,2) ...          a(150,150)              b(300,1)       

........................................         .......               

a(4351,1)   a(4351,2) ...           a(4351,150)                b(4351,1)                                   

.................................   .           ........   

a(4500,1)   a(4500,2) ...          a(4500,150)                 b(4500,1)