为了进一步测试，建立一个3X3的网格，并创建一个fvMatrix。

        fvVectorMatrix UEqn_lap
        (
          fvm::laplacian(nu, U)
        );

网格编号如下：

|6|7|8|
|3|4|5|
|0|1|2|

输出 upper 和lower 的address，没有问题

        Info<< "UEqn_lap upper add before sum   " << UEqn_lap.lduAddr().upperAddr() << nl << endl;
        Info<< "UEqn_lap lower add before sum   " << UEqn_lap.lduAddr().lowerAddr() << nl << endl;
UEqn_lap upper add before sum
12
(
1
3
2
4
5
4
6
5
7
8
7
8
)


UEqn_lap lower add before sum
12
(
0
0
1
1
2
3
3
4
4
5
6
7
)

下面把upper lower 和diag重新赋值并输出

        UEqn_lap.upper() = 100;
        UEqn_lap.diag() = 0;
        UEqn_lap.lower() = 1;

        Info<< "UEqn_lap diag after change   " << UEqn_lap.diag() << nl << endl;
        Info<< "UEqn_lap upper after change   " << UEqn_lap.upper() << nl << endl;
        Info<< "UEqn_lap lower after change   " << UEqn_lap.lower() << nl << endl;

输出为：

UEqn_lap diag after change   9{0}

UEqn_lap upper after change   12{100}

UEqn_lap lower after change   12{1}

调用

        UEqn_lap.negSumDiag();

        Info<< "UEqn_lap diag after negSumDiag   " << UEqn_lap.diag() << nl << endl;

输出结果为：

UEqn_lap diag after negSumDiag   9(-2 -102 -101 -102 -202 -201 -101 -201 -200)

我特意把upper里面的数值设的很大，从输出的结果看，diag 是把lower里面的数加和了，也就是把列中的数加和了。

这个结果让我有点儿弄不明白。

CFD中文网

mengweilm425