103 |
|
Mat3x3d tmpMat; |
104 |
|
tmpMat = outProduct(Rij, Rij) / rij2; |
105 |
|
RealType constant = 8.0 * NumericConstant::PI * viscosity * rij; |
106 |
< |
Tij = ((1.0 + sumSigma2OverRij2/3.0) * I + (1.0 - sumSigma2OverRij2) * tmpMat ) / constant; |
106 |
> |
RealType tmp1 = 1.0 + sumSigma2OverRij2/3.0; |
107 |
> |
RealType tmp2 = 1.0 - sumSigma2OverRij2; |
108 |
> |
Tij = (tmp1 * I + tmp2 * tmpMat ) / constant; |
109 |
|
}else { |
110 |
|
RealType constant = 1.0 / (6.0 * NumericConstant::PI * viscosity * beads[i].radius); |
111 |
|
Tij(0, 0) = constant; |
269 |
|
Mat3x3d tmpMat; |
270 |
|
tmpMat = outProduct(Rij, Rij) / rij2; |
271 |
|
RealType constant = 8.0 * NumericConstant::PI * viscosity * rij; |
272 |
< |
Tij = ((1.0 + sumSigma2OverRij2/3.0) * I + (1.0 - sumSigma2OverRij2) * tmpMat ) / constant; |
272 |
> |
RealType tmp1 = 1.0 + sumSigma2OverRij2/3.0; |
273 |
> |
RealType tmp2 = 1.0 - sumSigma2OverRij2; |
274 |
> |
Tij = (tmp1 * I + tmp2 * tmpMat ) / constant; |
275 |
|
}else { |
276 |
|
RealType constant = 1.0 / (6.0 * NumericConstant::PI * viscosity * beads[i].radius); |
277 |
|
Tij(0, 0) = constant; |