DTD_BConds< T, 3, PPNX > Struct Template Reference

specialization for a slab, general cell More...

Collaboration diagram for DTD_BConds< T, 3, PPNX >:

Public Member Functions
	DTD_BConds (const CrystalLattice< T, 3 > &lat)
T	get_min_distance (TinyVector< T, 3 > &a) const
	evaluate the minimum distance More...
T	apply_bc (TinyVector< T, 3 > &displ) const
	apply BC to a displacement vector a and return the minimum-image distance More...
void	apply_bc (std::vector< TinyVector< T, 3 >> &dr, std::vector< T > &r, std::vector< T > &rinv) const
void	apply_bc (std::vector< TinyVector< T, 3 >> &dr, std::vector< T > &r) const
void	evaluate_rsquared (TinyVector< T, 3 > *restrict dr, T *restrict rr, int n)

Public Attributes
T	r00
T	r10
T	r01
T	r11
T	g00
T	g10
T	g01
T	g11
T	r2max
std::vector< TinyVector< T, 3 > >	nextcells

Detailed Description

template<class T>
struct qmcplusplus::DTD_BConds< T, 3, PPNX >

specialization for a slab, general cell

Definition at line 638 of file ParticleBConds3D.h.

Constructor & Destructor Documentation

DTD_BConds()

DTD_BConds ( const CrystalLattice< T, 3 > &  lat )

inline

Definition at line 645 of file ParticleBConds3D.h.

      : r00(lat.R(0)),
        r10(lat.R(3)),
        r01(lat.R(1)),
        r11(lat.R(4)),
        g00(lat.G(0)),
        g10(lat.G(3)),
        g01(lat.G(1)),
        g11(lat.G(4)),
        r2max(lat.CellRadiusSq)
  {
    nextcells.resize(8);
    int ic = 0;
    for (int i = -1; i <= 1; ++i)
      for (int j = -1; j <= 1; ++j)
      {
        if (!(i || j))
          continue; //exclude zero
        nextcells[ic][0] = i * r00 + j * r10;
        nextcells[ic][1] = i * r01 + j * r11;
        nextcells[ic][2] = 0;
        ++ic;
      }
  }

Member Function Documentation

apply_bc() [1/3]

T apply_bc ( TinyVector< T, 3 > &  displ ) const

inline

apply BC to a displacement vector a and return the minimum-image distance

Parameters

lat	lattice
a	displacement vector

Returns

the minimum-image distance

Definition at line 708 of file ParticleBConds3D.h.

  {
    //cart2unit
    TinyVector<T, 2> ar(displ[0] * g00 + displ[1] * g10, displ[0] * g01 + displ[1] * g11);
    //put them in the box
    ar[0] -= round(ar[0]);
    ar[1] -= round(ar[1]);
    //unit2cart
    displ[0] = ar[0] * r00 + ar[1] * r10;
    displ[1] = ar[0] * r01 + ar[1] * r11;
    //return |displ|^2 after checking the ghost cells
    return get_min_distance(displ);
  }

apply_bc() [2/3]

void apply_bc ( std::vector< TinyVector< T, 3 >> &  dr, std::vector< T > &  r, std::vector< T > &  rinv  ) const

inline

Definition at line 722 of file ParticleBConds3D.h.

References DTD_BConds< T, D, SC >::apply_bc(), qmcplusplus::simd::inv(), qmcplusplus::n, and qmcplusplus::simd::sqrt().

  {
    const int n = dr.size();
    for (int i = 0; i < n; ++i)
      rinv[i] = apply_bc(dr[i]);
    simd::sqrt(&rinv[0], &r[0], n);
    simd::inv(&r[0], &rinv[0], n);
  }

apply_bc() [3/3]

void apply_bc ( std::vector< TinyVector< T, 3 >> &  dr, std::vector< T > &  r  ) const

inline

Definition at line 731 of file ParticleBConds3D.h.

References DTD_BConds< T, D, SC >::apply_bc().

  {
    for (int i = 0; i < dr.size(); ++i)
      r[i] = apply_bc(dr[i]);
  }

evaluate_rsquared()

void evaluate_rsquared ( TinyVector< T, 3 > *restrict  dr, T *restrict  rr, int  n  )

inline

Definition at line 737 of file ParticleBConds3D.h.

References DTD_BConds< T, D, SC >::apply_bc(), and qmcplusplus::n.

  {
    for (int i = 0; i < n; ++i)
      rr[i] = apply_bc(dr[i]);
  }

get_min_distance()

T get_min_distance ( TinyVector< T, 3 > &  a ) const

inline

evaluate the minimum distance

Parameters

lat	lattice
a	displacement vector
r2max	square of the maximum cutoff

Returns

square of the minimum-image distance

Search the ghost cells to match Wigner-Seitz cell

Definition at line 678 of file ParticleBConds3D.h.

  {
    T d2 = a[0] * a[0] + a[1] * a[1] + a[2] * a[2];
    if (d2 < r2max)
      return d2;
    else
    {
      T d2min = d2;
      int ic  = -1;
      for (int i = 0; i < 8; ++i)
      {
        TinyVector<T, 3> c(a + nextcells[i]);
        d2 = c[0] * c[0] + c[1] * c[1] + c[2] * c[2];
        if (d2 < d2min)
        {
          d2min = d2;
          ic    = i;
        }
      }
      if (ic >= 0)
        a += nextcells[ic];
      return d2min;
    }
  }

Member Data Documentation

g00

T g00

Definition at line 641 of file ParticleBConds3D.h.

g01

T g01

Definition at line 641 of file ParticleBConds3D.h.

g10

T g10

Definition at line 641 of file ParticleBConds3D.h.

g11

T g11

Definition at line 641 of file ParticleBConds3D.h.

nextcells

std::vector<TinyVector<T, 3> > nextcells

Definition at line 643 of file ParticleBConds3D.h.

r00

T r00

Definition at line 640 of file ParticleBConds3D.h.

r01

T r01

Definition at line 640 of file ParticleBConds3D.h.

r10

T r10

Definition at line 640 of file ParticleBConds3D.h.

r11

T r11

Definition at line 640 of file ParticleBConds3D.h.

r2max

T r2max

Definition at line 642 of file ParticleBConds3D.h.

---

The documentation for this struct was generated from the following file:

• /home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/Particle/Lattice/ParticleBConds3D.h