(zplane <- gell(span = diag(3)[,1:2])) # a plane
dual(zplane) # orthogonal line
(zplane2 <- gmult( cbind( c(1,1,1), c(1,-1,0), c(1,0,-1)), zplane))
# correctly wipes out one dimension
(zplane3 <- gmult( cbind( c(1,0,0), c(1,0,0), c(0,0,1)), zplane))
Run the code above in your browser using DataLab