| 781 | |
| 782 | |
| 783 | static void VLin1(real a, N_Vector x, N_Vector y, N_Vector z) |
| 784 | { |
| 785 | integer N; |
| 786 | real *xd, *yd, *zd; |
| 787 | |
| 788 | N = x->length; |
| 789 | xd = x->data; |
| 790 | yd = y->data; |
| 791 | zd = z->data; |
| 792 | |
| 793 | #pragma omp parallel for |
| 794 | for (integer i=0; i < N; i++) |
| 795 | zd[i] = a * xd[i] + yd[i]; |
| 796 | } |
| 797 | |
| 798 | |
| 799 | static void VLin2(real a, N_Vector x, N_Vector y, N_Vector z) |