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
/* -*- 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 FIX_CLASS
FixStyle(ave/correlate/long,FixAveCorrelateLong)
#else
#ifndef LMP_FIX_AVE_CORRELATE_LONG_H
#define LMP_FIX_AVE_CORRELATE_LONG_H
#include "fix.h"
namespace LAMMPS_NS {
class FixAveCorrelateLong : public Fix {
public:
FixAveCorrelateLong(class LAMMPS *, int, char **);
~FixAveCorrelateLong();
int setmask();
void init();
void setup(int);
void end_of_step();
void write_restart(FILE *);
void restart(char *);
double memory_usage();
double *t; // Time steps for result arrays
double **f; // Result arrays
unsigned int npcorr;
private:
// NOT OPTIMAL: shift2 and accumulator2 only needed in cross-correlations
double ***shift, *** shift2;
double ***correlation;
double **accumulator, **accumulator2;
unsigned long int **ncorrelation;
unsigned int *naccumulator;
unsigned int *insertindex;
int numcorrelators; // Recommended 20
unsigned int p; // Points per correlator (recommended 16)
unsigned int m; // Num points for average (recommended 2; p mod m = 0)
unsigned int dmin; // Min distance between ponts for correlators k>0; dmin=p/m
int length; // Length of result arrays
int kmax; // Maximum correlator attained during simulation
int me,nvalues;
int nfreq;
bigint nvalid,nvalid_last,last_accumulated_step;
int *which,*argindex,*value2index;
char **ids;
FILE *fp;
int type,startstep,overwrite;
long filepos;
int npair; // number of correlation pairs to calculate
double *values;
void accumulate();
void evaluate();
bigint nextvalid();
void add(const int i, const double w, const int k = 0);
void add(const int i, const double wA, const double wB, const int k = 0);
};
}
#endif
#endif
/* ERROR/WARNING messages:
E: Illegal ... command
Self-explanatory. Check the input script syntax and compare to the
documentation for the command. You can use -echo screen as a
command-line option when running LAMMPS to see the offending line.
E: Cannot open fix ave/correlate/long file %s
The specified file cannot be opened. Check that the path and name are
correct.
E: Compute ID for fix ave/correlate/long does not exist
Self-explanatory.
E: Fix ave/correlate/long compute does not calculate a scalar
Self-explanatory.
E: Fix ave/correlate/long compute does not calculate a vector
Self-explanatory.
E: Fix ave/correlate/long compute vector is accessed out-of-range
The index for the vector is out of bounds.
E: Fix ID for fix ave/correlate/long does not exist
Self-explanatory.
E: Fix ave/correlate/long fix does not calculate a scalar
Self-explanatory.
E: Fix ave/correlate/long fix does not calculate a vector
Self-explanatory.
E: Fix ave/correlate/long fix vector is accessed out-of-range
The index for the vector is out of bounds.
E: Fix for fix ave/correlate/long not computed at compatible time
Fixes generate their values on specific timesteps. Fix ave/correlate/long
is requesting a value on a non-allowed timestep.
E: Variable name for fix ave/correlate/long does not exist
Self-explanatory.
E: Fix ave/correlate/long variable is not equal-style variable
Self-explanatory.
E: Invalid timestep reset for fix ave/correlate/long
Resetting the timestep has invalidated the sequence of timesteps this
fix needs to process.
*/