[ooo]

Hi,

Interpolating a volVectorField to a vector, works fine in serial, but not in parallel.
Do you know what would be the solution?

dictionary interpolationDict = mesh.solutionDict().subDict("interpolationSchemes");
autoPtr<interpolation<vector> > Uinterp = interpolation<vector>::New(interpolationDict, U);
forAll (vector1, i)
vector1[i] = Uinterp->interpolate(coordinates[i], cellI) ;

Tnx

[tomislav_maric]

Don't roll your own for loops for field operations that are already implemented in OF. 99% of the field operations, such as the surface interpolations are already parallelized.

Use the member function:

virtual tmp< GeometricField
 < Type, fvsPatchField, 
 surfaceMesh > > interpolate (const GeometricField< Type, fvPatchField, volMesh > &) const

to interpolate from one field to another. If you interpolate a volVectorField to a volVector field, it will work fine in parallel as well.

[ooo]

Originally Posted by tomislav_maric

Don't roll your own for loops for field operations that are already implemented in OF. 99% of the field operations, such as the surface interpolations are already parallelized.

Use the member function:

virtual tmp< GeometricField
 < Type, fvsPatchField, 
 surfaceMesh > > interpolate (const GeometricField< Type, fvPatchField, volMesh > &) const

to interpolate from one field to another. If you interpolate a volVectorField to a volVector field, it will work fine in parallel as well.

Thank you for your reply.
I didn't get your point exactly, because i'm using the interpolation function of OF which is supposed to be parallelized...
If you mean merging the third and fourth line of my code to :

vector1 = Uinterp->interpolate(coordinates,mesh.findCell(coordinates));

it gives my the compilation error of "no matching function..."

[tomislav_maric]

Np. As allways the question is what you actually want to achieve.

My answer stated that if you want to interpolate a vector field onto a vector (it was not clear if your 'vector' is a position in space to me), you should do the interpolation from one field to the other like this:

volVectorField result = fvc::interpolate(input, "schemeName");

in which case you are using the member function I mentioned before.

However, if you want to interpolate a field onto a point that is located somewhere in the cell, use cell to point interpolation (see "interpolationCellPoint.H"). This gives you a field value at some point in space.

Both should be parallel automatically. In that you are right...

[cfd-bird]

I had the same problem,Did you solve the problem?