/home/eoseret/qaas/qaas_runs/178-547-6790/intel/miniqmc/build/miniqmc/src/Numerics/OhmmsPETE/VectorSoAContainer.h: 243 - 271
--------------------------------------------------------------------------------

243:       for (size_t i = 0; i < D; ++i)
244:         *(_base + M * i) = rhs[i];
[...]
271:   inline T* restrict data(size_t i) { return myData + i * nGhosts; }
/usr/include/c++/14/bits/stl_vector.h: 1128 - 1263
--------------------------------------------------------------------------------

1128:       operator[](size_type __n) _GLIBCXX_NOEXCEPT
1129:       {
1130: 	__glibcxx_requires_subscript(__n);
1131: 	return *(this->_M_impl._M_start + __n);
[...]
1150: 	return *(this->_M_impl._M_start + __n);
[...]
1263:       { return _M_data_ptr(this->_M_impl._M_start); }
/home/eoseret/qaas/qaas_runs/178-547-6790/intel/miniqmc/build/miniqmc/src/Numerics/OhmmsPETE/TinyVector.h: 61 - 62
--------------------------------------------------------------------------------

61:     for (size_t d = 0; d < D; ++d)
62:       X[d] = T(0);
/home/eoseret/qaas/qaas_runs/178-547-6790/intel/miniqmc/build/miniqmc/src/QMCWaveFunctions/Jastrow/TwoBodyJastrowRef.h: 310 - 355
--------------------------------------------------------------------------------

310: void TwoBodyJastrowRef<FT>::acceptMove(ParticleSet& P, int iat)
311: {
312:   // get the old u, du, d2u
313:   const auto& d_table = P.getDistTableAA(myTableID);
314:   computeU3(P, iat, d_table.getDistRow(iat).data(), old_u.data(), old_du.data(), old_d2u.data());
315:   if (UpdateMode == ORB_PBYP_RATIO)
316:   { // ratio-only during the move; need to compute derivatives
317:     computeU3(P, iat, d_table.getTempDists().data(), cur_u.data(), cur_du.data(), cur_d2u.data());
318:   }
319: 
320:   valT cur_d2Uat(0);
321:   const auto& new_dr    = d_table.getTempDispls();
322:   const auto& old_dr    = d_table.getDisplRow(iat);
323:   constexpr valT lapfac = OHMMS_DIM - RealType(1);
324:   for (int jat = 0; jat < N; jat++)
325:   {
326:     const valT du   = cur_u[jat] - old_u[jat];
327:     const valT newl = cur_d2u[jat] + lapfac * cur_du[jat];
328:     const valT dl   = old_d2u[jat] + lapfac * old_du[jat] - newl;
329:     Uat[jat] += du;
330:     d2Uat[jat] += dl;
331:     cur_d2Uat -= newl;
332:   }
333:   posT cur_dUat;
334:   for (int idim = 0; idim < OHMMS_DIM; ++idim)
[...]
341:     valT cur_g                     = cur_dUat[idim];
342:     for (int jat = 0; jat < N; jat++)
343:     {
344:       const valT newg = cur_du_pt[jat] * new_dX[jat];
345:       const valT dg   = newg - old_du_pt[jat] * old_dX[jat];
346:       save_g[jat] -= dg;
347:       cur_g += newg;
348:     }
349:     cur_dUat[idim] = cur_g;
350:   }
351:   LogValue += Uat[iat] - cur_Uat;
352:   Uat[iat]   = cur_Uat;
353:   dUat(iat)  = cur_dUat;
354:   d2Uat[iat] = cur_d2Uat;
355: }
/home/eoseret/qaas/qaas_runs/178-547-6790/intel/miniqmc/build/miniqmc/src/Numerics/OhmmsPETE/OhmmsVector.h: 223 - 249
--------------------------------------------------------------------------------

223:     return X[i];
[...]
249:   inline const_pointer data() const { return X; }
