p { identity r
C generate the identity matrix, e.g. identity 3
p { (Ir) J.= Ir