From 23f37e6f63a685f847e39c3c4cc6777768a9e17d Mon Sep 17 00:00:00 2001 From: alphataubio Date: Thu, 18 Jul 2024 03:53:44 -0400 Subject: [PATCH] minor fixes --- src/KOKKOS/angle_spica_kokkos.cpp | 3 --- .../pair_lj_charmmfsw_coul_long_kokkos.cpp | 2 +- src/KOKKOS/pair_lj_spica_coul_long_kokkos.cpp | 27 ++----------------- src/KOKKOS/pair_lj_spica_coul_long_kokkos.h | 12 ++++----- 4 files changed, 9 insertions(+), 35 deletions(-) diff --git a/src/KOKKOS/angle_spica_kokkos.cpp b/src/KOKKOS/angle_spica_kokkos.cpp index 6c8c8b8aba..05da415f32 100644 --- a/src/KOKKOS/angle_spica_kokkos.cpp +++ b/src/KOKKOS/angle_spica_kokkos.cpp @@ -417,9 +417,6 @@ void AngleSPICAKokkos::coeff(int narg, char **arg) k_theta0.template modify(); k_repscale.template modify(); k_setflag.template modify(); - - Kokkos::printf("*** AngleSPICAKokkos::coeff\n"); - } /* ---------------------------------------------------------------------- diff --git a/src/KOKKOS/pair_lj_charmmfsw_coul_long_kokkos.cpp b/src/KOKKOS/pair_lj_charmmfsw_coul_long_kokkos.cpp index 23247bf20d..1c769e82a6 100644 --- a/src/KOKKOS/pair_lj_charmmfsw_coul_long_kokkos.cpp +++ b/src/KOKKOS/pair_lj_charmmfsw_coul_long_kokkos.cpp @@ -13,7 +13,7 @@ ------------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- - Contributing author: Mitch Murphy (alphataubio) + Contributing author: Mitch Murphy (alphataubio@gmail.com) Based on serial kspace lj-fsw sections (force-switched) provided by Robert Meissner and Lucio Colombi Ciacchi of Bremen University, Germany, diff --git a/src/KOKKOS/pair_lj_spica_coul_long_kokkos.cpp b/src/KOKKOS/pair_lj_spica_coul_long_kokkos.cpp index a205218a82..c0861e2c77 100644 --- a/src/KOKKOS/pair_lj_spica_coul_long_kokkos.cpp +++ b/src/KOKKOS/pair_lj_spica_coul_long_kokkos.cpp @@ -65,8 +65,6 @@ PairLJSPICACoulLongKokkos::~PairLJSPICACoulLongKokkos() memoryKK->destroy_kokkos(k_eatom,eatom); memoryKK->destroy_kokkos(k_vatom,vatom); memoryKK->destroy_kokkos(k_cutsq,cutsq); - //memoryKK->destroy_kokkos(k_cut_lj,cut_lj); - //memoryKK->destroy_kokkos(k_cut_ljsq,cut_ljsq); } } @@ -311,26 +309,10 @@ void PairLJSPICACoulLongKokkos::allocate() memoryKK->create_kokkos(k_cutsq,cutsq,n+1,n+1,"pair:cutsq"); d_cutsq = k_cutsq.template view(); - //d_cut_lj = typename AT::t_ffloat_2d("pair:cut_lj",n+1,n+1); - //d_cut_ljsq = typename AT::t_ffloat_2d("pair:cut_ljsq",n+1,n+1); - //d_cut_coulsq = typename AT::t_ffloat_2d("pair:cut_coulsq",n+1,n+1); - - k_params = Kokkos::DualView("PairLJSPICACoulLong::params",n+1,n+1); + k_params = Kokkos::DualView("PairLJSPICACoulLong::params",n+1,n+1); params = k_params.template view(); } -/* ---------------------------------------------------------------------- - global settings -------------------------------------------------------------------------- */ - -template -void PairLJSPICACoulLongKokkos::settings(int narg, char **arg) -{ - if (narg > 2) error->all(FLERR,"Illegal pair_style command"); - - PairLJSPICACoulLong::settings(1,arg); -} - /* ---------------------------------------------------------------------- init tables ------------------------------------------------------------------------- */ @@ -447,9 +429,6 @@ void PairLJSPICACoulLongKokkos::init_style() { PairLJSPICACoulLong::init_style(); - //Kokkos::deep_copy(d_cut_ljsq,cut_ljsq); - Kokkos::deep_copy(d_cut_coulsq,cut_coulsq); - // error if rRESPA with inner levels if (update->whichflag == 1 && utils::strmatch(update->integrate_style,"^respa")) { @@ -484,16 +463,14 @@ double PairLJSPICACoulLongKokkos::init_one(int i, int j) k_params.h_view(i,j).lj3 = lj3[i][j]; k_params.h_view(i,j).lj4 = lj4[i][j]; k_params.h_view(i,j).offset = offset[i][j]; - //k_params.h_view(i,j).cutsq = cutone*cutone; - //k_params.h_view(i,j).cut_lj = cut_lj[i][j]; k_params.h_view(i,j).cut_ljsq = cut_ljsq[i][j]; k_params.h_view(i,j).cut_coulsq = cut_coulsq; k_params.h_view(i,j).lj_type = lj_type[i][j]; k_params.h_view(j,i) = k_params.h_view(i,j); + if (i k_params; - typename Kokkos::DualView::t_dev_const_um params; - params_lj_coul m_params[MAX_TYPES_STACKPARAMS+1][MAX_TYPES_STACKPARAMS+1]; // hardwired to space for 12 atom types + Kokkos::DualView k_params; + typename Kokkos::DualView::t_dev_const_um params; + params_lj_spica_coul m_params[MAX_TYPES_STACKPARAMS+1][MAX_TYPES_STACKPARAMS+1]; // hardwired to space for 12 atom types 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];