BrookIntegrateLangevinStepKernel.h 4.5 KB
Newer Older
Mark Friedrichs's avatar
Mark Friedrichs committed
1
2
#ifndef OPENMM_BROOK_INTEGRATE_LANGEVIN_STEP_KERNEL_H_
#define OPENMM_BROOK_INTEGRATE_LANGEVIN_STEP_KERNEL_H_
Mark Friedrichs's avatar
Mark Friedrichs committed
3
4
5
6
7
8
9
10
11

/* -------------------------------------------------------------------------- *
 *                                   OpenMM                                   *
 * -------------------------------------------------------------------------- *
 * This is part of the OpenMM molecular simulation toolkit originating from   *
 * Simbios, the NIH National Center for Physics-Based Simulation of           *
 * Biological Structures at Stanford, funded under the NIH Roadmap for        *
 * Medical Research, grant U54 GM072970. See https://simtk.org.               *
 *                                                                            *
12
13
 * Portions copyright (c) 2009 Stanford University and the Authors.           *
 * Authors: Mark Friedrichs, Mike Houston                                     *
Mark Friedrichs's avatar
Mark Friedrichs committed
14
15
 * Contributors:                                                              *
 *                                                                            *
16
17
18
19
 * This program 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.                                        *
Mark Friedrichs's avatar
Mark Friedrichs committed
20
 *                                                                            *
21
22
23
24
 * This program 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.                        *
Mark Friedrichs's avatar
Mark Friedrichs committed
25
 *                                                                            *
26
27
 * You should have received a copy of the GNU Lesser General Public License   *
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.      *
Mark Friedrichs's avatar
Mark Friedrichs committed
28
29
 * -------------------------------------------------------------------------- */

30
#include "openmm/kernels.h"
Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
31
#include "OpenMMBrookInterface.h"
32
#include "BrookLangevinDynamics.h"
Mark Friedrichs's avatar
Mark Friedrichs committed
33
#include "BrookShakeAlgorithm.h"
Mark Friedrichs's avatar
Mark Friedrichs committed
34
#include "BrookRandomNumberGenerator.h"
Mark Friedrichs's avatar
Mark Friedrichs committed
35
36
37
38

namespace OpenMM {

/**
Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
39
 * Performs Langevin integration step
Mark Friedrichs's avatar
Mark Friedrichs committed
40
41
 */

Mark Friedrichs's avatar
Mark Friedrichs committed
42
class BrookIntegrateLangevinStepKernel : public IntegrateLangevinStepKernel {
Mark Friedrichs's avatar
Mark Friedrichs committed
43

44
   public:
Mark Friedrichs's avatar
Mark Friedrichs committed
45

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
46
47
48
49
50
      // return values

      static const int DefaultReturnValue = 0;
      static const int ErrorReturnValue   = -1; 

Mark Friedrichs's avatar
Mark Friedrichs committed
51
      /**
Mark Friedrichs's avatar
Mark Friedrichs committed
52
       * BrookIntegrateLangevinStepKernel constructor
Mark Friedrichs's avatar
Mark Friedrichs committed
53
54
55
56
57
58
       * 
       * @param name        name of the stream to create
       * @param platform    platform
       *
       */
  
Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
59
      BrookIntegrateLangevinStepKernel( std::string name, const Platform& platform, OpenMMBrookInterface& openMMBrookInterface, System& system );
60

Mark Friedrichs's avatar
Mark Friedrichs committed
61
      /**
Mark Friedrichs's avatar
Mark Friedrichs committed
62
       * BrookIntegrateLangevinStepKernel destructor
Mark Friedrichs's avatar
Mark Friedrichs committed
63
64
65
       * 
       */
  
Mark Friedrichs's avatar
Mark Friedrichs committed
66
      ~BrookIntegrateLangevinStepKernel();
67
68

      /** 
Mark Friedrichs's avatar
Mark Friedrichs committed
69
       * Initialize the kernel, setting up all parameters related to integrator.
70
       * 
Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
71
72
       * @param system             System reference
       * @param integrator         LangevinIntegrator reference
Mark Friedrichs's avatar
Mark Friedrichs committed
73
       */
Mark Friedrichs's avatar
Mark Friedrichs committed
74

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
75
76
      void initialize( const System& system, const LangevinIntegrator& integrator );

Mark Friedrichs's avatar
Mark Friedrichs committed
77
      /** 
Mark Friedrichs's avatar
Mark Friedrichs committed
78
       * Execute kernel
Mark Friedrichs's avatar
Mark Friedrichs committed
79
       * 
80
       * @param context            ContextImpl reference
Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
81
       * @param integrator         LangevinIntegrator reference
Mark Friedrichs's avatar
Mark Friedrichs committed
82
       *
83
84
       */

85
      void execute( ContextImpl& context, const LangevinIntegrator& integrator );
86

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
      /** 
       * Set log file reference
       * 
       * @param  log file reference
       *
       * @return DefaultReturnValue
       *
       */
          
      int setLog( FILE* log );

      /*  
       * Get contents of object
       *
       * @param level of dump
       *
       * @return string containing contents
       *
       * */
          
      std::string getContents( int level ) const;

      /** 
       * Get log file reference
       * 
       * @return  log file reference
       *
       */
          
      FILE* getLog( void ) const;

Mark Friedrichs's avatar
Mark Friedrichs committed
118
   protected:
119

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
120
121
      FILE* _log;

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
122
      BrookLangevinDynamics*        _brookLangevinDynamics;
Mark Friedrichs's avatar
Mark Friedrichs committed
123
124
      BrookShakeAlgorithm*          _brookShakeAlgorithm;
      BrookRandomNumberGenerator*   _brookRandomNumberGenerator;
Mark Friedrichs's avatar
Mark Friedrichs committed
125

Mark Friedrichs's avatar
Mods  
Mark Friedrichs committed
126
127
128
129
130
131
132
133
      // interface

      OpenMMBrookInterface& _openMMBrookInterface;

      // System reference

      System& _system;

Mark Friedrichs's avatar
Mark Friedrichs committed
134
135
136
137
};

} // namespace OpenMM

Mark Friedrichs's avatar
Mark Friedrichs committed
138
#endif /* OPENMM_BROOK_INTEGRATE_LANGEVIN_STEP_KERNEL_H_ */