[willzyba]

I'm sure there must be a neater way to do this, but I can't find the documentation and nothing I've tried works. This is just ugly

const vector offset((endCorner_.x() - startCorner_.x())/(nPoints_.x() - 1),
                        (endCorner_.y() - startCorner_.y())/(nPoints_.y() - 1),
                        (endCorner_.z() - startCorner_.z())/(nPoints_.z() - 1) );

where startCorner and endCorner are both vectors that define opposite corners of a box.
And nPoints is an integer vector, which simply contains the number of points I want in each direction.

I really just want to say:
offset = (endCorner - startCorner)/(nPoints -1)

And hopefully C++ would figure that as nPoints is a vector, then when I subtract 1, that should remain a vector. And the subtraction of two vectors should be vector; and dividing one against the other give me my offsets in each direction. But it refused to compile.

So the final offset vector should contain the steps I need to take in each principal direction in order to fill a box with a lot of points. Just finishing off my development of LPT btw, that requires points throughout the domain which we're running through waveFOAM.



Finally function we're adding to solidParticleCloud.C looks like

void Foam::solidParticleCloud::seed(const dimensionedVector& g)
{
    // http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2011/OF_kurs_LPT_120911.pdf

    const vector offset((endCorner_.x() - startCorner_.x())/(nPoints_.x() - 1),
                        (endCorner_.y() - startCorner_.y())/(nPoints_.y() - 1),
                        (endCorner_.z() - startCorner_.z())/(nPoints_.z() - 1) );
  
    int Nz= int(nPoints_.z());
    int Ny= int(nPoints_.y());
    int Nx= int(nPoints_.x());

    // Hope this produces a line of particle points.
    for(int k=0; k<Nz; k++)
    {
        for(int j=0; j<Ny; j++)
        {
            for(int i=0; i<Nx; i++)
            { 
                vector pos(startCorner_.x() + offset.x()*i,
                           startCorner_.y() + offset.y()*j,
                           startCorner_.z() + offset.z()*k );

                label cellI =-1;
                label tetFaceI =-1;
                label tetPtI = -1;
                mesh_.findCellFacePt(pos, cellI, tetFaceI, tetPtI);
                if(cellI>0)
                {
                    solidParticle* pt= new solidParticle(mesh_, pos ,cellI, tetFaceI, tetPtI, 0.01, vector::zero );
                    Cloud<solidParticle>::addParticle(pt);
                }
                pos = pos + offset;
            }
        }
    }
}