1 /* blas/source_dot_r.h 2 * 3 * Copyright (C) 1996, 1997, 1998, 1999, 2000 Gerard Jungman 4 * 5 * This program is free software; you can redistribute it and/or modify 6 * it under the terms of the GNU General Public License as published by 7 * the Free Software Foundation; either version 3 of the License, or (at 8 * your option) any later version. 9 * 10 * This program is distributed in the hope that it will be useful, but 11 * WITHOUT ANY WARRANTY; without even the implied warranty of 12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 * General Public License for more details. 14 * 15 * You should have received a copy of the GNU General Public License 16 * along with this program; if not, write to the Free Software 17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. 18 */ 19 20 { 21 ACC_TYPE r = INIT_VAL; 22 INDEX i; 23 INDEX ix = OFFSET(N, incX); 24 INDEX iy = OFFSET(N, incY); 25 26 for (i = 0; i < N; i++) { 27 r += X[ix] * Y[iy]; 28 ix += incX; 29 iy += incY; 30 } 31 32 return r; 33 } 34