| /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
109 vb2z = x[i3][2] - x[i1][2]; in eval()
119 vb2x,vb2y,vb2z, in eval()
143 const double &vb2z, in add1_thr() argument
161 ax = vb1y*vb2z-vb1z*vb2y; in add1_thr()
162 ay = vb1z*vb2x-vb1x*vb2z; in add1_thr()
196 projhfg += (vb3x*vb2x+vb3y*vb2y+vb3z*vb2z) / in add1_thr()
197 sqrt(vb2x*vb2x+vb2y*vb2y+vb2z*vb2z); in add1_thr()
224 f3[0] = (-dhay*vb2z + dhaz*vb2y)*rar*a; in add1_thr()
225 f3[1] = (-dhaz*vb2x + dhax*vb2z)*rar*a; in add1_thr()
[all …]
|
| H A D | dihedral_helix_omp.cpp | 122 vb2z = x[i3].z - x[i2].z; in eval()
126 vb2zm = -vb2z; in eval()
137 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
149 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
154 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
179 cx = vb1y*vb2z - vb1z*vb2y; in eval()
180 cy = vb1z*vb2x - vb1x*vb2z; in eval()
217 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
221 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
229 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in eval()
[all …]
|
| H A D | dihedral_quadratic_omp.cpp | 122 vb2z = x[i3].z - x[i2].z; in eval()
126 vb2zm = -vb2z; in eval()
137 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
149 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
154 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
179 cx = vb1y*vb2z - vb1z*vb2y; in eval()
180 cy = vb1z*vb2x - vb1x*vb2z; in eval()
224 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
228 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
236 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in eval()
[all …]
|
| H A D | dihedral_opls_omp.cpp | 121 vb2z = x[i3].z - x[i2].z; in eval()
125 vb2zm = -vb2z; in eval()
136 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
148 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
153 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
178 cx = vb1y*vb2z - vb1z*vb2y; in eval()
179 cy = vb1z*vb2x - vb1x*vb2z; in eval()
221 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
225 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
233 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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
120 vb2z = x[i3].z - x[i2].z; in eval()
124 vb2zm = -vb2z; in eval()
135 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
147 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
152 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
206 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
210 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
218 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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
119 vb2z = x[i3].z - x[i2].z; in eval()
123 vb2zm = -vb2z; in eval()
134 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
146 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
151 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
210 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
214 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
222 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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
118 vb2z = x[i3].z - x[i2].z; in eval()
122 vb2zm = -vb2z; in eval()
133 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
145 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
150 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in eval()
234 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
238 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
246 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in eval()
280 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| H A D | dihedral_class2_omp.cpp | 124 vb2z = x[i3].z - x[i2].z; in eval()
128 vb2zm = -vb2z; in eval()
140 r2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in eval()
158 costh12 = (vb1x*vb2x + vb1y*vb2y + vb1z*vb2z) * r12c1; in eval()
199 double n123x = vb1y*vb2z - vb1z*vb2y; in eval()
200 double n123y = vb1z*vb2x - vb1x*vb2z; in eval()
221 sz1 = a11*vb1z + a12*vb2z + a13*vb3z; in eval()
222 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in eval()
223 sz12 = a13*vb1z + a23*vb2z + a33*vb3z; in eval()
284 dbonddr[1][1][2] = vb2z / r2; in eval()
[all …]
|
| H A D | improper_umbrella_omp.cpp | 88 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local
116 vb2z = x[i3].z - x[i1].z; in eval()
127 ax = vb1y*vb2z-vb1z*vb2y; in eval()
128 ay = vb1z*vb2x-vb1x*vb2z; in eval()
162 projhfg += (vb3x*vb2x+vb3y*vb2y+vb3z*vb2z) / in eval()
163 sqrt(vb2x*vb2x+vb2y*vb2y+vb2z*vb2z); in eval()
198 f3[0] = (-dhay*vb2z + dhaz*vb2y)*rar*a; in eval()
199 f3[1] = (-dhaz*vb2x + dhax*vb2z)*rar*a; in eval()
246 vb2z = x[i3].z - x[i2].z; 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
115 vb2z = x[i3].z - x[i2].z; in eval()
122 ss2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in eval()
132 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in eval()
133 c2 = -(vb3x * vb2x + vb3y * vb2y + vb3z * vb2z) * r3 * r2; in eval()
176 sz2 = a22*vb2z + a23*vb3z + a12*vb1z; in eval()
180 f1[2] = a12*vb2z + a13*vb3z + a11*vb1z; in eval()
188 f4[2] = a23*vb2z + a33*vb3z + a13*vb1z; in eval()
222 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/MOFFF/ |
| H A D | improper_inversion_harmonic.cpp | 92 vb2z = x[i3][2] - x[i1][2]; in compute()
93 rrvb2 = 1.0/sqrt(vb2x*vb2x+vb2y*vb2y+vb2z*vb2z); in compute()
107 vb2x, vb2y, vb2z, rrvb2, rr2vb2, in compute()
112 vb2x, vb2y, vb2z, rrvb2, rr2vb2); in compute()
114 vb2x, vb2y, vb2z, rrvb2, rr2vb2, in compute()
146 rjk=vb3x*vb2x+vb3y*vb2y+vb3z*vb2z; in invang()
152 upz=vb2z*rrvb2+vb1z*rrvb1; in invang()
162 umz=vb2z*rrvb2-vb1z*rrvb1; in invang()
186 rucb = rjk-rup*(vb2x*upx+vb2y*upy+vb2z*upz); in invang()
246 rb3z = vb1z - vb2z; in invang()
[all …]
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/EXTRA-MOLECULE/ |
| H A D | improper_fourier.cpp | 84 vb2z = x[i3][2] - x[i1][2]; in compute()
94 vb2x, vb2y, vb2z, in compute()
100 vb2x, vb2y, vb2z); in compute()
102 vb2x, vb2y, vb2z, in compute()
128 ax = vb1y*vb2z-vb1z*vb2y; in addone()
129 ay = vb1z*vb2x-vb1x*vb2z; in addone()
163 projhfg += (vb3x*vb2x+vb3y*vb2y+vb3z*vb2z) / in addone()
164 sqrt(vb2x*vb2x+vb2y*vb2y+vb2z*vb2z); in addone()
191 f3[0] = (-dhay*vb2z + dhaz*vb2y)*rar*a; in addone()
192 f3[1] = (-dhaz*vb2x + dhax*vb2z)*rar*a; in addone()
[all …]
|
| H A D | dihedral_quadratic.cpp | 97 vb2z = x[i3][2] - x[i2][2]; in compute()
101 vb2zm = -vb2z; in compute()
112 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
124 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
129 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
154 cx = vb1y*vb2z - vb1z*vb2y; in compute()
155 cy = vb1z*vb2x - vb1x*vb2z; in compute()
199 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
203 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
211 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in compute()
[all …]
|
| H A D | dihedral_helix.cpp | 98 vb2z = x[i3][2] - x[i2][2]; in compute()
102 vb2zm = -vb2z; in compute()
113 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
125 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
130 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
155 cx = vb1y*vb2z - vb1z*vb2y; in compute()
156 cy = vb1z*vb2x - vb1x*vb2z; in compute()
194 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
198 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
206 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in compute()
[all …]
|
| H A D | dihedral_table_cut.cpp | 200 vb2z = x[i3][2] - x[i2][2]; in compute()
204 vb2zm = -vb2z; in compute()
216 r2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
234 costh12 = (vb1x*vb2x + vb1y*vb2y + vb1z*vb2z) * r12c1; in compute()
270 double n123x = vb1y*vb2z - vb1z*vb2y; in compute()
271 double n123y = vb1z*vb2x - vb1x*vb2z; in compute()
293 sz1 = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
294 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
295 sz12 = a13*vb1z + a23*vb2z + a33*vb3z; in compute()
340 dthetadr[0][0][2] = sc1 * ((t1 * vb1z) - (vb2z * r12c1)); 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
97 vb2z = x[i3][2] - x[i2][2]; in compute()
101 vb2zm = -vb2z; in compute()
112 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
124 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
129 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
188 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
192 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
200 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in compute()
234 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/MOLECULE/ |
| H A D | dihedral_opls.cpp | 96 vb2z = x[i3][2] - x[i2][2]; in compute()
100 vb2zm = -vb2z; in compute()
111 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
123 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
128 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
153 cx = vb1y*vb2z - vb1z*vb2y; in compute()
154 cy = vb1z*vb2x - vb1x*vb2z; in compute()
196 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
200 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
208 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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
89 vb2z = x[i3][2] - x[i2][2]; in compute()
93 vb2zm = -vb2z; in compute()
104 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
116 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
121 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
205 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
209 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
217 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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
93 vb2z = x[i3][2] - x[i1][2]; in compute()
104 ax = vb1y*vb2z-vb1z*vb2y; in compute()
105 ay = vb1z*vb2x-vb1x*vb2z; in compute()
139 projhfg += (vb3x*vb2x+vb3y*vb2y+vb3z*vb2z) / in compute()
140 sqrt(vb2x*vb2x+vb2y*vb2y+vb2z*vb2z); in compute()
175 f3[0] = (-dhay*vb2z + dhaz*vb2y)*rar*a; in compute()
176 f3[1] = (-dhaz*vb2x + dhax*vb2z)*rar*a; in compute()
223 vb2z = x[i3][2] - x[i2][2]; 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
87 vb2z = x[i3][2] - x[i2][2]; in compute()
94 ss2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
104 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in compute()
105 c2 = -(vb3x * vb2x + vb3y * vb2y + vb3z * vb2z) * r3 * r2; in compute()
148 sz2 = a22*vb2z + a23*vb3z + a12*vb1z; in compute()
152 f1[2] = a12*vb2z + a13*vb3z + a11*vb1z; in compute()
160 f4[2] = a23*vb2z + a33*vb3z + a13*vb1z; 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
96 vb2z = x[i3][2] - x[i2][2]; in compute()
100 vb2zm = -vb2z; in compute()
111 sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute()
123 b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in compute()
128 ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in compute()
182 sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in compute()
186 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in compute()
194 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; 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 | 180 const F_FLOAT vb2z = x(i3,2) - x(i2,2); in operator ()() local
184 const F_FLOAT vb2zm = -vb2z; in operator ()()
195 const F_FLOAT sb2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in operator ()()
207 const F_FLOAT b2mag2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z; in operator ()()
212 F_FLOAT ctmp = vb1x*vb2x + vb1y*vb2y + vb1z*vb2z; in operator ()()
237 const F_FLOAT cx = vb1y*vb2z - vb1z*vb2y; in operator ()()
238 const F_FLOAT cy = vb1z*vb2x - vb1x*vb2z; in operator ()()
281 const F_FLOAT sz2 = a12*vb1z + a22*vb2z + a23*vb3z; in operator ()()
286 f1[2] = a11*vb1z + a12*vb2z + a13*vb3z; in operator ()()
294 f4[2] = a13*vb1z + a23*vb2z + a33*vb3z; in operator ()()
[all …]
|
| H A D | improper_harmonic_kokkos.cpp | 177 const F_FLOAT vb2z = x(i3,2) - x(i2,2); in operator ()() local
184 const F_FLOAT ss2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in operator ()()
194 const F_FLOAT c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in operator ()()
195 const F_FLOAT c2 = -(vb3x * vb2x + vb3y * vb2y + vb3z * vb2z) * r3 * r2; in operator ()()
239 const F_FLOAT sz2 = a22*vb2z + a23*vb3z + a12*vb1z; in operator ()()
244 f1[2] = a12*vb2z + a13*vb3z + a11*vb1z; in operator ()()
252 f4[2] = a23*vb2z + a33*vb3z + a13*vb1z; in operator ()()
286 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in operator ()()
363 const F_FLOAT &vb2x, const F_FLOAT &vb2y, const F_FLOAT &vb2z, in ev_tally() argument
393 v[2] = vb1z*f1[2] + vb2z*f3[2] + (vb3z+vb2z)*f4[2]; in ev_tally()
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/ |
| H A D | compute_improper_local.cpp | 113 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in compute_impropers() local
185 vb2z = x[atom3][2] - x[atom2][2]; in compute_impropers()
186 domain->minimum_image(vb2x,vb2y,vb2z); in compute_impropers()
194 ss2 = 1.0 / (vb2x*vb2x + vb2y*vb2y + vb2z*vb2z); in compute_impropers()
202 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in compute_impropers()
203 c2 = -(vb3x * vb2x + vb3y * vb2y + vb3z * vb2z) * r3 * r2; in compute_impropers()
|
| /dports/science/lammps/lammps-stable_29Sep2021/src/INTEL/ |
| H A D | improper_harmonic_intel.cpp | 193 const flt_t vb2z = x[i3].z - x[i2].z; local
200 flt_t ss2 = vb2x*vb2x + vb2y*vb2y + vb2z*vb2z;
214 const flt_t c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2;
215 const flt_t c2 = -(vb3x * vb2x + vb3y * vb2y + vb3z * vb2z) * r3 * r2;
262 const flt_t sz2 = a22*vb2z + a23*vb3z + a12*vb1z;
266 const flt_t f1z = a12*vb2z + a13*vb3z + a11*vb1z;
274 const flt_t f4z = a23*vb2z + a33*vb3z + a13*vb1z;
319 f4y, f4z, vb1x, vb1y, vb1z, vb2x, vb2y, vb2z,
325 f4y, f4z, vb1x, vb1y, vb1z, vb2x, vb2y, vb2z,
|