Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F78536563
pair_lj_charmm_coul_long_kokkos.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
Wed, Aug 21, 11:53
Size
7 KB
Mime Type
text/x-c
Expires
Fri, Aug 23, 11:53 (2 d)
Engine
blob
Format
Raw Data
Handle
20061309
Attached To
rLAMMPS lammps
pair_lj_charmm_coul_long_kokkos.h
View Options
/* -*- c++ -*- ----------------------------------------------------------
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
http://lammps.sandia.gov, Sandia National Laboratories
Steve Plimpton, sjplimp@sandia.gov
Copyright (2003) Sandia Corporation. Under the terms of Contract
DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
certain rights in this software. This software is distributed under
the GNU General Public License.
See the README file in the top-level LAMMPS directory.
------------------------------------------------------------------------- */
#ifdef PAIR_CLASS
PairStyle
(
lj
/
charmm
/
coul
/
long
/
kk
,
PairLJCharmmCoulLongKokkos
<
LMPDeviceType
>
)
PairStyle
(
lj
/
charmm
/
coul
/
long
/
kk
/
device
,
PairLJCharmmCoulLongKokkos
<
LMPDeviceType
>
)
PairStyle
(
lj
/
charmm
/
coul
/
long
/
kk
/
host
,
PairLJCharmmCoulLongKokkos
<
LMPHostType
>
)
#else
#ifndef LMP_PAIR_LJ_CHARMM_COUL_LONG_KOKKOS_H
#define LMP_PAIR_LJ_CHARMM_COUL_LONG_KOKKOS_H
#include "pair_kokkos.h"
#include "pair_lj_charmm_coul_long.h"
#include "neigh_list_kokkos.h"
namespace
LAMMPS_NS
{
template
<
class
DeviceType
>
class
PairLJCharmmCoulLongKokkos
:
public
PairLJCharmmCoulLong
{
public:
enum
{
EnabledNeighFlags
=
FULL
|
HALFTHREAD
|
HALF
};
enum
{
COUL_FLAG
=
1
};
typedef
DeviceType
device_type
;
PairLJCharmmCoulLongKokkos
(
class
LAMMPS
*
);
~
PairLJCharmmCoulLongKokkos
();
void
compute
(
int
,
int
);
void
init_tables
(
double
cut_coul
,
double
*
cut_respa
);
void
init_style
();
double
init_one
(
int
,
int
);
struct
params_lj_coul
{
params_lj_coul
(){
cut_ljsq
=
0
;
cut_coulsq
=
0
;
lj1
=
0
;
lj2
=
0
;
lj3
=
0
;
lj4
=
0
;
offset
=
0
;};
params_lj_coul
(
int
i
){
cut_ljsq
=
0
;
cut_coulsq
=
0
;
lj1
=
0
;
lj2
=
0
;
lj3
=
0
;
lj4
=
0
;
offset
=
0
;};
F_FLOAT
cut_ljsq
,
cut_coulsq
,
lj1
,
lj2
,
lj3
,
lj4
,
offset
;
};
protected:
void
cleanup_copy
();
template
<
bool
STACKPARAMS
,
class
Specialisation
>
KOKKOS_INLINE_FUNCTION
F_FLOAT
compute_fpair
(
const
F_FLOAT
&
rsq
,
const
int
&
i
,
const
int
&
j
,
const
int
&
itype
,
const
int
&
jtype
)
const
;
template
<
bool
STACKPARAMS
,
class
Specialisation
>
KOKKOS_INLINE_FUNCTION
F_FLOAT
compute_fcoul
(
const
F_FLOAT
&
rsq
,
const
int
&
i
,
const
int
&
j
,
const
int
&
itype
,
const
int
&
jtype
,
const
F_FLOAT
&
factor_coul
,
const
F_FLOAT
&
qtmp
)
const
;
template
<
bool
STACKPARAMS
,
class
Specialisation
>
KOKKOS_INLINE_FUNCTION
F_FLOAT
compute_evdwl
(
const
F_FLOAT
&
rsq
,
const
int
&
i
,
const
int
&
j
,
const
int
&
itype
,
const
int
&
jtype
)
const
;
template
<
bool
STACKPARAMS
,
class
Specialisation
>
KOKKOS_INLINE_FUNCTION
F_FLOAT
compute_ecoul
(
const
F_FLOAT
&
rsq
,
const
int
&
i
,
const
int
&
j
,
const
int
&
itype
,
const
int
&
jtype
,
const
F_FLOAT
&
factor_coul
,
const
F_FLOAT
&
qtmp
)
const
;
Kokkos
::
DualView
<
params_lj_coul
**
,
Kokkos
::
LayoutRight
,
DeviceType
>
k_params
;
typename
Kokkos
::
DualView
<
params_lj_coul
**
,
Kokkos
::
LayoutRight
,
DeviceType
>::
t_dev_const_um
params
;
// hardwired to space for 15 atom types
params_lj_coul
m_params
[
MAX_TYPES_STACKPARAMS
+
1
][
MAX_TYPES_STACKPARAMS
+
1
];
F_FLOAT
m_cutsq
[
MAX_TYPES_STACKPARAMS
+
1
][
MAX_TYPES_STACKPARAMS
+
1
];
F_FLOAT
m_cut_ljsq
[
MAX_TYPES_STACKPARAMS
+
1
][
MAX_TYPES_STACKPARAMS
+
1
];
F_FLOAT
m_cut_coulsq
[
MAX_TYPES_STACKPARAMS
+
1
][
MAX_TYPES_STACKPARAMS
+
1
];
typename
ArrayTypes
<
DeviceType
>::
t_x_array_randomread
x
;
typename
ArrayTypes
<
DeviceType
>::
t_x_array
c_x
;
typename
ArrayTypes
<
DeviceType
>::
t_f_array
f
;
typename
ArrayTypes
<
DeviceType
>::
t_int_1d_randomread
type
;
typename
ArrayTypes
<
DeviceType
>::
t_float_1d_randomread
q
;
DAT
::
tdual_efloat_1d
k_eatom
;
DAT
::
tdual_virial_array
k_vatom
;
typename
ArrayTypes
<
DeviceType
>::
t_efloat_1d
d_eatom
;
typename
ArrayTypes
<
DeviceType
>::
t_virial_array
d_vatom
;
int
newton_pair
;
typename
ArrayTypes
<
DeviceType
>::
tdual_ffloat_2d
k_cutsq
;
typename
ArrayTypes
<
DeviceType
>::
t_ffloat_2d
d_cutsq
;
typename
ArrayTypes
<
DeviceType
>::
tdual_ffloat_2d
k_cut_ljsq
;
typename
ArrayTypes
<
DeviceType
>::
t_ffloat_2d
d_cut_ljsq
;
typename
ArrayTypes
<
DeviceType
>::
tdual_ffloat_2d
k_cut_coulsq
;
typename
ArrayTypes
<
DeviceType
>::
t_ffloat_2d
d_cut_coulsq
;
typename
ArrayTypes
<
DeviceType
>::
t_ffloat_1d_randomread
d_rtable
,
d_drtable
,
d_ftable
,
d_dftable
,
d_ctable
,
d_dctable
,
d_etable
,
d_detable
;
int
neighflag
;
int
nlocal
,
nall
,
eflag
,
vflag
;
double
special_coul
[
4
];
double
special_lj
[
4
];
double
qqrd2e
;
void
allocate
();
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
FULL
,
true
,
CoulLongTable
<
1
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALF
,
true
,
CoulLongTable
<
1
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
true
,
CoulLongTable
<
1
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
FULL
,
false
,
CoulLongTable
<
1
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALF
,
false
,
CoulLongTable
<
1
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
false
,
CoulLongTable
<
1
>
>
;
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
FULL
,
CoulLongTable
<
1
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
HALF
,
CoulLongTable
<
1
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
CoulLongTable
<
1
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute
<
PairLJCharmmCoulLongKokkos
,
CoulLongTable
<
1
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
FULL
,
true
,
CoulLongTable
<
0
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALF
,
true
,
CoulLongTable
<
0
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
true
,
CoulLongTable
<
0
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
FULL
,
false
,
CoulLongTable
<
0
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALF
,
false
,
CoulLongTable
<
0
>
>
;
friend
class
PairComputeFunctor
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
false
,
CoulLongTable
<
0
>
>
;
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
FULL
,
CoulLongTable
<
0
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
HALF
,
CoulLongTable
<
0
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute_neighlist
<
PairLJCharmmCoulLongKokkos
,
HALFTHREAD
,
CoulLongTable
<
0
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
EV_FLOAT
pair_compute
<
PairLJCharmmCoulLongKokkos
,
CoulLongTable
<
0
>
>
(
PairLJCharmmCoulLongKokkos
*
,
NeighListKokkos
<
DeviceType
>*
);
friend
void
pair_virial_fdotr_compute
<
PairLJCharmmCoulLongKokkos
>
(
PairLJCharmmCoulLongKokkos
*
);
};
}
#endif
#endif
/* ERROR/WARNING messages:
E: Cannot use Kokkos pair style with rRESPA inner/middle
UNDOCUMENTED
E: Cannot use chosen neighbor list style with lj/charmm/coul/long/kk
UNDOCUMENTED
*/
Event Timeline
Log In to Comment