Yes, you are right, uij would just become uj. However, the pm.Deterministic call that I am assigning to ‘uj’ is attempting to multiply ‘invL’ by the entire ‘targets’ 2 by 60 matrix. What I want it to do is compute ‘uj’ for each column (trial) of that 2 by 60 matrix. Any idea how this can be implemented?