| 797 | |
| 798 | |
| 799 | static void VLin2(real a, N_Vector x, N_Vector y, N_Vector z) |
| 800 | { |
| 801 | integer N; |
| 802 | real *xd, *yd, *zd; |
| 803 | |
| 804 | N = x->length; |
| 805 | xd = x->data; |
| 806 | yd = y->data; |
| 807 | zd = z->data; |
| 808 | |
| 809 | #pragma omp parallel for |
| 810 | for (integer i=0; i < N; i++) |
| 811 | zd[i] = a * xd[i] - yd[i]; |
| 812 | } |
| 813 | |
| 814 | static void Vaxpy(real a, N_Vector x, N_Vector y) |
| 815 | { |