/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 103 vb1z = x[i2][2] - x[i1][2]; in eval() 118 vb1x,vb1y,vb1z, in eval() 140 const double &vb1z, in add1_thr() argument 161 ax = vb1y*vb2z-vb1z*vb2y; in add1_thr() 162 ay = vb1z*vb2x-vb1x*vb2z; 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() 220 f2[0] = (dhay*vb1z - dhaz*vb1y)*rar*a; in add1_thr() 221 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in add1_thr() [all …]
|
H A D | dihedral_helix_omp.cpp | 116 vb1z = x[i1].z - x[i2].z; 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() 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 | 116 vb1z = x[i1].z - x[i2].z; 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() 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 | 115 vb1z = x[i1].z - x[i2].z; 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() 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 114 vb1z = x[i1].z - x[i2].z; 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() 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 113 vb1z = x[i1].z - x[i2].z; 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() 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 112 vb1z = x[i1].z - x[i2].z; 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() 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 | improper_umbrella_omp.cpp | 88 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local 110 vb1z = x[i2].z - x[i1].z; in eval() 127 ax = vb1y*vb2z-vb1z*vb2y; in eval() 128 ay = vb1z*vb2x-vb1x*vb2z; in eval() 160 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in eval() 161 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in eval() 194 f2[0] = (dhay*vb1z - dhaz*vb1y)*rar*a; in eval() 195 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in eval() 242 vb1z = x[i1].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 111 vb1z = x[i1].z - x[i2].z; 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() 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()
|
H A D | improper_cossq_omp.cpp | 87 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in eval() local 109 vb1z = x[i2].z - x[i1].z; 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() 170 f1[2] = - vb3z + cfact3 * vb1z; in eval() 184 f3[2] = cfact2 * vb3z - vb1z; in eval() 248 -vb1x,-vb1y,-vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,thr); in eval()
|
H A D | dihedral_class2_omp.cpp | 118 vb1z = x[i1].z - x[i2].z; in eval() 138 r1mag2 = vb1x*vb1x + vb1y*vb1y + vb1z*vb1z; 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() 275 dbonddr[0][0][2] = vb1z / r1; in eval() 278 dbonddr[0][1][2] = -vb1z / r1; in eval() 484 fabcd[0][2] += tk2 * vb1z; in eval() [all …]
|
/dports/science/lammps/lammps-stable_29Sep2021/src/MOFFF/ |
H A D | improper_inversion_harmonic.cpp | 84 vb1z = x[i2][2] - x[i1][2]; 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() 152 upz=vb2z*rrvb2+vb1z*rrvb1; in invang() 162 umz=vb2z*rrvb2-vb1z*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 | 78 vb1z = x[i2][2] - x[i1][2]; in compute() 93 vb1x, vb1y, vb1z, in compute() 99 vb1x, vb1y, vb1z, in compute() 104 vb1x, vb1y, vb1z); in compute() 128 ax = vb1y*vb2z-vb1z*vb2y; in addone() 129 ay = vb1z*vb2x-vb1x*vb2z; in addone() 161 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in addone() 162 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in addone() 187 f2[0] = (dhay*vb1z - dhaz*vb1y)*rar*a; in addone() 188 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in addone() [all …]
|
H A D | dihedral_quadratic.cpp | 91 vb1z = x[i1][2] - x[i2][2]; 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() 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 | 92 vb1z = x[i1][2] - x[i2][2]; 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() 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_nharmonic.cpp | 62 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z,vb2xm,vb2ym,vb2zm; in compute() local 91 vb1z = x[i1][2] - x[i2][2]; 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() 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()
|
H A D | improper_cossq.cpp | 58 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z ; in compute() local 84 vb1z = x[i2][2] - x[i1][2]; 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() 145 f1[2] = - vb3z + cfact3 * vb1z; in compute() 159 f3[2] = cfact2 * vb3z - vb1z; in compute() 223 -vb1x,-vb1y,-vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in compute()
|
H A D | dihedral_table_cut.cpp | 194 vb1z = x[i1][2] - x[i2][2]; 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() 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 …]
|
/dports/science/lammps/lammps-stable_29Sep2021/src/MOLECULE/ |
H A D | dihedral_opls.cpp | 90 vb1z = x[i1][2] - x[i2][2]; 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() 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 83 vb1z = x[i1][2] - x[i2][2]; 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() 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 87 vb1z = x[i2][2] - x[i1][2]; in compute() 104 ax = vb1y*vb2z-vb1z*vb2y; in compute() 105 ay = vb1z*vb2x-vb1x*vb2z; in compute() 137 projhfg = (vb3x*vb1x+vb3y*vb1y+vb3z*vb1z) / in compute() 138 sqrt(vb1x*vb1x+vb1y*vb1y+vb1z*vb1z); in compute() 171 f2[0] = (dhay*vb1z - dhaz*vb1y)*rar*a; in compute() 172 f2[1] = (dhaz*vb1x - dhax*vb1z)*rar*a; in compute() 219 vb1z = x[i1][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 83 vb1z = x[i1][2] - x[i2][2]; 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() 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 90 vb1z = x[i1][2] - x[i2][2]; 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() 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/ |
H A D | compute_improper_local.cpp | 113 double vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z; in compute_impropers() local 180 vb1z = x[atom1][2] - x[atom2][2]; in compute_impropers() 181 domain->minimum_image(vb1x,vb1y,vb1z); in compute_impropers() 193 ss1 = 1.0 / (vb1x*vb1x + vb1y*vb1y + vb1z*vb1z); in compute_impropers() 201 c0 = (vb1x * vb3x + vb1y * vb3y + vb1z * vb3z) * r1 * r3; in compute_impropers() 202 c1 = (vb1x * vb2x + vb1y * vb2y + vb1z * vb2z) * r1 * r2; in compute_impropers()
|
/dports/science/lammps/lammps-stable_29Sep2021/src/KOKKOS/ |
H A D | dihedral_opls_kokkos.cpp | 174 const F_FLOAT vb1z = x(i1,2) - x(i2,2); 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 ()() 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 ()() 328 vb1x,vb1y,vb1z,vb2x,vb2y,vb2z,vb3x,vb3y,vb3z); in operator ()() [all …]
|