| 702 | |
| 703 | |
| 704 | static void VSum(N_Vector x, N_Vector y, N_Vector z) |
| 705 | { |
| 706 | integer N; |
| 707 | real *xd, *yd, *zd; |
| 708 | |
| 709 | N = x->length; |
| 710 | xd = x->data; |
| 711 | yd = y->data; |
| 712 | zd = z->data; |
| 713 | |
| 714 | #pragma omp parallel for |
| 715 | for (integer i=0; i < N; i++) |
| 716 | zd[i] = xd[i] + yd[i]; |
| 717 | } |
| 718 | |
| 719 | |
| 720 | static void VDiff(N_Vector x, N_Vector y, N_Vector z) |