Subroutine */
| 30716 | } /* dlasq5_ */ |
| 30717 | |
| 30718 | /* Subroutine */ int dlasq6_(integer *i0, integer *n0, doublereal *z__, |
| 30719 | integer *pp, doublereal *dmin__, doublereal *dmin1, doublereal *dmin2, |
| 30720 | doublereal *dn, doublereal *dnm1, doublereal *dnm2) |
| 30721 | { |
| 30722 | /* System generated locals */ |
| 30723 | integer i__1; |
| 30724 | doublereal d__1, d__2; |
| 30725 | |
| 30726 | /* Local variables */ |
| 30727 | static doublereal d__; |
| 30728 | static integer j4, j4p2; |
| 30729 | static doublereal emin, temp; |
| 30730 | |
| 30731 | static doublereal safmin; |
| 30732 | |
| 30733 | |
| 30734 | /* |
| 30735 | -- LAPACK routine (version 3.2) -- |
| 30736 | |
| 30737 | -- Contributed by Osni Marques of the Lawrence Berkeley National -- |
| 30738 | -- Laboratory and Beresford Parlett of the Univ. of California at -- |
| 30739 | -- Berkeley -- |
| 30740 | -- November 2008 -- |
| 30741 | |
| 30742 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 30743 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 30744 | |
| 30745 | |
| 30746 | Purpose |
| 30747 | ======= |
| 30748 | |
| 30749 | DLASQ6 computes one dqd (shift equal to zero) transform in |
| 30750 | ping-pong form, with protection against underflow and overflow. |
| 30751 | |
| 30752 | Arguments |
| 30753 | ========= |
| 30754 | |
| 30755 | I0 (input) INTEGER |
| 30756 | First index. |
| 30757 | |
| 30758 | N0 (input) INTEGER |
| 30759 | Last index. |
| 30760 | |
| 30761 | Z (input) DOUBLE PRECISION array, dimension ( 4*N ) |
| 30762 | Z holds the qd array. EMIN is stored in Z(4*N0) to avoid |
| 30763 | an extra argument. |
| 30764 | |
| 30765 | PP (input) INTEGER |
| 30766 | PP=0 for ping, PP=1 for pong. |
| 30767 | |
| 30768 | DMIN (output) DOUBLE PRECISION |
| 30769 | Minimum value of d. |
| 30770 | |
| 30771 | DMIN1 (output) DOUBLE PRECISION |
| 30772 | Minimum value of d, excluding D( N0 ). |
| 30773 | |
| 30774 | DMIN2 (output) DOUBLE PRECISION |
| 30775 | Minimum value of d, excluding D( N0 ) and D( N0-1 ). |