00001 /* dgennd.f -- translated by f2c (version 20061008). 00002 You must link the resulting object file with libf2c: 00003 on Microsoft Windows system, link with libf2c.lib; 00004 on Linux or Unix systems, link with .../path/to/libf2c.a -lm 00005 or, if you install libf2c.a in a standard place, with -lf2c -lm 00006 -- in that order, at the end of the command line, as in 00007 cc *.o -lf2c -lm 00008 Source for libf2c is in /netlib/f2c/libf2c.zip, e.g., 00009 00010 http://www.netlib.org/f2c/libf2c.zip 00011 */ 00012 00013 #include "f2c.h" 00014 #include "blaswrap.h" 00015 00016 logical dgennd_(integer *m, integer *n, doublereal *a, integer *lda) 00017 { 00018 /* System generated locals */ 00019 integer a_dim1, a_offset, i__1; 00020 logical ret_val; 00021 00022 /* Local variables */ 00023 integer i__, k; 00024 00025 00026 /* -- LAPACK test routine (version 3.1) -- */ 00027 /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */ 00028 /* February 2008 */ 00029 00030 /* .. Scalar Arguments .. */ 00031 /* .. */ 00032 /* .. Array Arguments .. */ 00033 /* .. */ 00034 00035 /* Purpose */ 00036 /* ======= */ 00037 00038 /* DGENND tests that its argument has a non-negative diagonal. */ 00039 00040 /* Arguments */ 00041 /* ========= */ 00042 00043 /* M (input) INTEGER */ 00044 /* The number of rows in A. */ 00045 00046 /* N (input) INTEGER */ 00047 /* The number of columns in A. */ 00048 00049 /* A (input) DOUBLE PRECISION array, dimension (LDA, N) */ 00050 /* The matrix. */ 00051 00052 /* LDA (input) INTEGER */ 00053 /* Leading dimension of A. */ 00054 00055 /* ===================================================================== */ 00056 00057 /* .. Parameters .. */ 00058 /* .. */ 00059 /* .. Local Scalars .. */ 00060 /* .. */ 00061 /* .. Intrinsics .. */ 00062 /* .. */ 00063 /* .. Executable Statements .. */ 00064 /* Parameter adjustments */ 00065 a_dim1 = *lda; 00066 a_offset = 1 + a_dim1; 00067 a -= a_offset; 00068 00069 /* Function Body */ 00070 k = min(*m,*n); 00071 i__1 = k; 00072 for (i__ = 1; i__ <= i__1; ++i__) { 00073 if (a[i__ + i__ * a_dim1] < 0.) { 00074 ret_val = FALSE_; 00075 return ret_val; 00076 } 00077 } 00078 ret_val = TRUE_; 00079 return ret_val; 00080 } /* dgennd_ */