[toboto]

Dear all,

To the best of my knowledge the lduMatrix class consists of 5 arrays: diag(), lower(), upper(), lowerAddr() and upperAddr().

in the lower() array the lower triangle coefficients are stored. My question is how are they stores, i.e. in which order?

In the post here is said that the coefficient is stored row wised. Although this is inconsistent in my point of view because for any face k the owner and neighbour of that face should be:
lowerAddr(k), upperAddr(k).
which means that this coefficients should be arranged in an ascending order based on the face they belong to.

Any clarification?
Best regards.