def partialDomainMatrix(path, angle): return matrixProduct([transmitMatrix(path), rotationMatrix(angle)])
def domainMatrix(path, angle, radius): return matrixProduct([transmitMatrix(path), rotationMatrix(angle), reflectionMatrix(incidenceAngle(angle), radius)])