Implement a strategy to allow operators to not support the 'entry=(r,c)' syntax

Note: The solution is neither pretty nor efficient.
We can not fix this 'the right way' for now.
Fix the matrix builder call to work in homogeneous and inhomogeneous case
Fix some documentation
Fix observable computations