35 #ifndef TEMPLATE_LAPACK_LACON_HEADER
36 #define TEMPLATE_LAPACK_LACON_HEADER
105 double d_sign(Treal *, Treal *);
112 Treal altsgn, estold;
122 for (i__ = 1; i__ <= i__1; ++i__) {
123 x[i__] = 1. / (Treal) (*n);
152 for (i__ = 1; i__ <= i__1; ++i__) {
153 x[i__] = d_sign(&c_b11, &x[i__]);
154 isgn[i__] = i_dnnt(&x[i__]);
172 for (i__ = 1; i__ <= i__1; ++i__) {
189 for (i__ = 1; i__ <= i__1; ++i__) {
190 d__1 = d_sign(&c_b11, &x[i__]);
191 if (i_dnnt(&d__1) != isgn[i__]) {
201 if (*est <= estold) {
206 for (i__ = 1; i__ <= i__1; ++i__) {
207 x[i__] = d_sign(&c_b11, &x[i__]);
208 isgn[i__] = i_dnnt(&x[i__]);
221 if (x[jlast] != (d__1 = x[j],
absMACRO(d__1)) && iter < 5) {
231 for (i__ = 1; i__ <= i__1; ++i__) {
232 x[i__] = altsgn * ((Treal) (i__ - 1) / (Treal) (*n - 1) +
#define absMACRO(x)
Definition: template_blas_common.h:45
integer template_blas_idamax(const integer *n, const Treal *dx, const integer *incx)
Definition: template_blas_idamax.h:40
int integer
Definition: template_blas_common.h:38
int template_lapack_lacon(const integer *n, Treal *v, Treal *x, integer *isgn, Treal *est, integer *kase)
Definition: template_lapack_lacon.h:40
Treal template_blas_asum(const integer *n, const Treal *dx, const integer *incx)
Definition: template_blas_asum.h:40
int template_blas_copy(const integer *n, const Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_copy.h:40