Différences

Ci-dessous, les différences entre deux révisions de la page.

--- arpack [2009/01/12 08:55] – gerard
+++ arpack [2009/01/12 09:40] – gerard
@@ Ligne 42: / Ligne 42: @@
 ---
 > DIRS         = $(UTILdir) $(SRCdir)
+---
+> DIRS         = $(UTILdir) $(SRCdir)
+c61
+< LAPACKLIB =
+---
+> LAPACKLIB = -xlic_lib=sunperf
 c105
 < FFLAGS        = -O -cg89
 ---
-> FFLAGS        = -fast
+> FFLAGS        = -dalign
 c107
 < LDFLAGS =
@@ Ligne 62: / Ligne 68: @@
 </code>
-mais l'exemple ne passe pas en 32bits (FIXME: ne pas compiler avec l'option fast!), il faut compiler en 64 bits
 <code>
 > file dssimp
@@ Ligne 68: / Ligne 74: @@
 > ./dssimp
-  Error with _saupd, info =   -9999
+ _saupd: number of update iterations taken
-  Check documentation in _saupd
+ -----------------------------------------
+-    1:     8
- Note: IEEE floating-point exception flags raised:
-    Inexact;  Invalid Operation;
- See the Numerical Computation Guide, ieee_flags(3M)
+ _saupd: number of "converged" Ritz values
+ -----------------------------------------
+-    1:     4
+ _saupd: final Ritz values
+ -------------------------
+-    4:   8.912D+02   9.198D+02   9.198D+02   9.484D+02
+ _saupd: corresponding error bounds
+ ----------------------------------
+-    4:   6.399D-24   1.169D-14   7.116D-17   1.637D-35
+     ==========================================
+     = Symmetric implicit Arnoldi update code =
+     = Version Number: 2.4                    =
+     = Version Date:   07/31/96               =
+     ==========================================
+     = Summary of timing statistics           =
+     ==========================================
+     Total number update iterations             =     8
+     Total number of OP*x operations            =   125
+     Total number of B*x operations             =     0
+     Total number of reorthogonalization steps  =   125
+     Total number of iterative refinement steps =     0
+     Total number of restart steps              =     0
+     Total time in user OP*x operation          =     0.000467
+     Total time in user B*x operation           =     0.000000
+     Total time in Arnoldi update routine       =     0.005470
+     Total time in saup2 routine                =     0.004773
+     Total time in basic Arnoldi iteration loop =     0.001817
+     Total time in reorthogonalization phase    =     0.000646
+     Total time in (re)start vector generation  =     0.000054
+     Total time in trid eigenvalue subproblem   =     0.000684
+     Total time in getting the shifts           =     0.000054
+     Total time in applying the shifts          =     0.002106
+     Total time in convergence testing          =     0.000007
+ Ritz values and relative residuals
+ ----------------------------------
+               Col   1       Col   2
+  Row   1:    8.91167D+02   1.78240D-15
+  Row   2:    9.19781D+02   1.44880D-15
+  Row   3:    9.19781D+02   2.07175D-15
+  Row   4:    9.48395D+02   2.57489D-15
+  _SSIMP
+  ======
+  Size of the matrix is   100
+  The number of Ritz values requested is   4
+  The number of Arnoldi vectors generated (NCV) is   20
+  What portion of the spectrum: LM
+  The number of converged Ritz values is   4
+  The number of Implicit Arnoldi update iterations taken is   8
+  The number of OP*x is   125
+  The convergence criterion is     1.1102230246252D-16
+ Note: IEEE floating-point exception flags raised:
+    Inexact;  Underflow;
+ See the Numerical Computation Guide, ieee_flags(3M)
 </code>