Ainv(iA,jA) = Ainv(iA,jA) + Dinv(iA,kA)*CB(kA,jA)
_RL Ainv(3,3), Binv(3,3), Cinv(3,3), Dinv(3,3), CB(3,3)
Dinv(1,1) = COS(phiRad)
Dinv(1,2) = -SIN(phiRad)
Dinv(1,3) = 0. _d 0
Dinv(2,1) = SIN(phiRad)
Dinv(2,2) = COS(phiRad)
Dinv(2,3) = 0. _d 0
Dinv(3,1) = 0. _d 0
Dinv(3,2) = 0. _d 0
Dinv(3,3) = 1. _d 0