(zsph <- gell(Sigma = diag(3))) # unit sphere in R^3
(zplane <- gell(span = diag(3)[,1:2])) # a plane
dual(zplane) # line orthogonal to that plane
(zhplane <- gell(center = c(0,0,2), span = diag(3)[,1:2])) # a hyperplane
dual(zhplane) # orthogonal line through same center (note that the 'gell'
# object with a center contains more information than the geometric plane)
zorigin <- gell(span = cbind(c(0,0,0)))
dual( zorigin )
# signatures of these ellipsoids
signature(zsph)
signature(zhplane)
signature(dual(zhplane))
Run the code above in your browser using DataLab