[jaming]

Hi Foamers,
I have a problem setting in the forAll() loop the right condition:
##########################
forAll(alpha1,celli)
{
if (alpha1[celli] > 1 )
{

Pout<<"alpha1= "<<alpha1[celli]<<endl;
}
}
##########################

I also get

alpha1= 1

in the results what I think is not logic.

So if someone has any idea I will be thankful
Yamine

[ngj]

Hallo,

It is probably related to your write precision in your controlDict, e.g. if the write precision is 6, then 1.000000007 is written as 1.

You could test this by using the following code:

Pout << alpha1[celli] << tab << 1 - alpha1[celli] << endl;

Kind regards,

Niels

[jaming]

Many thanks for you quick replay,
but look for example hier:

forAll(alpha1, celli)
    {
     if ( alpha1[celli] < scalar(1.))
       {
    
         Pout<< "Alpha1=  "<<alpha1[celli]<< endl;
       }

I get the flow:

PHP Code:

...........   ..................
Alpha1=  0.0144193383
Alpha1=  0.9575155878
Alpha1=  1
Alpha1=  0.1574193345
Alpha1=  5.502415325e-11
Alpha1=  5.162033835e-11
Alpha1=  8.173407791e-13
...........    ....................... 


The idea is to set some conditions with alpha1 but I found that alpha1[celli] didn't pass the wished value at all.
And like you can see unfortunately 1 is not filtered here

With kind regards,
Yamine

[ngj]

With your recent code, you turned around the inequality, and the value being 1, can actually be 0.9999999999, i.e. to me the code still does the rigth thing.

Kind regards,

Niels

[jaming]

Hi ngj

Many Thanks for your unpayable help.

Regards
Yamine