[cliffoi]

Hi all,
I am looking into a possible block solver implementation for a set of closely coupled equations in FOAM 1.3. For the time being I am trying to understand how the fvMatrix, lduMatrix and lduAddressing classes work. I have started with a simple problem, trying to reconstruct the coefficient matrix of an existing fvMatrix object using its lduAddressing object. I then pack this new matrix into a simpleMatrix object, copy the source term across, and solve for the field. However, I find that the solution does not match the fvMatrix::solve() solution.
Am I missing something here; the model is not decomposed, the matrix is symmetric and the mesh is symmetric.

Thanks