[David*]

Hi Foamers,

I have some problems with a vector assignment. To do some further computations, I need the unit normal vector of a patch stored as a vector. My command for this is

vector unitNormalVector;
...
forAll(wallshearstress4Rho.boundaryField(), patchi)
{
unitNormalVector(-mesh.Sf().boundaryField()[patchi]/mesh.magSf().boundaryField()[patchi]);
...
}

I tried also with an equal (=) symbol. I assume, that on both sides of the equation is a vector with three components. But on the right hand side, it's something more complicated, as wmake's error message shows:

error: no match for call to ‘(Foam::vector) (Foam::tmp<Foam::Field<Foam::Vector<double> > >)’

Which looks like a template vector field (?). But how can there be a field, if I choose a specific patch by using [patchi] within a forAll-loop?
I'm a beginner with C++, so I ask someone more experienced to tell me what I'm doing wrong.

Regards,
David