| 1 | 
tim | 
1423 | 
#include "BaseLattice.hpp" | 
| 2 | 
  | 
  | 
void BaseLattice::getLatticePoints(vector<double>& latticePosX, vector<double>& latticePosY, vector<double>& latticePosZ,  | 
| 3 | 
  | 
  | 
                                                           int nx, int ny, int nz){ | 
| 4 | 
  | 
  | 
                                                            | 
| 5 | 
  | 
  | 
  for( int i=0;i < nCellSites;i++){ | 
| 6 | 
  | 
  | 
   | 
| 7 | 
  | 
  | 
    latticePosX[i] = startX + thisUnitCell.sx[i] + cellLengthX * (double( ix ) - 0.5); | 
| 8 | 
  | 
  | 
    latticePosY[i] = startY + thisUnitCell.sy[i] + cellLengthY * (double( iy ) - 0.5); | 
| 9 | 
  | 
  | 
    latticePosZ[i] = startZ + thisUnitCell.sz[i] + cellLengthZ * (double( iz ) - 0.5); | 
| 10 | 
  | 
  | 
     | 
| 11 | 
  | 
  | 
  } | 
| 12 | 
  | 
  | 
 | 
| 13 | 
  | 
  | 
} | 
| 14 | 
  | 
  | 
 |