
[Sponsors] 
September 15, 2014, 12:52 
Vector Mathematics

#1 
New Member
Will
Join Date: Dec 2011
Posts: 15
Rep Power: 5 
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
Code:
const vector offset((endCorner_.x()  startCorner_.x())/(nPoints_.x()  1), (endCorner_.y()  startCorner_.y())/(nPoints_.y()  1), (endCorner_.z()  startCorner_.z())/(nPoints_.z()  1) ); 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 Code:
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; } } } } 

September 15, 2014, 14:25 

#2 
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,641
Rep Power: 25 
Hallo,
I think that the only thing you lack is the Foam::cmptDivide functionality. The vector minus vector and vector minus scalar does already work. Have fun playing with waves2Foam; curious to know how the Lagrangian particles will behave across the free surface. Kind regards, Niels
__________________
Please note that I do not use the Friendfeature, so do not be offended, if I do not accept a request. 

September 16, 2014, 06:57 

#3  
New Member
Will
Join Date: Dec 2011
Posts: 15
Rep Power: 5 
Thanks.
Just for the record, I replace: Code:
const vector offset((endCorner_.x()  startCorner_.x())/(nPoints_.x()  1), (endCorner_.y()  startCorner_.y())/(nPoints_.y()  1), (endCorner_.z()  startCorner_.z())/(nPoints_.z()  1) ); Code:
const vector offset(Foam::cmptDivide(endCorner_  startCorner_, nPoints_  vector(1,1,1))); Code:
vector pos(startCorner_.x() + offset.x()*i, startCorner_.y() + offset.y()*j, startCorner_.z() + offset.z()*k ); Code:
vector pos(startCorner_ + Foam::cmptMultiply(offset, vector(i,j,k))); Quote:
As for WaveFoam. Preliminary results look amazing and I think our physical understanding of what we're doing will be greatly enhanced. We're loosing a few particles to the atmosphere, mostly due to wave splashing onto a surface, but that aside it seems acceptable. Will post some pictures when we're happier. 

Tags 
vector mathematics 
Thread Tools  
Display Modes  

