Doing %o% that operates on columns instead of atomics

my.foo <- function( a, b ) { c <- a - b; }  #really more complex, but just
to illustrate