| /dports/science/lammps/lammps-stable_29Sep2021/src/OPENMP/ |
| H A D | improper_fourier_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local
101 vb1x = x[i2][0] - x[i1][0]; in eval()
118 vb1x,vb1y,vb1z, in eval()
138 const double &vb1x, in add1_thr() argument
162 ay = vb1z*vb2x-vb1x*vb2z; in add1_thr()
163 az = vb1x*vb2y-vb1y*vb2x; in add1_thr()
194 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in add1_thr()
195 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in add1_thr()
221 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in add1_thr()
222 f2[2] = (dhax*vb1y - dhay*vb1x)*rar*a; in add1_thr()
[all …]
|
| H A D | dihedral_helix_omp.cpp | 114 vb1x = x[i1].x - x[i2].x; in eval()
136 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
143 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
147 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
154 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
180 cy = vb1z*vb2x - vb1x*vb2z; in eval()
181 cz = vb1x*vb2y - vb1y*vb2x; in eval()
215 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
219 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
227 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
[all …]
|
| H A D | dihedral_quadratic_omp.cpp | 114 vb1x = x[i1].x - x[i2].x; in eval()
136 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
143 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
147 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
154 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
180 cy = vb1z*vb2x - vb1x*vb2z; in eval()
181 cz = vb1x*vb2y - vb1y*vb2x; in eval()
222 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
226 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
234 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
[all …]
|
| H A D | dihedral_opls_omp.cpp | 113 vb1x = x[i1].x - x[i2].x; in eval()
135 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
142 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
146 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
153 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
179 cy = vb1z*vb2x - vb1x*vb2z; in eval()
180 cz = vb1x*vb2y - vb1y*vb2x; in eval()
219 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
223 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
231 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
[all …]
|
| H A D | dihedral_multi_harmonic_omp.cpp | 88 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in eval() local
112 vb1x = x[i1].x - x[i2].x; in eval()
134 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
141 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
145 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
152 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
204 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
208 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
216 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
252 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | dihedral_nharmonic_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in eval() local
111 vb1x = x[i1].x - x[i2].x; in eval()
133 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
140 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
144 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
151 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
208 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
212 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
220 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
256 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | improper_cvff_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in eval() local
110 vb1x = x[i1].x - x[i2].x; in eval()
132 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
139 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in eval()
143 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
150 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
232 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
236 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
244 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
280 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | improper_umbrella_omp.cpp | 88 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local
108 vb1x = x[i2].x - x[i1].x; in eval()
128 ay = vb1z*vb2x-vb1x*vb2z; in eval()
129 az = vb1x*vb2y-vb1y*vb2x; in eval()
160 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in eval()
161 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in eval()
195 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in eval()
196 f2[2] = (dhax*vb1y - dhay*vb1x)*rar*a; in eval()
240 vb1x = x[i1].x - x[i2].x; in eval()
253 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | improper_harmonic_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local
109 vb1x = x[i1].x - x[i2].x; in eval()
121 ss1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in eval()
131 c0 = (vb1x * vb3x + vb1y * vb3y + vb1z * vb3z) * r1 * r3; in eval()
132 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in eval()
174 sx2 = a22*vb2x + a23*vb3x + a12*vb1x; in eval()
178 f1[0] = a12*vb2x + a13*vb3x + a11*vb1x; in eval()
186 f4[0] = a23*vb2x + a33*vb3x + a13*vb1x; in eval()
222 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | improper_cossq_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local
107 vb1x = x[i2].x - x[i1].x; in eval()
110 rjisq = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z ; in eval()
125 cosphi = (vb3x*vb1x + vb3y*vb1y + vb3z*vb1z)/(rji * rlk); in eval()
154 clkji = vb3x*vb1x + vb3y*vb1y + vb3z*vb1z; in eval()
168 f1[0] = - vb3x + cfact3 * vb1x; in eval()
182 f3[0] = cfact2 * vb3x - vb1x; in eval()
248 -vb1x,-vb1y,-vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | dihedral_class2_omp.cpp | 116 vb1x = x[i1].x - x[i2].x; in eval()
138 r1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in eval()
200 double n123y = vb1z*vb2x - vb1x*vb2z; in eval()
201 double n123z = vb1x*vb2y - vb1y*vb2x; in eval()
215 sx1 = a11*vb1x + a12*vb2x + a13*vb3x; in eval()
216 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in eval()
217 sx12 = a13*vb1x + a23*vb2x + a33*vb3x; in eval()
273 dbonddr[0][0][0] = vb1x / r1; in eval()
276 dbonddr[0][1][0] = -vb1x / r1; in eval()
482 fabcd[0][0] += tk2 * vb1x; in eval()
[all …]
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/MOFFF/ |
| H A D | improper_inversion_harmonic.cpp | 82 vb1x = x[i2][0] - x[i1][0]; in compute()
85 rrvb1 = 1.0/sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in compute()
108 vb1x, vb1y, vb1z, rrvb1, rr2vb1); in compute()
110 vb1x, vb1y, vb1z, rrvb1, rr2vb1, in compute()
115 vb1x, vb1y, vb1z, rrvb1, rr2vb1, in compute()
147 rjl=vb1x*vb3x+vb1y*vb3y+vb1z*vb3z; in invang()
150 upx=vb2x*rrvb2+vb1x*rrvb1; in invang()
160 umx=vb2x*rrvb2-vb1x*rrvb1; in invang()
187 rudb = rjl-rup*(vb1x*upx+vb1y*upy+vb1z*upz); in invang()
189 rvdb = rjl-rum*(vb1x*umx+vb1y*umy+vb1z*umz); in invang()
[all …]
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/EXTRA-MOLECULE/ |
| H A D | improper_fourier.cpp | 76 vb1x = x[i2][0] - x[i1][0]; in compute()
93 vb1x, vb1y, vb1z, in compute()
99 vb1x, vb1y, vb1z, in compute()
104 vb1x, vb1y, vb1z); in compute()
129 ay = vb1z*vb2x-vb1x*vb2z; in addone()
130 az = vb1x*vb2y-vb1y*vb2x; in addone()
161 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in addone()
162 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in addone()
188 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in addone()
189 f2[2] = (dhax*vb1y - dhay*vb1x)*rar*a; in addone()
[all …]
|
| H A D | dihedral_quadratic.cpp | 89 vb1x = x[i1][0] - x[i2][0]; in compute()
111 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
118 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
122 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
129 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
155 cy = vb1z*vb2x - vb1x*vb2z; in compute()
156 cz = vb1x*vb2y - vb1y*vb2x; in compute()
197 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
201 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
209 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
[all …]
|
| H A D | dihedral_helix.cpp | 90 vb1x = x[i1][0] - x[i2][0]; in compute()
112 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
119 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
123 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
130 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
156 cy = vb1z*vb2x - vb1x*vb2z; in compute()
157 cz = vb1x*vb2y - vb1y*vb2x; in compute()
192 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
196 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
204 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
[all …]
|
| H A D | dihedral_nharmonic.cpp | 62 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in compute() local
89 vb1x = x[i1][0] - x[i2][0]; in compute()
111 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
118 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
122 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
129 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
186 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
190 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
198 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
234 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| H A D | improper_cossq.cpp | 58 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z ; in compute() local
82 vb1x = x[i2][0] - x[i1][0]; in compute()
85 rjisq = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z ; in compute()
100 cosphi = (vb3x*vb1x + vb3y*vb1y + vb3z*vb1z)/(rji * rlk); in compute()
129 clkji = vb3x*vb1x + vb3y*vb1y + vb3z*vb1z; in compute()
143 f1[0] = - vb3x + cfact3 * vb1x; in compute()
157 f3[0] = cfact2 * vb3x - vb1x; in compute()
223 -vb1x,-vb1y,-vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| H A D | dihedral_table_cut.cpp | 192 vb1x = x[i1][0] - x[i2][0]; in compute()
214 r1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
228 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
234 costh12 = (vb1x*vb2x + vb1y*vb2y + vb1z*vb2z) * r12c1; in compute()
271 double n123y = vb1z*vb2x - vb1x*vb2z; in compute()
272 double n123z = vb1x*vb2y - vb1y*vb2x; in compute()
287 sx1 = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
288 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
289 sx12 = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
338 dthetadr[0][0][0] = sc1 * ((t1 * vb1x) - (vb2x * r12c1)); in compute()
[all …]
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/MOLECULE/ |
| H A D | dihedral_opls.cpp | 88 vb1x = x[i1][0] - x[i2][0]; in compute()
110 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
117 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
121 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
128 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
154 cy = vb1z*vb2x - vb1x*vb2z; in compute()
155 cz = vb1x*vb2y - vb1y*vb2x; in compute()
194 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
198 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
206 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
[all …]
|
| H A D | improper_cvff.cpp | 55 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in compute() local
81 vb1x = x[i1][0] - x[i2][0]; in compute()
103 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
110 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
114 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
121 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
203 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
207 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
215 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
251 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| H A D | improper_umbrella.cpp | 62 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in compute() local
85 vb1x = x[i2][0] - x[i1][0]; in compute()
105 ay = vb1z*vb2x-vb1x*vb2z; in compute()
106 az = vb1x*vb2y-vb1y*vb2x; in compute()
137 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in compute()
138 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in compute()
172 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in compute()
173 f2[2] = (dhax*vb1y - dhay*vb1x)*rar*a; in compute()
217 vb1x = x[i1][0] - x[i2][0]; in compute()
230 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| H A D | improper_harmonic.cpp | 56 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in compute() local
81 vb1x = x[i1][0] - x[i2][0]; in compute()
93 ss1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
103 c0 = (vb1x * vb3x + vb1y * vb3y + vb1z * vb3z) * r1 * r3; in compute()
104 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in compute()
146 sx2 = a22*vb2x + a23*vb3x + a12*vb1x; in compute()
150 f1[0] = a12*vb2x + a13*vb3x + a11*vb1x; in compute()
158 f4[0] = a23*vb2x + a33*vb3x + a13*vb1x; in compute()
194 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| H A D | dihedral_multi_harmonic.cpp | 61 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in compute() local
88 vb1x = x[i1][0] - x[i2][0]; in compute()
110 sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute()
117 c0 = (vb1x*vb3x + vb1y*vb3y + vb1z*vb3z) * rb1*rb3; in compute()
121 b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in compute()
128 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
180 sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in compute()
184 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in compute()
192 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in compute()
228 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/KOKKOS/ |
| H A D | dihedral_opls_kokkos.cpp | 172 const F_FLOAT vb1x = x(i1,0) - x(i2,0); in operator ()() local
194 const F_FLOAT sb1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in operator ()()
205 const F_FLOAT b1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; in operator ()()
212 F_FLOAT ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in operator ()()
238 const F_FLOAT cy = vb1z*vb2x - vb1x*vb2z; in operator ()()
239 const F_FLOAT cz = vb1x*vb2y - vb1y*vb2x; in operator ()()
279 const F_FLOAT sx2 = a12*vb1x + a22*vb2x + a23*vb3x; in operator ()()
284 f1[0] = a11*vb1x + a12*vb2x + a13*vb3x; in operator ()()
292 f4[0] = a13*vb1x + a23*vb2x + a33*vb3x; in operator ()()
448 v[0] = vb1x*f1[0] + vb2x*f3[0] + (vb3x+vb2x)*f4[0]; in ev_tally()
[all …]
|
| H A D | improper_harmonic_kokkos.cpp | 171 const F_FLOAT vb1x = x(i1,0) - x(i2,0); in operator ()() local
183 const F_FLOAT ss1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in operator ()()
193 const F_FLOAT c0 = (vb1x * vb3x + vb1y * vb3y + vb1z * vb3z) * r1 * r3; in operator ()()
194 const F_FLOAT c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in operator ()()
237 const F_FLOAT sx2 = a22*vb2x + a23*vb3x + a12*vb1x; in operator ()()
242 f1[0] = a12*vb2x + a13*vb3x + a11*vb1x; in operator ()()
250 f4[0] = a23*vb2x + a33*vb3x + a13*vb1x; in operator ()()
286 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in operator ()()
391 v[0] = vb1x*f1[0] + vb2x*f3[0] + (vb3x+vb2x)*f4[0]; in ev_tally()
394 v[3] = vb1x*f1[1] + vb2x*f3[1] + (vb3x+vb2x)*f4[1]; in ev_tally()
[all …]
|