| 1477 | return oOr; |
| 1478 | } |
| 1479 | inline double* getrw2(double* r, int rsize) |
| 1480 | { |
| 1481 | double* rw2 = (double*) malloc(sd*rsize); |
| 1482 | for(int w = 0; w < rsize; w++){ |
| 1483 | rw2[w] = r[w]*r[w]; |
| 1484 | } |
| 1485 | return rw2; |
| 1486 | } |
| 1487 | inline void expMs(double* rExpDiff, double eta, double* r, double* ri, int isize, int rsize) |
| 1488 | { |
| 1489 | double rDiff; |