def _threeDdot_simple(M,a):
result = np.empty(a.shape,dtype=a.dtype)<EOL>for i in range(a.shape[<NUM_LIT:0>]):<EOL><INDENT>for j in range(a.shape[<NUM_LIT:1>]):<EOL><INDENT>A = np.array([a[i,j,<NUM_LIT:0>],a[i,j,<NUM_LIT:1>],a[i,j,<NUM_LIT:2>]]).reshape((<NUM_LIT:3>,<NUM_LIT:1>))<EOL>L = np.dot(M,A)<EOL>result[i,j,<NUM_LIT:0>] = L[<NUM_LIT:0>]<EOL>result[i,j,<NUM_LIT:1>] = L[<NUM_LIT:1>]<EOL>result[i,j,<NUM_LIT:2>] = L[<NUM_LIT:2>]<EOL><DEDENT><DEDENT>return result<EOL>
Return Ma, where M is a 3x3 transformation matrix, for each pixel