Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F62442151
lal_tersoff_mod.h
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Mon, May 13, 05:43
Size
3 KB
Mime Type
text/x-c++
Expires
Wed, May 15, 05:43 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
17647195
Attached To
rLAMMPS lammps
lal_tersoff_mod.h
View Options
/***************************************************************************
tersoff_mod.h
-------------------
Trung Dac Nguyen
Class for acceleration of the tersoff/mod pair style.
__________________________________________________________________________
This file is part of the LAMMPS Accelerator Library (LAMMPS_AL)
__________________________________________________________________________
begin :
email : ndactrung@gmail.com
***************************************************************************/
#ifndef LAL_TERSOFF_MOD_H
#define LAL_TERSOFF_MOD_H
#include "lal_base_three.h"
namespace
LAMMPS_AL
{
template
<
class
numtyp
,
class
acctyp
>
class
TersoffMod
:
public
BaseThree
<
numtyp
,
acctyp
>
{
public:
TersoffMod
();
~
TersoffMod
();
/// Clear any previous data and set up for a new LAMMPS run for generic systems
/** \param max_nbors initial number of rows in the neighbor matrix
* \param cell_size cutoff + skin
* \param gpu_split fraction of particles handled by device
*
* Returns:
* - 0 if successfull
* - -1 if fix gpu not found
* - -3 if there is an out of memory error
* - -4 if the GPU library was not compiled for GPU
* - -5 Double precision is not supported on card **/
int
init
(
const
int
ntypes
,
const
int
nlocal
,
const
int
nall
,
const
int
max_nbors
,
const
double
cell_size
,
const
double
gpu_split
,
FILE
*
screen
,
int
*
host_map
,
const
int
nelements
,
int
***
host_elem2param
,
const
int
nparams
,
const
double
*
lam1
,
const
double
*
lam2
,
const
double
*
lam3
,
const
double
*
powermint
,
const
double
*
biga
,
const
double
*
bigb
,
const
double
*
bigr
,
const
double
*
bigd
,
const
double
*
c1
,
const
double
*
c2
,
const
double
*
c3
,
const
double
*
c4
,
const
double
*
c5
,
const
double
*
h
,
const
double
*
beta
,
const
double
*
powern
,
const
double
*
powern_del
,
const
double
*
ca1
,
const
double
*
cutsq
);
/// Clear all host and device data
/** \note This is called at the beginning of the init() routine **/
void
clear
();
/// Returns memory usage on device per atom
int
bytes_per_atom
(
const
int
max_nbors
)
const
;
/// Total host memory used by library for pair style
double
host_memory_usage
()
const
;
// --------------------------- TYPE DATA --------------------------
/// If atom type constants fit in shared memory, use fast kernels
bool
shared_types
;
/// Number of atom types
int
_lj_types
;
/// ts1.x = lam1, ts1.y = lam2, ts1.z = lam3, ts1.w = powermint
UCL_D_Vec
<
numtyp4
>
ts1
;
/// ts2.x = biga, ts2.y = bigb, ts2.z = bigr, ts2.w = bigd
UCL_D_Vec
<
numtyp4
>
ts2
;
/// ts3.x = beta, ts3.y = powern, ts3.z = powern_del, ts3.w = ca1
UCL_D_Vec
<
numtyp4
>
ts3
;
/// ts4.x = c1, ts4.y = c2, ts4.z = c3, ts4.w = c4
UCL_D_Vec
<
numtyp4
>
ts4
;
/// ts5.x = c5, ts5.y = h
UCL_D_Vec
<
numtyp4
>
ts5
;
UCL_D_Vec
<
numtyp
>
cutsq
;
UCL_D_Vec
<
int
>
elem2param
;
UCL_D_Vec
<
int
>
map
;
int
_nparams
,
_nelements
;
/// Per-atom arrays:
/// zetaij.x = force, zetaij.y = prefactor, zetaij.z = evdwl,
/// zetaij.w = zetaij
UCL_D_Vec
<
acctyp4
>
_zetaij
;
UCL_Kernel
k_zeta
;
UCL_Texture
ts1_tex
,
ts2_tex
,
ts3_tex
,
ts4_tex
,
ts5_tex
;
numtyp
_cutshortsq
;
private:
bool
_allocated
;
void
loop
(
const
bool
_eflag
,
const
bool
_vflag
,
const
int
evatom
);
};
}
#endif
Event Timeline
Log In to Comment