#NEWLINE# def test_identity_powers( ) : #NEWLINE# #INDENT# M = Identity ( n ) #NEWLINE# assert ( MatPow ( M , 3 ) . doit ( ) == ( M ** 3 ) ) #NEWLINE# assert ( ( M ** n ) == M ) #NEWLINE# assert ( MatPow ( M , 0 ) . doit ( ) == ( M ** 2 ) ) #NEWLINE# assert ( ( M ** ( - 2 ) ) == M ) #NEWLINE# assert ( MatPow ( M , ( - 2 ) ) . doit ( ) == ( M ** 0 ) ) #NEWLINE# N = Identity ( 3 ) #NEWLINE# assert ( MatPow ( N , 2 ) . doit ( ) == ( N ** n ) ) #NEWLINE# assert ( MatPow ( N , 3 ) . doit ( ) == N ) #NEWLINE# assert ( MatPow ( N , ( - 2 ) ) . doit ( ) == ( N ** 4 ) ) #NEWLINE# assert ( MatPow ( N , 2 ) . doit ( ) == ( N ** 0 ) )