|
[Sponsors] |
November 25, 2009, 12:23 |
EigenVectors and volSymmTensorField
|
#1 |
New Member
Sachin Shanbhag
Join Date: Nov 2009
Posts: 3
Rep Power: 17 |
Hi, I am an OpenFOAM newbie using v1.6. I am modifying a solver I found on this forum and have a particular question, regarding the following code snippet.
volSymmTensorField A2 ( IOobject ( "A2", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); // // Investigate why I cant say just eigenVectors(A2) later. // volTensorField eigenVectorsA2 ( IOobject ( "eigenVectorsA2", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), // eigenVectors(A2) eigenVectors(A2 & UnitTensor) ); Where UnitTensor is (1 0 0 0 1 0 0 0 1). I have the following observations: 1. When I try to define the eigenVectors without "multiplying" by the UnitTensor ( as eigenVectors(A2) ), I get an error while compiling. 2. If I specify my A2 as volTensorField instead of volSymmTensorField, then I don't have the problem while compiling, even when I don't multiply by the UnitTensor. What is going on? I'll be happy to furnish compiler errors, if required. Sachin |
|
|
|