61a62,64
>   double* hptr;
>   double* gptr;
>   double h, g;
80c83,85
<       for (ii = 0, i = 0; i < n; i += 2, ii++)
---
>       hptr = work->scratch;
>       gptr = &(work->scratch[nh]);
>       for (i = 0; i < n; i += 2)
82a88,89
> 	  h = *hptr;
> 	  g = *gptr;
86,87c93,94
<               work->scratch[ii] += w->h1[k] * ELEMENT (a, stride, jf);
<               work->scratch[ii + nh] += w->g1[k] * ELEMENT (a, stride, jf);
---
>               h += w->h1[k] * ELEMENT (a, stride, jf);
>               g += w->g1[k] * ELEMENT (a, stride, jf);
88a96,97
> 	  *hptr++ = h;
> 	  *gptr++ = g;
97a107
>