[kmooney]

Hello All,

In a recent attempt to implement a tensor based viscosity field I ran into a problem with taking a Laplacian involving two volTensorFields. To perform some more simplified testing I altered laplacianFoam with the following:

        for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
        {
            solve
            (
                fvm::ddt(tau) - fvm::laplacian(dTau, tau)
            );
        }

where dTau and tau are both volTensorFields. I get a segfault within the HashTable lib. Below is the stack trace I get from my development code which is a little more comprehensive than the one that laplacianFoam gives me but crashes at the same point.

#0  0x00007ffff76e8850 in Foam::HashTable<Foam::tmp<Foam::fv::laplacianScheme<Foam::Vector<double>, Foam::Tensor<double> > > (*)(Foam::fvMesh const&, Foam::Istream&), Foam::word, Foam::string::hash>::find (this=0x0, key=...) at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/HashTable.C:177
#1  0x00007ffff76de993 in Foam::fv::laplacianScheme<Foam::Vector<double>, Foam::Tensor<double> >::New (mesh=..., schemeData=...)
    at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/laplacianScheme.C:73
#2  0x00007ffff76d646b in Foam::fvm::laplacian<Foam::Vector<double>, Foam::Tensor<double> > (gamma=..., vf=..., name=...)
    at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/fvmLaplacian.C:220
#3  0x00007ffff76cf35d in Foam::fvm::laplacian<Foam::Vector<double>, Foam::Tensor<double> > (tgamma=..., vf=..., name=...)
    at /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/fvmLaplacian.C:233
#4  0x00007ffff76c99be in Foam::LPTT::divTau (this=0x4c9c430, U=...) at viscoelasticLaws/LPTT/LPTT.C:146
#5  0x00007ffff7790ad6 in Foam::multiMode::divTau(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const ()
   from /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/lib/linux64GccDPOpt/libviscoelasticTransportModels.so
#6  0x00007ffff76c6e89 in Foam::viscoelasticModel::divTau(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const ()
   from /home/kmooney/OpenFOAM/OpenFOAM-1.6-ext/lib/linux64GccDPOpt/libviscoelasticTransportModels.so
#7  0x000000000041bbf1 in main ()

I'm a tad confused because the code compiles just fine and the Programmer's Guide says that those fields are accepted by FOAM in both implicit and explicit laplacian schemes. Also note that this error occurs even if all of the tensor fields are zeros. Does anyone have any idea as to why this is happening? Thanks a lot for any insight you guys might have for me.