1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
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
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
|
/*
* $Id: update.h,v 1.45 2005/08/04 13:31:03 spoel Exp $
*
* This source code is part of
*
* G R O M A C S
*
* GROningen MAchine for Chemical Simulations
*
* VERSION 3.2.0
* Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.
* Copyright (c) 1991-2000, University of Groningen, The Netherlands.
* Copyright (c) 2001-2004, The GROMACS development team,
* check out http://www.gromacs.org for more information.
* This program is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License
* as published by the Free Software Foundation; either version 2
* of the License, or (at your option) any later version.
*
* If you want to redistribute modifications, please consider that
* scientific software is very special. Version control is crucial -
* bugs must be traceable. We will be happy to consider code for
* inclusion in the official distribution, but derived work must not
* be called official GROMACS. Details are found in the README & COPYING
* files - if they are missing, get the official version at www.gromacs.org.
*
* To help us fund GROMACS development, we humbly ask that you cite
* the papers on the package - you can find them in the top README file.
*
* For more info, check our website at http://www.gromacs.org
*
* And Hey:
* Gromacs Runs On Most of All Computer Systems
*/
#ifndef _update_h
#define _update_h
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include "typedefs.h"
#include "mshift.h"
#include "tgroup.h"
#include "network.h"
#include "force.h"
#include "pull.h"
extern void update(int natoms, /* number of atoms in simulation */
int start,
int homenr, /* number of home particles */
int step,
real *dvdlambda, /* FEP stuff */
t_inputrec *inputrec, /* input record and box stuff */
t_mdatoms *md,
t_state *state,
t_graph *graph,
rvec force[],/* forces on home particles */
rvec vold[], /* Old velocities */
t_topology *top,
t_groups *grps,
tensor vir_part,
t_commrec *cr,
t_nrnb *nrnb,
t_edsamyn *edyn,
t_pull *pulldata,
bool bNEMD,
bool bDoUpdate,
bool bFirstStep,
rvec *shakefirst_x,
tensor pres);
/* Return TRUE if OK, FALSE in case of Shake Error */
extern void calc_ke_part(int start,int homenr,rvec v[],
t_grpopts *opts,t_mdatoms *md,
t_groups *grps,t_nrnb *nrnb,
real lambda);
/*
* Compute the partial kinetic energy for home particles;
* will be accumulated in the calling routine.
* The tensor is
*
* Ekin = SUM(i) 0.5 m[i] v[i] (x) v[i]
*
* use v[i] = v[i] - u[i] when calculating temperature
*
* u must be accumulated already.
*
* Now also computes the contribution of the kinetic energy to the
* free energy
*
*/
extern void calc_ke_part_visc(int start,int homenr,
matrix box,rvec x[],rvec v[],
t_grpopts *opts,t_mdatoms *md,
t_groups *grps,t_nrnb *nrnb,
real lambda);
/* The same as calc_ke_part, but for viscosity calculations.
* The cosine velocity profile is excluded from the kinetic energy.
* The new amplitude of the velocity profile is calculated for this
* node and stored in grps->cosacc.mvcos.
*/
extern void init_sd_consts(int ngtc,real tau_t[],real dt);
/* Initialization of the SD constants (obviously). */
/* Routines from coupling.c to do with Temperature, Pressure and coupling
* algorithms.
*/
extern real run_aver(real old,real cur,int step,int nmem);
extern void berendsen_tcoupl(t_grpopts *opts,t_groups *grps,real dt);
extern void nosehoover_tcoupl(t_grpopts *opts,t_groups *grps,real dt,
real xi[]);
/* Compute temperature scaling. For Nose-Hoover it is done in update. */
/* Set reference temp for simulated annealing at time t*/
extern void update_annealing_target_temp(t_grpopts *opts,real t);
extern real calc_temp(real ekin,real nrdf);
/* Calculate the temperature */
extern void calc_pres(int ePBC,matrix box,
tensor ekin,tensor vir,tensor pres,real Elr);
/* Calculate the pressure. Unit of pressure is bar, If Elr != 0
* a long range correction based on Ewald/PPPM is made (see c-code)
*/
extern void parrinellorahman_pcoupl(t_inputrec *ir,int step,tensor pres,
tensor box,tensor boxv,tensor M,
bool bFirstStep);
extern void berendsen_pcoupl(t_inputrec *ir,int step,tensor pres,matrix box,
matrix mu);
extern void berendsen_pscale(matrix mu,
matrix box,int start,int nr_atoms,
rvec x[],unsigned short cFREEZE[],
t_nrnb *nrnb,ivec nFreeze[]);
extern void correct_ekin(FILE *log,int start,int end,rvec v[],
rvec vcm,real mass[],real tmass,tensor ekin);
/* Correct ekin for vcm */
#endif /* _update_h */
|
