It was my misunderstanding coming from Matrix.norm():
        ord    norm for matrices             norm for vectors
        =====  ============================  ==========================
        None   Frobenius norm                2-norm
        'fro'  Frobenius norm                - does not exist
        inf    --                            max(abs(x))

I will replace us with S.Infinity in test_P8()