lemon-project-template-glpk

view deps/glpk/src/amd/amd_info.c @ 11:4fc6ad2fb8a6

Test GLPK in src/main.cc
author Alpar Juttner <alpar@cs.elte.hu>
date Sun, 06 Nov 2011 21:43:29 +0100
parents
children
line source
1 /* ========================================================================= */
2 /* === AMD_info ============================================================ */
3 /* ========================================================================= */
5 /* ------------------------------------------------------------------------- */
6 /* AMD, Copyright (c) Timothy A. Davis, */
7 /* Patrick R. Amestoy, and Iain S. Duff. See ../README.txt for License. */
8 /* email: davis at cise.ufl.edu CISE Department, Univ. of Florida. */
9 /* web: http://www.cise.ufl.edu/research/sparse/amd */
10 /* ------------------------------------------------------------------------- */
12 /* User-callable. Prints the output statistics for AMD. See amd.h
13 * for details. If the Info array is not present, nothing is printed.
14 */
16 #include "amd_internal.h"
18 #define PRI(format,x) { if (x >= 0) { PRINTF ((format, x)) ; }}
20 GLOBAL void AMD_info
21 (
22 double Info [ ]
23 )
24 {
25 double n, ndiv, nmultsubs_ldl, nmultsubs_lu, lnz, lnzd ;
27 PRINTF (("\nAMD version %d.%d.%d, %s, results:\n",
28 AMD_MAIN_VERSION, AMD_SUB_VERSION, AMD_SUBSUB_VERSION, AMD_DATE)) ;
30 if (!Info)
31 {
32 return ;
33 }
35 n = Info [AMD_N] ;
36 ndiv = Info [AMD_NDIV] ;
37 nmultsubs_ldl = Info [AMD_NMULTSUBS_LDL] ;
38 nmultsubs_lu = Info [AMD_NMULTSUBS_LU] ;
39 lnz = Info [AMD_LNZ] ;
40 lnzd = (n >= 0 && lnz >= 0) ? (n + lnz) : (-1) ;
42 /* AMD return status */
43 PRINTF ((" status: ")) ;
44 if (Info [AMD_STATUS] == AMD_OK)
45 {
46 PRINTF (("OK\n")) ;
47 }
48 else if (Info [AMD_STATUS] == AMD_OUT_OF_MEMORY)
49 {
50 PRINTF (("out of memory\n")) ;
51 }
52 else if (Info [AMD_STATUS] == AMD_INVALID)
53 {
54 PRINTF (("invalid matrix\n")) ;
55 }
56 else if (Info [AMD_STATUS] == AMD_OK_BUT_JUMBLED)
57 {
58 PRINTF (("OK, but jumbled\n")) ;
59 }
60 else
61 {
62 PRINTF (("unknown\n")) ;
63 }
65 /* statistics about the input matrix */
66 PRI (" n, dimension of A: %.20g\n", n);
67 PRI (" nz, number of nonzeros in A: %.20g\n",
68 Info [AMD_NZ]) ;
69 PRI (" symmetry of A: %.4f\n",
70 Info [AMD_SYMMETRY]) ;
71 PRI (" number of nonzeros on diagonal: %.20g\n",
72 Info [AMD_NZDIAG]) ;
73 PRI (" nonzeros in pattern of A+A' (excl. diagonal): %.20g\n",
74 Info [AMD_NZ_A_PLUS_AT]) ;
75 PRI (" # dense rows/columns of A+A': %.20g\n",
76 Info [AMD_NDENSE]) ;
78 /* statistics about AMD's behavior */
79 PRI (" memory used, in bytes: %.20g\n",
80 Info [AMD_MEMORY]) ;
81 PRI (" # of memory compactions: %.20g\n",
82 Info [AMD_NCMPA]) ;
84 /* statistics about the ordering quality */
85 PRINTF (("\n"
86 " The following approximate statistics are for a subsequent\n"
87 " factorization of A(P,P) + A(P,P)'. They are slight upper\n"
88 " bounds if there are no dense rows/columns in A+A', and become\n"
89 " looser if dense rows/columns exist.\n\n")) ;
91 PRI (" nonzeros in L (excluding diagonal): %.20g\n",
92 lnz) ;
93 PRI (" nonzeros in L (including diagonal): %.20g\n",
94 lnzd) ;
95 PRI (" # divide operations for LDL' or LU: %.20g\n",
96 ndiv) ;
97 PRI (" # multiply-subtract operations for LDL': %.20g\n",
98 nmultsubs_ldl) ;
99 PRI (" # multiply-subtract operations for LU: %.20g\n",
100 nmultsubs_lu) ;
101 PRI (" max nz. in any column of L (incl. diagonal): %.20g\n",
102 Info [AMD_DMAX]) ;
104 /* total flop counts for various factorizations */
106 if (n >= 0 && ndiv >= 0 && nmultsubs_ldl >= 0 && nmultsubs_lu >= 0)
107 {
108 PRINTF (("\n"
109 " chol flop count for real A, sqrt counted as 1 flop: %.20g\n"
110 " LDL' flop count for real A: %.20g\n"
111 " LDL' flop count for complex A: %.20g\n"
112 " LU flop count for real A (with no pivoting): %.20g\n"
113 " LU flop count for complex A (with no pivoting): %.20g\n\n",
114 n + ndiv + 2*nmultsubs_ldl,
115 ndiv + 2*nmultsubs_ldl,
116 9*ndiv + 8*nmultsubs_ldl,
117 ndiv + 2*nmultsubs_lu,
118 9*ndiv + 8*nmultsubs_lu)) ;
119 }
120 }