matrixWrite.H

Go to the documentation of this file.

// Investigate the Coefficient Matrix
label Nc = mesh.nCells(); //Total number of cells
simpleMatrix<scalar> A(Nc); //Coefficient matrix
// Initialization of matrix
for(label i=0; i<Nc; i++)
{
    A.source()[i] = 0.0;
    for(label j=0; j<Nc; j++)
    {
        A[i][j] = 0.0;
    }
}
// Assigning diagonal coefficients
for(label i=0; i<Nc; i++)
{
    A[i][i] = TEqn.diag()[i];
}
// Assigning off-diagonal coefficients
for(label faceI=0; faceI<TEqn.lduAddr().lowerAddr().size(); faceI++)
{
    label l = TEqn.lduAddr().lowerAddr()[faceI];
    label u = TEqn.lduAddr().upperAddr()[faceI];
    A[l][u] = TEqn.upper()[faceI];
    A[u][l] = TEqn.upper()[faceI];
}
// Assigning contribution from BC
forAll(T.boundaryField(), patchI)
{
    const fvPatch &pp = T.boundaryField()[patchI].patch();
    forAll(pp, faceI)
    {
        label cellI = pp.faceCells()[faceI];
        A[cellI][cellI] += TEqn.internalCoeffs()[patchI][faceI];
        A.source()[cellI] += TEqn.boundaryCoeffs()[patchI][faceI];
    }
}
 
Info<< "¥n==Coefficient Matrix==" << endl;
for(label i=0; i<Nc; i++)
{
for(label j=0; j<Nc; j++)
{
Info<< A[i][j] << " ";
}
Info<< A.source()[i] << endl;
}
//Info<< A << endl;
Info<< "Solution: " << A.solve() << endl;