| 116 |
|
} |
| 117 |
|
} |
| 118 |
|
|
| 119 |
< |
template<typename T> void NPTxyz<T>::getVelScaleA(double sc[3], double vel[3]) { |
| 119 |
> |
template<typename T> void NPTxyz<T>::calcVelScale(void) { |
| 120 |
|
int i,j; |
| 121 |
– |
double vScale[3][3]; |
| 121 |
|
|
| 122 |
|
for (i = 0; i < 3; i++ ) { |
| 123 |
|
for (j = 0; j < 3; j++ ) { |
| 128 |
|
} |
| 129 |
|
} |
| 130 |
|
} |
| 131 |
+ |
} |
| 132 |
|
|
| 133 |
+ |
template<typename T> void NPTxyz<T>::getVelScaleA(double sc[3], double vel[3]) { |
| 134 |
|
info->matVecMul3( vScale, vel, sc ); |
| 135 |
|
} |
| 136 |
|
|
| 137 |
|
template<typename T> void NPTxyz<T>::getVelScaleB(double sc[3], int index ){ |
| 138 |
< |
int i,j; |
| 138 |
> |
int j; |
| 139 |
|
double myVel[3]; |
| 139 |
– |
double vScale[3][3]; |
| 140 |
– |
|
| 141 |
– |
for (i = 0; i < 3; i++ ) { |
| 142 |
– |
for (j = 0; j < 3; j++ ) { |
| 143 |
– |
vScale[i][j] = eta[i][j]; |
| 140 |
|
|
| 145 |
– |
if (i == j) { |
| 146 |
– |
vScale[i][j] += chi; |
| 147 |
– |
} |
| 148 |
– |
} |
| 149 |
– |
} |
| 150 |
– |
|
| 141 |
|
for (j = 0; j < 3; j++) |
| 142 |
|
myVel[j] = oldVel[3*index + j]; |
| 143 |
|
|