#include "global.h"

void SetTemperature(v,N,m,T)
struct vector  *v;
int  N;
double T,m;
{  

  //
  // Subroutine sets the temperature of the atoms in the system
  // to T assuming a Maxwell-Boltzmann velocity distribution
  //


  // stuff removed

}


void GetTemperature(v,N,m,T)

struct vector  *v;
int  N;
double  *T,m;
{
  //
  // Subroutine gets the temperature of the atoms in the system
  //

  double  Ekinsum,gaussianrand();
  int  i;

  // Get sum of kinetic energies
  Ekinsum=0;
  for (i=1;i<=N;i++) {

     Ekinsum=Ekinsum+v[i].x*v[i].x+v[i].y*v[i].y+v[i].z*v[i].z;

  }
  // Do necessary unit transformations to get E in J
  Ekinsum=0.5*m*u*Ekinsum*vunit*vunit;
  // and get temperature using E=3/2 N k T
  *T = Ekinsum/(1.5*N*kB);
  
}

void GetEnergies(v,N,m,Ekinsum,Epotsum,Etot,Ekin,Epot)

struct vector  *v;
double  *Ekinsum,*Epotsum,*Etot,m,*Ekin,*Epot;
int  N;
{
  double  help1;
  int  i;
  
  *Ekinsum=0.0; *Epotsum=0.0;
  help1=0.5*m*u*vunit*vunit/e;
  for (i=1;i<=N;i++) {
     Ekin[i]=help1*(v[i].x*v[i].x+v[i].y*v[i].y+v[i].z*v[i].z);
     *Ekinsum = *Ekinsum+Ekin[i];
     *Epotsum = *Epotsum+Epot[i];
  }
  *Etot = *Ekinsum+*Epotsum;

}
  



//---------------------------------------------------------------------------

double uniformrand(seed)
long *seed;  

  // -----------------------------------------------------------
  // Park-Miller "minimal" Random number generator
  // uniform distribution ]0,1[ . See Numerical Recipes ch. 7.0
  // -----------------------------------------------------------
{
  long IA,IM,IQ,IR,MASK;
  double AM,ans;
  long  k;
  
  IA=16807; IM=2147483647; AM=1.0/IM; IQ=127773; IR=2836; MASK=123459876;

  *seed= (*seed)^MASK;
  k=(*seed)/IQ;
  *seed = IA*((*seed)-k*IQ)-IR*k;
  if (*seed < 0) *seed = *seed+IM;
  ans=AM*(*seed);
  *seed= (*seed)^MASK;
  
  return(ans);
}



double gaussianrand(seed)
int  *seed;

  // ---------------------------------------------------------
  // Random numbers with normal (Gaussian) distribution.
  // Mean is 0 and standard deviation is 1
  // ---------------------------------------------------------
{

  // stuff removed

}