pq.h

#pragma once
#include "ff/atom.h"
 
namespace tinker {
constexpr int BOUNDS_EVERY_X_STEPS = 500;
 
void mdPos(time_prec, pos_prec*, pos_prec*, pos_prec*, //
   const vel_prec*, const vel_prec*, const vel_prec*);
void mdPos(time_prec dt);
 
void mdPosAxbv(pos_prec a, pos_prec b);
 
void mdPosAxbvAn(pos_prec (*a)[3], pos_prec (*b)[3]);
 
void mdDebugPosNorm(pos_prec poseps, time_prec dt, //
   const vel_prec* vx, const vel_prec* vy, const vel_prec* vz);
 
void mdVel(time_prec dt, const grad_prec* grx, const grad_prec* gry,
   const grad_prec* grz);
 
void mdVel2(time_prec dt, const grad_prec* grx, const grad_prec* gry,
   const grad_prec* grz, time_prec dt2, const grad_prec* grx2,
   const grad_prec* gry2, const grad_prec* grz2);
 
void mdVelScale(vel_prec scal, int nelem, vel_prec* vx0, vel_prec* vy0,
   vel_prec* vz0);
 
void mdVelAvbf(int nrespa, vel_prec a, vel_prec b,                   //
   const grad_prec* gx1, const grad_prec* gy1, const grad_prec* gz1, //
   const grad_prec* gx2, const grad_prec* gy2, const grad_prec* gz2);
 
void mdVelAvbfAn(int nrespa, vel_prec a[3][3], vel_prec b[3][3],     //
   const grad_prec* gx1, const grad_prec* gy1, const grad_prec* gz1, //
   const grad_prec* gx2, const grad_prec* gy2, const grad_prec* gz2);
 
void mdVelData(RcOp);
}
 
//====================================================================//
//                                                                    //
//                          Global Variables                          //
//                                                                    //
//====================================================================//
 
namespace tinker {
TINKER_EXTERN vel_prec *vx, *vy, *vz;
}