[cheng1988sjtu]

Hi Foamers,

I recently compiled OpenFOAM-2.1.0 with gcc version 4.8.2, I can compile OpenFOAM without error messages, except that I keep getting the warning:

In file included from lnInclude/Field.H:360:0,
from lnInclude/labelField.H:39,
from lnInclude/primitiveFields.H:37,
from lnInclude/pointField.H:36,
from lnInclude/edge.H:40,
from lnInclude/edgeList.H:32,
from lnInclude/primitiveMesh.H:57,
from lnInclude/polyMesh.H:44,
from lnInclude/regIOobject.C:28,
from Make/linux64GccDPOpt/global.C:71:
lnInclude/Field.C: In member function 'void Foam::Field<Type>:perator=(const Foam::VectorSpace<Form, Cmpt, nCmpt>&)':
lnInclude/Field.C:680:42: warning: typedef 'VSType' locally defined but not used [-Wunused-local-typedefs]
typedef VectorSpace<Form,Cmpt,nCmpt> VSType;
^

I've never encountered this warning message when I compile OpenFOAM-2.1.0 with gcc-4.6.1, but this time I have to use gcc-4.8.2, because it's the only available version. Anyway, I compiled OpenFOAM-2.1.0 with gcc-4.8.2 by ignoring the warning message.

However, When I want to identify certain boundary type such as symmetryFvPatchScalarField, I try the following lines, which used to work on the OF2.1.0 compiled with gcc4.6.1:

forAll(p.boundaryField(), patchi)
{
if (isA<symmetryFvPatchScalarField>(p.boundaryField()[patchi]))
{
Info<<"symmetry correct"<<endl;
phiDraga.boundaryField()[patchi] = 0.0;
phiDragb.boundaryField()[patchi] = 0.0;
}
}

This time, it doesn't do anything ( don't print out "symmetry correct") even if I define p boundary as "symmetryPlane", and in blockMesh I've defined corresponding "symmetryPlane" boundary.

Do you have any ideas about how to solve this problem without changing gcc version? Because I don't have the permission to compile my own gcc on that cluster.

Any hints will be highly appreciated, Thanks!

[chrisb2244]

Does

forAll(p.boundaryField(), patchI)
{
    if (p.boundaryField().types()[patchI] == symmetryPolyPatch)
    {
        // Do stuff
    }
}

work?

[cheng1988sjtu]

Thanks, Christian, your method works, and you saved my life!!!

    forAll(p.boundaryField(), patchi)
    {
//        if (isA<symmetryFvPatchScalarField>(p.boundaryField()[patchi]))
        if (p.boundaryField().types()[patchi] == "symmetryPlane")
        {
//            Info<<p.boundaryField().types()[patchi]<<endl;
            phiDraga.boundaryField()[patchi] = 0.0;
            phiDragb.boundaryField()[patchi] = 0.0;
        }
    }