Now that each operator is verified for single column vectors, all that remains is to verify proper action on multiple columns of input. To check cI()'s action on multiple columns, perform the following test at the octave prompt:
in=randn(prod(S),3)+i*randn(prod(S),3); %# Form random input with 3 columns out1=cI(in); %# Output of new operator out2=[cI(in(:,1)), cI(in(:,2)), cI(in(:,3))]; %# Output using debugged case max(abs(out2-out1)) %# Check maximum value for discrepancy in each columnFinally, you should repeat the test, substituting each of the remaining operators for cI in the code fragment above.