Lines Matching refs:j1
29 static integer j1, j2, n1, n2; in dlaqtr_() local
247 j1 = j; in dlaqtr_()
252 j1 = j - 1; in dlaqtr_()
257 if (j1 == j2) { in dlaqtr_()
264 xj = (d__1 = x[j1], abs(d__1)); in dlaqtr_()
265 tjj = (d__1 = t[j1 + j1 * t_dim1], abs(d__1)); in dlaqtr_()
266 tmp = t[j1 + j1 * t_dim1]; in dlaqtr_()
285 x[j1] /= tmp; in dlaqtr_()
286 xj = (d__1 = x[j1], abs(d__1)); in dlaqtr_()
293 if (work[j1] > (bignum - xmax) * rec) { in dlaqtr_()
298 if (j1 > 1) { in dlaqtr_()
299 i__1 = j1 - 1; in dlaqtr_()
300 d__1 = -x[j1]; in dlaqtr_()
301 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in dlaqtr_()
303 i__1 = j1 - 1; in dlaqtr_()
315 d__[0] = x[j1]; in dlaqtr_()
317 dlaln2_(&c_false, &c__2, &c__1, &smin, &c_b21, &t[j1 + j1 in dlaqtr_()
328 x[j1] = v[0]; in dlaqtr_()
340 d__1 = work[j1], d__2 = work[j2]; in dlaqtr_()
349 if (j1 > 1) { in dlaqtr_()
350 i__1 = j1 - 1; in dlaqtr_()
351 d__1 = -x[j1]; in dlaqtr_()
352 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in dlaqtr_()
354 i__1 = j1 - 1; in dlaqtr_()
358 i__1 = j1 - 1; in dlaqtr_()
379 j1 = j; in dlaqtr_()
389 if (j1 == j2) { in dlaqtr_()
396 xj = (d__1 = x[j1], abs(d__1)); in dlaqtr_()
399 if (work[j1] > (bignum - xj) * rec) { in dlaqtr_()
406 i__2 = j1 - 1; in dlaqtr_()
407 x[j1] -= ddot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[1], & in dlaqtr_()
410 xj = (d__1 = x[j1], abs(d__1)); in dlaqtr_()
411 tjj = (d__1 = t[j1 + j1 * t_dim1], abs(d__1)); in dlaqtr_()
412 tmp = t[j1 + j1 * t_dim1]; in dlaqtr_()
427 x[j1] /= tmp; in dlaqtr_()
429 d__2 = xmax, d__3 = (d__1 = x[j1], abs(d__1)); in dlaqtr_()
440 d__3 = (d__1 = x[j1], abs(d__1)), d__4 = (d__2 = x[j2], in dlaqtr_()
446 d__1 = work[j2], d__2 = work[j1]; in dlaqtr_()
454 i__2 = j1 - 1; in dlaqtr_()
455 d__[0] = x[j1] - ddot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, in dlaqtr_()
457 i__2 = j1 - 1; in dlaqtr_()
461 dlaln2_(&c_true, &c__2, &c__1, &smin, &c_b21, &t[j1 + j1 * in dlaqtr_()
472 x[j1] = v[0]; in dlaqtr_()
475 d__3 = (d__1 = x[j1], abs(d__1)), d__4 = (d__2 = x[j2], in dlaqtr_()
499 j1 = j; in dlaqtr_()
504 j1 = j - 1; in dlaqtr_()
509 if (j1 == j2) { in dlaqtr_()
516 if (j1 == 1) { in dlaqtr_()
519 xj = (d__1 = x[j1], abs(d__1)) + (d__2 = x[*n + j1], abs( in dlaqtr_()
521 tjj = (d__1 = t[j1 + j1 * t_dim1], abs(d__1)) + abs(z__); in dlaqtr_()
522 tmp = t[j1 + j1 * t_dim1]; in dlaqtr_()
541 dladiv_(&x[j1], &x[*n + j1], &tmp, &z__, &sr, &si); in dlaqtr_()
542 x[j1] = sr; in dlaqtr_()
543 x[*n + j1] = si; in dlaqtr_()
544 xj = (d__1 = x[j1], abs(d__1)) + (d__2 = x[*n + j1], abs( in dlaqtr_()
552 if (work[j1] > (bignum - xmax) * rec) { in dlaqtr_()
558 if (j1 > 1) { in dlaqtr_()
559 i__1 = j1 - 1; in dlaqtr_()
560 d__1 = -x[j1]; in dlaqtr_()
561 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in dlaqtr_()
563 i__1 = j1 - 1; in dlaqtr_()
564 d__1 = -x[*n + j1]; in dlaqtr_()
565 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[* in dlaqtr_()
568 x[1] += b[j1] * x[*n + j1]; in dlaqtr_()
569 x[*n + 1] -= b[j1] * x[j1]; in dlaqtr_()
572 i__1 = j1 - 1; in dlaqtr_()
586 d__[0] = x[j1]; in dlaqtr_()
588 d__[2] = x[*n + j1]; in dlaqtr_()
591 dlaln2_(&c_false, &c__2, &c__2, &sminw, &c_b21, &t[j1 + in dlaqtr_()
592 j1 * t_dim1], ldt, &c_b21, &c_b21, d__, &c__2, & in dlaqtr_()
603 x[j1] = v[0]; in dlaqtr_()
605 x[*n + j1] = v[2]; in dlaqtr_()
618 d__1 = work[j1], d__2 = work[j2]; in dlaqtr_()
627 if (j1 > 1) { in dlaqtr_()
628 i__1 = j1 - 1; in dlaqtr_()
629 d__1 = -x[j1]; in dlaqtr_()
630 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in dlaqtr_()
632 i__1 = j1 - 1; in dlaqtr_()
637 i__1 = j1 - 1; in dlaqtr_()
638 d__1 = -x[*n + j1]; in dlaqtr_()
639 daxpy_(&i__1, &d__1, &t[j1 * t_dim1 + 1], &c__1, &x[* in dlaqtr_()
641 i__1 = j1 - 1; in dlaqtr_()
646 x[1] = x[1] + b[j1] * x[*n + j1] + b[j2] * x[*n + j2]; in dlaqtr_()
647 x[*n + 1] = x[*n + 1] - b[j1] * x[j1] - b[j2] * x[j2]; in dlaqtr_()
650 i__1 = j1 - 1; in dlaqtr_()
675 j1 = j; in dlaqtr_()
685 if (j1 == j2) { in dlaqtr_()
692 xj = (d__1 = x[j1], abs(d__1)) + (d__2 = x[j1 + *n], abs( in dlaqtr_()
696 if (work[j1] > (bignum - xj) * rec) { in dlaqtr_()
703 i__2 = j1 - 1; in dlaqtr_()
704 x[j1] -= ddot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[1], & in dlaqtr_()
706 i__2 = j1 - 1; in dlaqtr_()
707 x[*n + j1] -= ddot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[ in dlaqtr_()
709 if (j1 > 1) { in dlaqtr_()
710 x[j1] -= b[j1] * x[*n + 1]; in dlaqtr_()
711 x[*n + j1] += b[j1] * x[1]; in dlaqtr_()
713 xj = (d__1 = x[j1], abs(d__1)) + (d__2 = x[j1 + *n], abs( in dlaqtr_()
717 if (j1 == 1) { in dlaqtr_()
724 tjj = (d__1 = t[j1 + j1 * t_dim1], abs(d__1)) + abs(z__); in dlaqtr_()
725 tmp = t[j1 + j1 * t_dim1]; in dlaqtr_()
741 dladiv_(&x[j1], &x[*n + j1], &tmp, &d__1, &sr, &si); in dlaqtr_()
742 x[j1] = sr; in dlaqtr_()
743 x[j1 + *n] = si; in dlaqtr_()
745 d__3 = (d__1 = x[j1], abs(d__1)) + (d__2 = x[j1 + *n], in dlaqtr_()
757 d__5 = (d__1 = x[j1], abs(d__1)) + (d__2 = x[*n + j1], in dlaqtr_()
764 d__1 = work[j1], d__2 = work[j2]; in dlaqtr_()
772 i__2 = j1 - 1; in dlaqtr_()
773 d__[0] = x[j1] - ddot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, in dlaqtr_()
775 i__2 = j1 - 1; in dlaqtr_()
778 i__2 = j1 - 1; in dlaqtr_()
779 d__[2] = x[*n + j1] - ddot_(&i__2, &t[j1 * t_dim1 + 1], & in dlaqtr_()
781 i__2 = j1 - 1; in dlaqtr_()
784 d__[0] -= b[j1] * x[*n + 1]; in dlaqtr_()
786 d__[2] += b[j1] * x[1]; in dlaqtr_()
789 dlaln2_(&c_true, &c__2, &c__2, &sminw, &c_b21, &t[j1 + j1 in dlaqtr_()
800 x[j1] = v[0]; in dlaqtr_()
802 x[*n + j1] = v[2]; in dlaqtr_()
805 d__5 = (d__1 = x[j1], abs(d__1)) + (d__2 = x[*n + j1], in dlaqtr_()