Re: velocity component

[robin.trunk]

Hi,

what is your goal? Do you want to know the velocity at a specific point? In this case the code below should work:

SuperLatticePhysVelocity3D<T, DESCRIPTOR> velocity(sLattice, converter);
AnalyticalFfromSuperF3D<T> intpolateVelocity( velocity, true );
T v1[3];
intpolateVelocity( v1,point1 );

Here point1 is an array of T, for example
T point1[3] = {1., 0., 0.};

The code for 2D should be analogue.

Best
Robin