[romant]

I have a question. Let's say I have a field

const volScalarField& k

and I have a some labels, let's call them

label faceI        // surface between two cells
label faceCellI   // cell center

when I now call

k[faceCellI]

I get the value of k at the cell center, as the label points to a cell center. When I call

k[faceI]

do I get the interpolated value from the face? I always thought that I would need to interpolate the field first and then call the interpolated field in order to get my value at the interface between two cells. However, after implementing it this way, I still get a value, which also seems reasonable.

Now the question, does OpenFOAM automatically interpolate in these cases? I wonder if the same also works for volVectorFields, when retrieving a value from the interface between two neighboring cells (through a label) that it then interpolates.