int Atom::unpack_border(int i, double *buf)
{
if (i == nmax) growarray();
int m = 0;
x[i][0] = buf[m++];
x[i][1] = buf[m++];
x[i][2] = buf[m++];
return m;
}
int Atom::unpack_border(MMD_int i, MMD_float* buf)
{
if(i == nmax) growarray();
int m = 0;
x[i * PAD + 0] = buf[m++];
x[i * PAD + 1] = buf[m++];
x[i * PAD + 2] = buf[m++];
type[i] = buf[m++];
return m;
}
void Atom::addatom(double x_in, double y_in, double z_in,
double vx_in, double vy_in, double vz_in)
{
if (nlocal == nmax) growarray();
x[nlocal][0] = x_in;
x[nlocal][1] = y_in;
x[nlocal][2] = z_in;
v[nlocal][0] = vx_in;
v[nlocal][1] = vy_in;
v[nlocal][2] = vz_in;
nlocal++;
}
void Atom::addatom(MMD_float x_in, MMD_float y_in, MMD_float z_in,
MMD_float vx_in, MMD_float vy_in, MMD_float vz_in)
{
if(nlocal == nmax) growarray();
x[nlocal*PAD + 0] = x_in;
x[nlocal*PAD + 1] = y_in;
x[nlocal*PAD + 2] = z_in;
v[nlocal*PAD + 0] = vx_in;
v[nlocal*PAD + 1] = vy_in;
v[nlocal*PAD + 2] = vz_in;
type[nlocal] = rand()%ntypes;
nlocal++;
}
