| 1485 | return rw2; |
| 1486 | } |
| 1487 | inline void expMs(double* rExpDiff, double eta, double* r, double* ri, int isize, int rsize) |
| 1488 | { |
| 1489 | double rDiff; |
| 1490 | for (int i = 0; i < isize; i++) { |
| 1491 | for (int w = 0; w < rsize; w++) { |
| 1492 | rDiff = r[w] - ri[i]; |
| 1493 | if (rDiff > 5.0 ) { |
| 1494 | rExpDiff[rsize*i + w] = 0.0; |
| 1495 | } else { |
| 1496 | rExpDiff[rsize*i + w] = exp(-eta*rDiff*rDiff); |
| 1497 | } |
| 1498 | } |
| 1499 | } |
| 1500 | } |
| 1501 | inline void expPs(double* rExpSum, double eta, double* r, double* ri, int isize, int rsize) |
| 1502 | { |
| 1503 | double rSum; |