| 718 | |
| 719 | |
| 720 | static void VDiff(N_Vector x, N_Vector y, N_Vector z) |
| 721 | { |
| 722 | integer N; |
| 723 | real *xd, *yd, *zd; |
| 724 | |
| 725 | N = x->length; |
| 726 | xd = x->data; |
| 727 | yd = y->data; |
| 728 | zd = z->data; |
| 729 | |
| 730 | #pragma omp parallel for |
| 731 | for (integer i=0; i < N; i++) |
| 732 | zd[i] = xd[i] - yd[i]; |
| 733 | } |
| 734 | |
| 735 | |
| 736 | static void VNeg(N_Vector x, N_Vector z) |