[pablodecastillo]

hello,

I found that findMax not working in parallel with scalarField.

I want to find the min coordinate_X for my pacht and at the same time with coordinate_Z <= 0.

Any idea?

Thanks

[Phicau]

I've never tried to use findMax, but this should work:

scalar whatever = findMax...;

reduce(whatever, maxOp<scalar>());

There are other options depending on your needs: minOp, sumOp...

Regards

[pablodecastillo]

So it would be;

label whatever = findMax...;

reduce(whatever, maxOp<label>());

I am using gMax, so get the max value, but now how can i get what label is the max value (and working in parallel)

Thanks

[Phicau]

Just loop the entire patch:

    vectorField coord = patch().Cf();
    scalar auxiliar = 999999.0;

    forAll(coord.component(2), index) 
    {
        if(coord.component(2)[index]<=0.0)
        {
            if(coord.component(0)[index]<auxiliar)
            {
                auxiliar = coord.component(0)[index];
            } 
        } 
    }

reduce(auxiliar, minOp<scalar>());

[pablodecastillo]

thank you, but still i can not get the label for the value, any idea?

Pablo

[Phicau]

Another loop...

    label auxiliarLabel = -999999;

    forAll(coord.component(2), index2) 
    {
             if(coord.component(2)[index2]<=0.0)
             {
                      if(coord.component(0)[index2]==auxiliar)
                      {
                               auxiliarLabel = index2;
                               break;
                      } 
             } 
    }

reduce(auxiliarLabel, maxOp<label>());

[pablodecastillo]

the second loop to find the cell is not working; i mean it is not giving the correct label

label celdabuscada = -99999;
scalar pto_X = -10000.0;

forAll(mesh.boundary()[patchK], faceI)
{
if (NN2[faceI].component(vector::Z) <= 0.0) {
if (NN2[faceI].component(vector::X) >= pto_X){
pto_X=NN2[faceI].component(vector::X);
}
}
}
reduce(pto_X, maxOp<scalar>());
forAll(mesh.boundary()[patchK], face2)
{
if (NN2[face2].component(vector::Z) <= 0.0) {
if (NN2[face2].component(vector::X)==pto_X){
celdabuscada=face2;
Info <<"////celda buscada " << celdabuscada << nl;
break;
}
}
}
reduce(celdabuscada, maxOp<label>());

more ideas??