[Jerryfan]

Hi foamers,



I have some difficulties in understanding the losort addressing in the lduMatrix. Since we already have the lower, upper and diagonal addresses in the ldu address system, what is the purpose of defining this losort address? Can anyone help explain this? Thanks.

[Jerryfan]

I got it. It's a list that holds the face index in which the cell appears as the neighbor cell for that face.

[Jerryfan]

A good example of using this is in the DILU preconditioner.

for (register label face=0; face<nFaces; face++)
{
sface = losortPtr[face];
wAPtr[uPtr[sface]] -=
rDPtr[uPtr[sface]]*lowerPtr[sface]*wAPtr[lPtr[sface]];
}