/*
 * Author:  G. Jungman
 * RCS:     $Id: source_spr.h,v 1.4 1999/10/18 00:45:12 jungman Exp $
 */

  size_t j, k;
  size_t ix, jx;
  size_t kk = 0;

  if(Uplo == CblasUpper) {
    jx = 0;
    for(j=0; j<N; j++) {
      BASE_TYPE tmp = alpha * X[jx];
      ix = jx;
      for(k=kk; k<kk+N-j; k++) {
        Ap[k] += X[ix]*tmp;
	ix += incX;
      }
      jx += incX;
      kk += N - j;
    }
  }
  else {
    jx = 0;
    for(j=0; j<N; j++) {
      BASE_TYPE tmp = alpha * X[jx];
      ix = 0;
      for(k=kk; k<kk+j+1; k++) {
        Ap[k] += X[ix] * tmp;
        ix += incX;
      }
      jx += incX;
      kk += j+1;
    }
  }