Lines Matching +full:- +full:r
1 /* zhpmv.f -- translated by f2c (version 20100827).
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
41 /* ZHPMV performs the matrix-vector operation */ in zhpmv_()
51 /* UPLO - CHARACTER*1. */ in zhpmv_()
64 /* N - INTEGER. */ in zhpmv_()
69 /* ALPHA - COMPLEX*16 . */ in zhpmv_()
73 /* AP - COMPLEX*16 array of DIMENSION at least */ in zhpmv_()
89 /* X - COMPLEX*16 array of dimension at least */ in zhpmv_()
90 /* ( 1 + ( n - 1 )*abs( INCX ) ). */ in zhpmv_()
95 /* INCX - INTEGER. */ in zhpmv_()
100 /* BETA - COMPLEX*16 . */ in zhpmv_()
105 /* Y - COMPLEX*16 array of dimension at least */ in zhpmv_()
106 /* ( 1 + ( n - 1 )*abs( INCY ) ). */ in zhpmv_()
111 /* INCY - INTEGER. */ in zhpmv_()
121 /* -- Written on 22-October-1986. */ in zhpmv_()
143 --y; in zhpmv_()
144 --x; in zhpmv_()
145 --ap; in zhpmv_()
166 if (*n == 0 || (alpha->r == 0. && alpha->i == 0. && (beta->r == 1. && in zhpmv_()
167 beta->i == 0.))) { in zhpmv_()
176 kx = 1 - (*n - 1) * *incx; in zhpmv_()
181 ky = 1 - (*n - 1) * *incy; in zhpmv_()
189 if (beta->r != 1. || beta->i != 0.) { in zhpmv_()
191 if (beta->r == 0. && beta->i == 0.) { in zhpmv_()
195 y[i__2].r = 0., y[i__2].i = 0.; in zhpmv_()
203 z__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in zhpmv_()
204 z__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in zhpmv_()
205 .r; in zhpmv_()
206 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
212 if (beta->r == 0. && beta->i == 0.) { in zhpmv_()
216 y[i__2].r = 0., y[i__2].i = 0.; in zhpmv_()
225 z__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in zhpmv_()
226 z__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in zhpmv_()
227 .r; in zhpmv_()
228 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
235 if (alpha->r == 0. && alpha->i == 0.) { in zhpmv_()
247 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
248 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
249 temp1.r = z__1.r, temp1.i = z__1.i; in zhpmv_()
250 temp2.r = 0., temp2.i = 0.; in zhpmv_()
252 i__2 = j - 1; in zhpmv_()
257 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
258 z__2.i = temp1.r * ap[i__5].i + temp1.i * ap[i__5] in zhpmv_()
259 .r; in zhpmv_()
260 z__1.r = y[i__4].r + z__2.r, z__1.i = y[i__4].i + z__2.i; in zhpmv_()
261 y[i__3].r = z__1.r, y[i__3].i = z__1.i; in zhpmv_()
264 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
265 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
266 z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i; in zhpmv_()
267 temp2.r = z__1.r, temp2.i = z__1.i; in zhpmv_()
273 i__4 = kk + j - 1; in zhpmv_()
274 d__1 = ap[i__4].r; in zhpmv_()
275 z__3.r = d__1 * temp1.r, z__3.i = d__1 * temp1.i; in zhpmv_()
276 z__2.r = y[i__3].r + z__3.r, z__2.i = y[i__3].i + z__3.i; in zhpmv_()
277 z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i = in zhpmv_()
278 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
279 z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i; in zhpmv_()
280 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
290 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
291 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
292 temp1.r = z__1.r, temp1.i = z__1.i; in zhpmv_()
293 temp2.r = 0., temp2.i = 0.; in zhpmv_()
296 i__2 = kk + j - 2; in zhpmv_()
301 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
302 z__2.i = temp1.r * ap[i__5].i + temp1.i * ap[i__5] in zhpmv_()
303 .r; in zhpmv_()
304 z__1.r = y[i__4].r + z__2.r, z__1.i = y[i__4].i + z__2.i; in zhpmv_()
305 y[i__3].r = z__1.r, y[i__3].i = z__1.i; in zhpmv_()
308 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
309 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
310 z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i; in zhpmv_()
311 temp2.r = z__1.r, temp2.i = z__1.i; in zhpmv_()
318 i__4 = kk + j - 1; in zhpmv_()
319 d__1 = ap[i__4].r; in zhpmv_()
320 z__3.r = d__1 * temp1.r, z__3.i = d__1 * temp1.i; in zhpmv_()
321 z__2.r = y[i__3].r + z__3.r, z__2.i = y[i__3].i + z__3.i; in zhpmv_()
322 z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i = in zhpmv_()
323 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
324 z__1.r = z__2.r + z__4.r, z__1.i = z__2.i + z__4.i; in zhpmv_()
325 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
340 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
341 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
342 temp1.r = z__1.r, temp1.i = z__1.i; in zhpmv_()
343 temp2.r = 0., temp2.i = 0.; in zhpmv_()
347 d__1 = ap[i__4].r; in zhpmv_()
348 z__2.r = d__1 * temp1.r, z__2.i = d__1 * temp1.i; in zhpmv_()
349 z__1.r = y[i__3].r + z__2.r, z__1.i = y[i__3].i + z__2.i; in zhpmv_()
350 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
357 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
358 z__2.i = temp1.r * ap[i__5].i + temp1.i * ap[i__5] in zhpmv_()
359 .r; in zhpmv_()
360 z__1.r = y[i__4].r + z__2.r, z__1.i = y[i__4].i + z__2.i; in zhpmv_()
361 y[i__3].r = z__1.r, y[i__3].i = z__1.i; in zhpmv_()
364 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
365 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
366 z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i; in zhpmv_()
367 temp2.r = z__1.r, temp2.i = z__1.i; in zhpmv_()
373 z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i = in zhpmv_()
374 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
375 z__1.r = y[i__3].r + z__2.r, z__1.i = y[i__3].i + z__2.i; in zhpmv_()
376 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
377 kk += *n - j + 1; in zhpmv_()
386 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
387 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
388 temp1.r = z__1.r, temp1.i = z__1.i; in zhpmv_()
389 temp2.r = 0., temp2.i = 0.; in zhpmv_()
393 d__1 = ap[i__4].r; in zhpmv_()
394 z__2.r = d__1 * temp1.r, z__2.i = d__1 * temp1.i; in zhpmv_()
395 z__1.r = y[i__3].r + z__2.r, z__1.i = y[i__3].i + z__2.i; in zhpmv_()
396 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
399 i__2 = kk + *n - j; in zhpmv_()
406 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
407 z__2.i = temp1.r * ap[i__5].i + temp1.i * ap[i__5] in zhpmv_()
408 .r; in zhpmv_()
409 z__1.r = y[i__4].r + z__2.r, z__1.i = y[i__4].i + z__2.i; in zhpmv_()
410 y[i__3].r = z__1.r, y[i__3].i = z__1.i; in zhpmv_()
413 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
414 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
415 z__1.r = temp2.r + z__2.r, z__1.i = temp2.i + z__2.i; in zhpmv_()
416 temp2.r = z__1.r, temp2.i = z__1.i; in zhpmv_()
421 z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i = in zhpmv_()
422 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
423 z__1.r = y[i__3].r + z__2.r, z__1.i = y[i__3].i + z__2.i; in zhpmv_()
424 y[i__2].r = z__1.r, y[i__2].i = z__1.i; in zhpmv_()
427 kk += *n - j + 1; in zhpmv_()