Lines Matching refs:j1
29 static integer j1, j2, n1, n2; in slaqtr_() local
244 j1 = j; in slaqtr_()
249 j1 = j - 1; in slaqtr_()
254 if (j1 == j2) { in slaqtr_()
261 xj = (r__1 = x[j1], dabs(r__1)); in slaqtr_()
262 tjj = (r__1 = t[j1 + j1 * t_dim1], dabs(r__1)); in slaqtr_()
263 tmp = t[j1 + j1 * t_dim1]; in slaqtr_()
282 x[j1] /= tmp; in slaqtr_()
283 xj = (r__1 = x[j1], dabs(r__1)); in slaqtr_()
290 if (work[j1] > (bignum - xmax) * rec) { in slaqtr_()
295 if (j1 > 1) { in slaqtr_()
296 i__1 = j1 - 1; in slaqtr_()
297 r__1 = -x[j1]; in slaqtr_()
298 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in slaqtr_()
300 i__1 = j1 - 1; in slaqtr_()
312 d__[0] = x[j1]; in slaqtr_()
314 slaln2_(&c_false, &c__2, &c__1, &smin, &c_b21, &t[j1 + j1 in slaqtr_()
325 x[j1] = v[0]; in slaqtr_()
337 r__1 = work[j1], r__2 = work[j2]; in slaqtr_()
346 if (j1 > 1) { in slaqtr_()
347 i__1 = j1 - 1; in slaqtr_()
348 r__1 = -x[j1]; in slaqtr_()
349 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in slaqtr_()
351 i__1 = j1 - 1; in slaqtr_()
355 i__1 = j1 - 1; in slaqtr_()
376 j1 = j; in slaqtr_()
386 if (j1 == j2) { in slaqtr_()
393 xj = (r__1 = x[j1], dabs(r__1)); in slaqtr_()
396 if (work[j1] > (bignum - xj) * rec) { in slaqtr_()
403 i__2 = j1 - 1; in slaqtr_()
404 x[j1] -= sdot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[1], & in slaqtr_()
407 xj = (r__1 = x[j1], dabs(r__1)); in slaqtr_()
408 tjj = (r__1 = t[j1 + j1 * t_dim1], dabs(r__1)); in slaqtr_()
409 tmp = t[j1 + j1 * t_dim1]; in slaqtr_()
424 x[j1] /= tmp; in slaqtr_()
426 r__2 = xmax, r__3 = (r__1 = x[j1], dabs(r__1)); in slaqtr_()
437 r__3 = (r__1 = x[j1], dabs(r__1)), r__4 = (r__2 = x[j2], in slaqtr_()
443 r__1 = work[j2], r__2 = work[j1]; in slaqtr_()
451 i__2 = j1 - 1; in slaqtr_()
452 d__[0] = x[j1] - sdot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, in slaqtr_()
454 i__2 = j1 - 1; in slaqtr_()
458 slaln2_(&c_true, &c__2, &c__1, &smin, &c_b21, &t[j1 + j1 * in slaqtr_()
469 x[j1] = v[0]; in slaqtr_()
472 r__3 = (r__1 = x[j1], dabs(r__1)), r__4 = (r__2 = x[j2], in slaqtr_()
496 j1 = j; in slaqtr_()
501 j1 = j - 1; in slaqtr_()
506 if (j1 == j2) { in slaqtr_()
513 if (j1 == 1) { in slaqtr_()
516 xj = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[*n + j1], in slaqtr_()
518 tjj = (r__1 = t[j1 + j1 * t_dim1], dabs(r__1)) + dabs(z__) in slaqtr_()
520 tmp = t[j1 + j1 * t_dim1]; in slaqtr_()
539 sladiv_(&x[j1], &x[*n + j1], &tmp, &z__, &sr, &si); in slaqtr_()
540 x[j1] = sr; in slaqtr_()
541 x[*n + j1] = si; in slaqtr_()
542 xj = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[*n + j1], in slaqtr_()
550 if (work[j1] > (bignum - xmax) * rec) { in slaqtr_()
556 if (j1 > 1) { in slaqtr_()
557 i__1 = j1 - 1; in slaqtr_()
558 r__1 = -x[j1]; in slaqtr_()
559 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in slaqtr_()
561 i__1 = j1 - 1; in slaqtr_()
562 r__1 = -x[*n + j1]; in slaqtr_()
563 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[* in slaqtr_()
566 x[1] += b[j1] * x[*n + j1]; in slaqtr_()
567 x[*n + 1] -= b[j1] * x[j1]; in slaqtr_()
570 i__1 = j1 - 1; in slaqtr_()
584 d__[0] = x[j1]; in slaqtr_()
586 d__[2] = x[*n + j1]; in slaqtr_()
589 slaln2_(&c_false, &c__2, &c__2, &sminw, &c_b21, &t[j1 + in slaqtr_()
590 j1 * t_dim1], ldt, &c_b21, &c_b21, d__, &c__2, & in slaqtr_()
601 x[j1] = v[0]; in slaqtr_()
603 x[*n + j1] = v[2]; in slaqtr_()
616 r__1 = work[j1], r__2 = work[j2]; in slaqtr_()
625 if (j1 > 1) { in slaqtr_()
626 i__1 = j1 - 1; in slaqtr_()
627 r__1 = -x[j1]; in slaqtr_()
628 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[1] in slaqtr_()
630 i__1 = j1 - 1; in slaqtr_()
635 i__1 = j1 - 1; in slaqtr_()
636 r__1 = -x[*n + j1]; in slaqtr_()
637 saxpy_(&i__1, &r__1, &t[j1 * t_dim1 + 1], &c__1, &x[* in slaqtr_()
639 i__1 = j1 - 1; in slaqtr_()
644 x[1] = x[1] + b[j1] * x[*n + j1] + b[j2] * x[*n + j2]; in slaqtr_()
645 x[*n + 1] = x[*n + 1] - b[j1] * x[j1] - b[j2] * x[j2]; in slaqtr_()
648 i__1 = j1 - 1; in slaqtr_()
673 j1 = j; in slaqtr_()
683 if (j1 == j2) { in slaqtr_()
690 xj = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[j1 + *n], in slaqtr_()
694 if (work[j1] > (bignum - xj) * rec) { in slaqtr_()
701 i__2 = j1 - 1; in slaqtr_()
702 x[j1] -= sdot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[1], & in slaqtr_()
704 i__2 = j1 - 1; in slaqtr_()
705 x[*n + j1] -= sdot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, &x[ in slaqtr_()
707 if (j1 > 1) { in slaqtr_()
708 x[j1] -= b[j1] * x[*n + 1]; in slaqtr_()
709 x[*n + j1] += b[j1] * x[1]; in slaqtr_()
711 xj = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[j1 + *n], in slaqtr_()
715 if (j1 == 1) { in slaqtr_()
722 tjj = (r__1 = t[j1 + j1 * t_dim1], dabs(r__1)) + dabs(z__) in slaqtr_()
724 tmp = t[j1 + j1 * t_dim1]; in slaqtr_()
740 sladiv_(&x[j1], &x[*n + j1], &tmp, &r__1, &sr, &si); in slaqtr_()
741 x[j1] = sr; in slaqtr_()
742 x[j1 + *n] = si; in slaqtr_()
744 r__3 = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[j1 + *n], in slaqtr_()
756 r__5 = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[*n + j1], in slaqtr_()
763 r__1 = work[j1], r__2 = work[j2]; in slaqtr_()
771 i__2 = j1 - 1; in slaqtr_()
772 d__[0] = x[j1] - sdot_(&i__2, &t[j1 * t_dim1 + 1], &c__1, in slaqtr_()
774 i__2 = j1 - 1; in slaqtr_()
777 i__2 = j1 - 1; in slaqtr_()
778 d__[2] = x[*n + j1] - sdot_(&i__2, &t[j1 * t_dim1 + 1], & in slaqtr_()
780 i__2 = j1 - 1; in slaqtr_()
783 d__[0] -= b[j1] * x[*n + 1]; in slaqtr_()
785 d__[2] += b[j1] * x[1]; in slaqtr_()
788 slaln2_(&c_true, &c__2, &c__2, &sminw, &c_b21, &t[j1 + j1 in slaqtr_()
799 x[j1] = v[0]; in slaqtr_()
801 x[*n + j1] = v[2]; in slaqtr_()
804 r__5 = (r__1 = x[j1], dabs(r__1)) + (r__2 = x[*n + j1], in slaqtr_()