[hawkeye321]

Hi FOAMERS

I am trying to access to elements of the coefficient matrix of an equation using lduAddresing

Foloowing the procedure in "OpenFOAM guide/Matrices in OpenFOAM" http://openfoamwiki.net/index.php/Matrices_in_OpenFOAM

I use

scalarField *lowerPtr_;
scalarField *diagPtr_;
scalarField *upperPtr_;

for (register label face=0; face<l.size(); face++)
{
Diag[l[face]] -= Lower[face];
Diag[u[face]] -= Upper[face];
}

But I get the following errors

error: ‘l’ was not declared in this scope
error: ‘Diag’ was not declared in this scope
error: ‘Lower’ was not declared in this scope
error: ‘u’ was not declared in this scope
error: ‘Upper’ was not declared in this scope

Do you have any idea what these errors are?
PS: I have included these header files

#include "fvCFD.H"
#include "wallFvPatch.H"
#include <fstream>
#include <string>
#include <stdio.h>
#include <time.h>
#include <iostream>

[fumiya]

Hi,

If your equation is UEqn, you can access the diag, upper and lower
coefficients of UEqn in the following manner:

    label Nc = mesh.nCells(); //Total number of cells

    //Diagonal coefficients
    for(label i=0; i<Nc; i++)
    {
        Info<< UEqn.diag()[i] << endl;
    }

    //Off-diagonal coefficients
    for(label faceI=0; faceI<UEqn.lduAddr().lowerAddr().size(); faceI++)
    {
        Info<< UEqn.upper()[faceI] << endl; //Lower coefficients
        Info<< UEqn.lower()[faceI] << endl; //Upper coefficients
    }

Hope that helps,
Fumiya

[hawkeye321]

Fumiya Thanks for the answer

Would you please also let me know what does "lduAddr()" and "lowerAddr()" mean here?

[fumiya]

Hi mahdiiowa,

I don't understand what exactly "lduAddr()" means,
but you can access the "owner" cell label by UEqn.lduAddr().lowerAddr()[faceI]
and the "neighbour" cell label by UEqn.lduAddr().upperAddr()[faceI], respectively.

Hope that helps,
Fumiya

[Yahoo]

Thanks for the comment.

[EmadTandis]

Originally Posted by mahdiiowa

Hi FOAMERS

I am trying to access to elements of the coefficient matrix of an equation using lduAddresing

Foloowing the procedure in "OpenFOAM guide/Matrices in OpenFOAM" http://openfoamwiki.net/index.php/Matrices_in_OpenFOAM

I use

scalarField *lowerPtr_;
scalarField *diagPtr_;
scalarField *upperPtr_;

for (register label face=0; face<l.size(); face++)
{
Diag[l[face]] -= Lower[face];
Diag[u[face]] -= Upper[face];
}

But I get the following errors

error: ‘l’ was not declared in this scope
error: ‘Diag’ was not declared in this scope
error: ‘Lower’ was not declared in this scope
error: ‘u’ was not declared in this scope
error: ‘Upper’ was not declared in this scope

Do you have any idea what these errors are?
PS: I have included these header files

#include "fvCFD.H"
#include "wallFvPatch.H"
#include <fstream>
#include <string>
#include <stdio.h>
#include <time.h>
#include <iostream>

Thanks mahdiiowa for your good link
this error is because the Diag, upper, ... is defined as a function in special class (I think lduMatrix or fvMatrix). If you want to call this function, you have to construct object from those class. example:
ldumatrix A
... (Defenition A)
now, you can call A.diag() and A.upper() , ...
I Hope this help you.
best Regards
Emad Tandis

[santiagomarquezd]

Hi, you can use gdbOF to access the matrix coefficients at debugging time.

Regards.

[thang_cegroup]

We use lduaddressing to access non-zero coefficient matrix. But how about the zero coefficient ? Can we access them too ? I need to change some zero coefficient into non-zero coefficient to implement cyclic bc for block matrix in foam extend 4-0 since the cyclic is not implemented for block matrix ?

Do you guys have any suggest to solve this problem ? Even a few information might help a lot ?