| 749 | |
| 750 | |
| 751 | static void VScaleSum(real c, N_Vector x, N_Vector y, N_Vector z) |
| 752 | { |
| 753 | integer N; |
| 754 | real *xd, *yd, *zd; |
| 755 | |
| 756 | N = x->length; |
| 757 | xd = x->data; |
| 758 | yd = y->data; |
| 759 | zd = z->data; |
| 760 | |
| 761 | #pragma omp parallel for |
| 762 | for (integer i=0; i < N; i++) |
| 763 | zd[i] = c * (xd[i] + yd[i]); |
| 764 | } |
| 765 | |
| 766 | |
| 767 | static void VScaleDiff(real c, N_Vector x, N_Vector y, N_Vector z) |