| 149 | |
| 150 | |
| 151 | double ScalN (int n, const double * pa, const double * pb) |
| 152 | { |
| 153 | double sum = 0; |
| 154 | int i; |
| 155 | for (i = 0; i < n; i++) |
| 156 | sum += pa[i] * pb[i]; |
| 157 | return sum; |
| 158 | } |
| 159 | |
| 160 | template <int N> |
| 161 | double Scal (const double * pa, const double * pb) |