integrate_kernel.h

/* ----------------------------------------------------------------------
   miniMD is a simple, parallel molecular dynamics (MD) code.   miniMD is
   an MD microapplication in the Mantevo project at Sandia National 
   Laboratories ( http://www.mantevo.org ). The primary 
   authors of miniMD are Steve Plimpton (sjplimp@sandia.gov) , Paul Crozier 
   (pscrozi@sandia.gov) and Christian Trott (crtrott@sandia.gov).

   Copyright (2008) Sandia Corporation.  Under the terms of Contract
   DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
   certain rights in this software.  This library is free software; you 
   can redistribute it and/or modify it under the terms of the GNU Lesser 
   General Public License as published by the Free Software Foundation; 
   either version 3 of the License, or (at your option) any later 
   version.
  
   This library is distributed in the hope that it will be useful, but
   WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
   Lesser General Public License for more details.
    
   You should have received a copy of the GNU Lesser General Public 
   License along with this software; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
   USA.  See also: http://www.gnu.org/licenses/lgpl.txt .

   For questions, contact Paul S. Crozier (pscrozi@sandia.gov) or
   Christian Trott (crtrott@sandia.gov). 

   Please read the accompanying README and LICENSE files.
---------------------------------------------------------------------- */

//#define MMD_floatK float
//#define MMD_floatK3 float3
//#define MMD_floatK3 float3;
//#define MMD_float float;

#include "precision.h"

__kernel void integrate_initial(__global MMD_floatK3* x, __global MMD_floatK3* v,__global MMD_floatK3* f, int nlocal, MMD_float dt,MMD_float dtforce,int nmax)
{
  int i = get_global_id(0);
  if(i<nlocal)
  {
	  v[i] += dtforce*f[i];
	  x[i] += dt*v[i];
  }
}

__kernel void integrate_final(__global MMD_floatK3* v, __global MMD_floatK3* f, int nlocal, MMD_float dt,int nmax)
{
  int i = get_global_id(0);
  if(i<nlocal)
  {
	  v[i] += dt*f[i];
  }

}