[bieshuxuhe]

hi!
look down the piece of code of "pEqn.H" in interFoam (openfoam 2.2.0):

00049     #include "continuityErrs.H"
00050 
00051     p == p_rgh + rho*gh;
00052 
00053     if (p_rgh.needReference())
00054     {
00055         p += dimensionedScalar
00056         (
00057             "p",
00058             p.dimensions(),
00059             pRefValue - getRefCellValue(p, pRefCell)
00060         );
00061         p_rgh = p - rho*gh;
00062     }
00063 }

the red code "p == p_rgh + rho*gh;" use the operator "==" , and the green code "p_rgh = p - rho*gh;" use the operator "=" !
Then I look up the source guide .
I could only understand from the source guide that the difference between "=" and "==" is that:
operator "=" will check first to make sure the variable couldn't be assigned to itself , while "==" won't check that ! ;
I wonder is there any other difference between "==" and "=" ?

could you help me ?
thanks!

[kmooney]

I "think" that the == will force an overwrite of the current boundary values even if they are set to a non-Dirichlet type. This is done in addition to altering the internalField values. I recall seeing a discussion on this somewhere in this forum but I can't seem to find it.

I might be wrong about this so don't hold me accountable!

[bieshuxuhe]

thanks for your reply!

[bieshuxuhe]

it is there:

http://www.cfd-online.com/Forums/openfoam-programming-development/85013-operator.html

Nope: operator== will FORCE assignment, even if patch field type says it should not be so. For example for a fixedValue patch, operator== will change its value and operator= will do nothing.

Hrv
__________________
Hrvoje Jasak