MatMult: norm difference < sqrt(eps_machine)
              MatMultTranspose: norm difference < sqrt(eps_machine)
             MatMultAdd v2==v3: norm difference < sqrt(eps_machine)
             MatMultAdd v2!=v3: norm difference < sqrt(eps_machine)
    MatMultTransposeAdd v2==v3: norm difference < sqrt(eps_machine)
    MatMultTransposeAdd v2!=v3: norm difference < sqrt(eps_machine)