37 #ifndef TEMPLATE_LAPACK_LACON_HEADER 38 #define TEMPLATE_LAPACK_LACON_HEADER 107 double d_sign(Treal *, Treal *);
114 Treal altsgn, estold;
124 for (i__ = 1; i__ <= i__1; ++i__) {
125 x[i__] = 1. / (Treal) (*n);
154 for (i__ = 1; i__ <= i__1; ++i__) {
155 x[i__] = d_sign(&c_b11, &x[i__]);
156 isgn[i__] = i_dnnt(&x[i__]);
174 for (i__ = 1; i__ <= i__1; ++i__) {
191 for (i__ = 1; i__ <= i__1; ++i__) {
192 d__1 = d_sign(&c_b11, &x[i__]);
193 if (i_dnnt(&d__1) != isgn[i__]) {
203 if (*est <= estold) {
208 for (i__ = 1; i__ <= i__1; ++i__) {
209 x[i__] = d_sign(&c_b11, &x[i__]);
210 isgn[i__] = i_dnnt(&x[i__]);
223 if (x[jlast] != (d__1 = x[j],
absMACRO(d__1)) && iter < 5) {
233 for (i__ = 1; i__ <= i__1; ++i__) {
234 x[i__] = altsgn * ((Treal) (i__ - 1) / (Treal) (*n - 1) +
#define absMACRO(x)
Definition: template_blas_common.h:47
integer template_blas_idamax(const integer *n, const Treal *dx, const integer *incx)
Definition: template_blas_idamax.h:42
int integer
Definition: template_blas_common.h:40
int template_lapack_lacon(const integer *n, Treal *v, Treal *x, integer *isgn, Treal *est, integer *kase)
Definition: template_lapack_lacon.h:42
Treal template_blas_asum(const integer *n, const Treal *dx, const integer *incx)
Definition: template_blas_asum.h:42
int template_blas_copy(const integer *n, const Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_copy.h:42