PairwiseInteractionEnergies.h
上传用户:chinafayin
上传日期:2022-04-05
资源大小:153k
文件大小:2k
源码类别:

并行计算

开发平台:

Visual C++

  1. /*
  2.     FastGrid (formerly AutoGrid)
  3.     Copyright (C) 2009 The Scripps Research Institute. All rights reserved.
  4.     Copyright (C) 2009 Masaryk University. All rights reserved.
  5.     AutoGrid is a Trade Mark of The Scripps Research Institute.
  6.     This program is free software; you can redistribute it and/or
  7.     modify it under the terms of the GNU General Public License
  8.     as published by the Free Software Foundation; either version 2
  9.     of the License, or (at your option) any later version.
  10.     This program is distributed in the hope that it will be useful,
  11.     but WITHOUT ANY WARRANTY; without even the implied warranty of
  12.     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13.     GNU General Public License for more details.
  14.     You should have received a copy of the GNU General Public License
  15.     along with this program; if not, write to the Free Software
  16.     Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  17. */
  18. #pragma once
  19. #include "GridMap.h"
  20. // Calculates the lookup table of the pairwise interaction energies
  21. class PairwiseInteractionEnergies
  22. {
  23. public:
  24.     struct LookupTable
  25.     {
  26.         // MAX_DIST is not NBCUTOFF times 100 as it should be, it's a power of two for a little faster memory access,
  27.         // values beyond this threshold are unused
  28.         double table[NUM_RECEPTOR_TYPES][MAX_DIST][MAX_MAPS];
  29.     };
  30.     PairwiseInteractionEnergies();
  31.     ~PairwiseInteractionEnergies();
  32.     // calculates the table of pair-wise interaction energies
  33.     void calculate(const GridMapList &gridmaps, LogFile &logFile,
  34.                    int numReceptorTypes, const char (&receptorTypes)[NUM_RECEPTOR_TYPES][3], double rSmooth);
  35.     // returns the precalculated value
  36.     double operator ()(int atomType, int indexR, int mapIndex) const
  37.     {
  38.         return energyLookup->table[atomType][indexR][mapIndex];
  39.     }
  40. private:
  41.     LookupTable *energyLookup;
  42. };