Use matrix provided by computeLinearDecompositionA() for inner product
in case of using EIM or non symmetric bilinear form a()
Linked to #339