Lines Matching +refs:eta +refs:temp
34 doublereal eta, phi, eps, tau, psi; in dlasd4_() local
38 doublereal temp, prew, sg2lb, sg2ub, temp1, temp2, dtiim, delsq, dtiip; in dlasd4_() local
200 temp = *rho / 2.; in dlasd4_()
205 temp1 = temp / (d__[*n] + sqrt(d__[*n] * d__[*n] + temp)); in dlasd4_()
226 temp = z__[*n - 1] * z__[*n - 1] / ((d__[*n - 1] + temp1) * (d__[* in dlasd4_()
233 if (c__ <= temp) { in dlasd4_()
271 eta = tau / (d__[*n] + sqrt(d__[*n] * d__[*n] + tau)); in dlasd4_()
273 *sigma = d__[*n] + eta; in dlasd4_()
276 delta[j] = d__[j] - d__[*i__] - eta; in dlasd4_()
277 work[j] = d__[j] + d__[*i__] + eta; in dlasd4_()
288 temp = z__[j] / (delta[j] * work[j]); in dlasd4_()
289 psi += z__[j] * temp; in dlasd4_()
290 dpsi += temp * temp; in dlasd4_()
298 temp = z__[*n] / (delta[*n] * work[*n]); in dlasd4_()
299 phi = z__[*n] * temp; in dlasd4_()
300 dphi = temp * temp; in dlasd4_()
324 eta = *rho - *sigma * *sigma; in dlasd4_()
326 eta = (a + sqrt((d__1 = a * a - b * 4. * c__, abs(d__1)))) / (c__ in dlasd4_()
329 eta = b * 2. / (a - sqrt((d__1 = a * a - b * 4. * c__, abs(d__1))) in dlasd4_()
339 if (w * eta > 0.) { in dlasd4_()
340 eta = -w / (dpsi + dphi); in dlasd4_()
342 temp = eta - dtnsq; in dlasd4_()
343 if (temp > *rho) { in dlasd4_()
344 eta = *rho + dtnsq; in dlasd4_()
347 tau += eta; in dlasd4_()
348 eta /= *sigma + sqrt(eta + *sigma * *sigma); in dlasd4_()
351 delta[j] -= eta; in dlasd4_()
352 work[j] += eta; in dlasd4_()
356 *sigma += eta; in dlasd4_()
365 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
366 psi += z__[j] * temp; in dlasd4_()
367 dpsi += temp * temp; in dlasd4_()
375 temp = z__[*n] / (work[*n] * delta[*n]); in dlasd4_()
376 phi = z__[*n] * temp; in dlasd4_()
377 dphi = temp * temp; in dlasd4_()
403 eta = (a + sqrt((d__1 = a * a - b * 4. * c__, abs(d__1)))) / ( in dlasd4_()
406 eta = b * 2. / (a - sqrt((d__1 = a * a - b * 4. * c__, abs( in dlasd4_()
416 if (w * eta > 0.) { in dlasd4_()
417 eta = -w / (dpsi + dphi); in dlasd4_()
419 temp = eta - dtnsq; in dlasd4_()
420 if (temp <= 0.) { in dlasd4_()
421 eta /= 2.; in dlasd4_()
424 tau += eta; in dlasd4_()
425 eta /= *sigma + sqrt(eta + *sigma * *sigma); in dlasd4_()
428 delta[j] -= eta; in dlasd4_()
429 work[j] += eta; in dlasd4_()
433 *sigma += eta; in dlasd4_()
442 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
443 psi += z__[j] * temp; in dlasd4_()
444 dpsi += temp * temp; in dlasd4_()
452 temp = z__[*n] / (work[*n] * delta[*n]); in dlasd4_()
453 phi = z__[*n] * temp; in dlasd4_()
454 dphi = temp * temp; in dlasd4_()
480 temp = delsq2 / (d__[*i__] + sqrt(d__[*i__] * d__[*i__] + delsq2)); in dlasd4_()
483 work[j] = d__[j] + d__[*i__] + temp; in dlasd4_()
484 delta[j] = d__[j] - d__[*i__] - temp; in dlasd4_()
528 eta = tau / (d__[*i__] + sqrt(d__[*i__] * d__[*i__] + tau)); in dlasd4_()
552 eta = tau / (d__[ip1] + sqrt((d__1 = d__[ip1] * d__[ip1] + tau, in dlasd4_()
558 *sigma = d__[*i__] + eta; in dlasd4_()
561 work[j] = d__[j] + d__[*i__] + eta; in dlasd4_()
562 delta[j] = d__[j] - d__[*i__] - eta; in dlasd4_()
567 *sigma = d__[ip1] + eta; in dlasd4_()
570 work[j] = d__[j] + d__[ip1] + eta; in dlasd4_()
571 delta[j] = d__[j] - d__[ip1] - eta; in dlasd4_()
585 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
586 psi += z__[j] * temp; in dlasd4_()
587 dpsi += temp * temp; in dlasd4_()
599 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
600 phi += z__[j] * temp; in dlasd4_()
601 dphi += temp * temp; in dlasd4_()
625 temp = z__[ii] / (work[ii] * delta[ii]); in dlasd4_()
626 dw = dpsi + dphi + temp * temp; in dlasd4_()
627 temp = z__[ii] * temp; in dlasd4_()
628 w += temp; in dlasd4_()
629 erretm = (phi - psi) * 8. + erretm + rhoinv * 2. + abs(temp) * 3. + in dlasd4_()
671 eta = b / a; in dlasd4_()
673 eta = (a - sqrt((d__1 = a * a - b * 4. * c__, abs(d__1)))) / ( in dlasd4_()
676 eta = b * 2. / (a + sqrt((d__1 = a * a - b * 4. * c__, abs( in dlasd4_()
685 temp = rhoinv + psi + phi; in dlasd4_()
689 c__ = temp - dtiip * (dpsi + dphi) - (d__[iim1] - d__[iip1]) * in dlasd4_()
700 c__ = temp - dtiim * (dpsi + dphi) - (d__[iip1] - d__[iim1]) * in dlasd4_()
713 dlaed6_(&niter, &orgati, &c__, dd, zz, &w, &eta, info); in dlasd4_()
725 if (w * eta >= 0.) { in dlasd4_()
726 eta = -w / dw; in dlasd4_()
730 temp = eta - temp1; in dlasd4_()
733 temp = eta - temp1; in dlasd4_()
735 if (temp > sg2ub || temp < sg2lb) { in dlasd4_()
737 eta = (sg2ub - tau) / 2.; in dlasd4_()
739 eta = (sg2lb - tau) / 2.; in dlasd4_()
743 tau += eta; in dlasd4_()
744 eta /= *sigma + sqrt(*sigma * *sigma + eta); in dlasd4_()
748 *sigma += eta; in dlasd4_()
751 work[j] += eta; in dlasd4_()
752 delta[j] -= eta; in dlasd4_()
763 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
764 psi += z__[j] * temp; in dlasd4_()
765 dpsi += temp * temp; in dlasd4_()
777 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
778 phi += z__[j] * temp; in dlasd4_()
779 dphi += temp * temp; in dlasd4_()
784 temp = z__[ii] / (work[ii] * delta[ii]); in dlasd4_()
785 dw = dpsi + dphi + temp * temp; in dlasd4_()
786 temp = z__[ii] * temp; in dlasd4_()
787 w = rhoinv + phi + psi + temp; in dlasd4_()
788 erretm = (phi - psi) * 8. + erretm + rhoinv * 2. + abs(temp) * 3. + in dlasd4_()
836 temp = z__[ii] / (work[ii] * delta[ii]); in dlasd4_()
838 dpsi += temp * temp; in dlasd4_()
840 dphi += temp * temp; in dlasd4_()
860 eta = b / a; in dlasd4_()
862 eta = (a - sqrt((d__1 = a * a - b * 4. * c__, abs(d__1)))) in dlasd4_()
865 eta = b * 2. / (a + sqrt((d__1 = a * a - b * 4. * c__, in dlasd4_()
874 temp = rhoinv + psi + phi; in dlasd4_()
876 c__ = temp - dtiim * dpsi - dtiip * dphi; in dlasd4_()
885 c__ = temp - dtiip * (dpsi + dphi) - temp2; in dlasd4_()
897 c__ = temp - dtiim * (dpsi + dphi) - temp2; in dlasd4_()
909 dlaed6_(&niter, &orgati, &c__, dd, zz, &w, &eta, info); in dlasd4_()
921 if (w * eta >= 0.) { in dlasd4_()
922 eta = -w / dw; in dlasd4_()
926 temp = eta - temp1; in dlasd4_()
929 temp = eta - temp1; in dlasd4_()
931 if (temp > sg2ub || temp < sg2lb) { in dlasd4_()
933 eta = (sg2ub - tau) / 2.; in dlasd4_()
935 eta = (sg2lb - tau) / 2.; in dlasd4_()
939 tau += eta; in dlasd4_()
940 eta /= *sigma + sqrt(*sigma * *sigma + eta); in dlasd4_()
942 *sigma += eta; in dlasd4_()
945 work[j] += eta; in dlasd4_()
946 delta[j] -= eta; in dlasd4_()
959 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
960 psi += z__[j] * temp; in dlasd4_()
961 dpsi += temp * temp; in dlasd4_()
973 temp = z__[j] / (work[j] * delta[j]); in dlasd4_()
974 phi += z__[j] * temp; in dlasd4_()
975 dphi += temp * temp; in dlasd4_()
980 temp = z__[ii] / (work[ii] * delta[ii]); in dlasd4_()
981 dw = dpsi + dphi + temp * temp; in dlasd4_()
982 temp = z__[ii] * temp; in dlasd4_()
983 w = rhoinv + phi + psi + temp; in dlasd4_()
984 erretm = (phi - psi) * 8. + erretm + rhoinv * 2. + abs(temp) * 3. in dlasd4_()