I am planning to use lapack/blas to attempt a least-squares fit of a
time series consisting of two million numerically calculated function
values, where each calculation has a relative numerical precision of
order 1.d-10. Many coefficients will likely be required to obtain a
good fit of this time series. I plan to use singular value decomposition
and statistical tests to identify which linear combinations of the
coefficients can be safely dropped from the fit.
Although I haven't yet programmed this fit, from my experience with
many other least-squares fits I know that ill-conditioning often
limits the number of terms that can be fit and therefore the quality
of the fit. To reduce the effect of that issue as much as possible I
am attempting to transform lapack/blas into quadruple- precision
libraries using the recently implemented quadruple-precision support
in gfortran-4.6.1.
One of the nice features of that support is the -fdefault-real-8
gfortran option which interprets the indefinite types of "real",
"complex", "double precision", and "double complex" as respectively
64-bit real, 128-bit complex, 128-bit real, and 256-bit complex. And
similarly for the constants 1.,(3.0,-1.5), 1.d0, and (3.0d0,-1.5d0).
The -fdefault-real-8 option has no effect on definite precision types
such as real*4, real*8, or complex*16 so you have to replace
those types by the corresponding real, double precision, and double
complex types to get the benefit of the -fdefault-real-8 option.
I used find and grep to scan through the lapack-3.1.1 source code and
found no instances of definite types other than COMPLEX*16. I
converted all of those to DOUBLE COMPLEX using the following commands:
cp -a lapack-3.3.1 lapack-3.3.1_double_complex
for FILE in $(find lapack-3.3.1_double_complex -name "*\.f*"); \
do echo $FILE; NAME=$(echo $FILE|sed "s?^.*/??"); echo $NAME; \
sed 's?COMPLEX\*16?DOUBLE COMPLEX?' <$FILE >| /tmp/$NAME; \
mv -f /tmp/$NAME $FILE; done
I then used the following method of building and testing
that changed code for the default double precision case
(as a comparison for the corresponding -fdefault-real-8 case
further below).
export FC="gfortran-4.6"
export FFLAGS="-O3 -fPIC -ffixed-line-length-132"
# The above -fPIC option is so I can link to the lapack/blas static
# libraries from a shared library, and the -ffixed-line-length-132
# option is to work around the issue that some lines have length
# greater than 72 bytes because of the replacement of "COMPLEX*16" by
# the longer "DOUBLE COMPLEX". N.B. The above FFLAGS does not contain
# the -fdefault-real-8 option so the result will have real, complex,
# double precision and double complex types that correspond to the
# usual 32-bit float, 64-bit complex, 64-bit float, and 128-bit
# complex types.
mkdir -p build_double_dir
cd build_double_dir
cmake ../lapack-3.3.1_double_complex >& cmake.out
make VERBOSE=1 -j4 >& make.out
ctest >& ctest_double.out
As expected, the cmake and make steps had no warnings or errors and ctest
sailed through the 98 tests of the resulting default precision
lapack/blas in 100 seconds with all tests
passed (see attached ctest_double.out file).
I then repeated the above steps with everything the same
except for using a different initially empty build tree, a
ctest_quadruple.out filename for the ctest output, and
using the -fdefault-real-8 gfortran option, i.e.,
export FFLAGS="-O3 -fPIC -fdefault-real-8 -ffixed-line-length-132"
# N.B. The above FFLAGS contains the -fdefault-real-8 option so the
# result will have real, complex, double precision and double complex
# types which correspond to the 64-bit float, 128-bit complex, 128-bit
# float, and 256-bit complex types.
As in the default precision case, the cmake and make steps completed
without any warning or error messages, but the ctest stage took 6000
seconds with 9 tests taking more than 100 seconds each and two of
those 9 timing out at 1500 seconds (!) On the other hand, many of the
test results completed in times comparable to the corresponding
default test times. (See the attached ctest_quadruple.out results).
So I think this approach for implementing a quadruple-precision
lapack/blas might be close to success.
Assuming that no errors are being introduced by gfortran-4.6 (as
supported by the excellent ctest results for the default precision
case), then I suspect the issue causing the above ctest test problems
for the -fdefault-real-8 precision case is caused by some assumptions
buried in the lapack/blas library code OR the lapack/blas test code
that are just not compatible with 128-bit floating point or 256-bit
complex calculations.
For example, there is an intriguing reference in
http://www.cs.berkeley.edu/~demmel/MS_Intel_ParLab_Apr08/Hida_high-precision-lapack.pdf
to providing a special safe minimum value returned for (d)lamch for
higher numerical precision results. So I have written a test routine
to see what is returned by dlamch for the 10 possible character
arguments. I have attached that test routine, and results returned
when that routine is compiled (with -fdefault-real-8 option for the
quadruple-precision case) and linked respectively against the default
and quadruple-precision lapack/blas libraries created above.
Do you see any results in dlamch_quadruple.out that are problematic
for quadruple-precision arithmetic?
The current approach is very similar to the "high-precision lapack"
approach other than my much simpler implementation because I am using
the gfortran -fdefault-real-8 option. So I think I might be close to
success for quadruple-precision lapack/blas, and I would appreciate
any advice that would allow me to achieve complete success (i.e., all
ctest test results passing within a time that may be larger than the
corresponding test time for the default precision tests, but not by
the huge factors shown for a small minority of the attached
ctest_quadruple.out test results when compared to the corresponding
attached ctest_double.out results).
Alan
__________________________
Alan W. Irwin
Astronomical research affiliation with Department of Physics and Astronomy,
University of Victoria (astrowww.phys.uvic.ca).
Programming affiliations with the FreeEOS equation-of-state
implementation for stellar interiors (freeeos.sf.net); the Time
Ephemerides project (timeephem.sf.net); PLplot scientific plotting
software package (plplot.sf.net); the libLASi project
(unifont.org/lasi); the Loads of Linux Links project (loll.sf.net);
and the Linux Brochure Project (lbproject.sf.net).
__________________________
Linux-powered Science
__________________________
-------------- next part --------------
Test project /home/software/lapack/build_double_dir
Start 1: xblat1s
1/98 Test #1: xblat1s .......................... Passed 0.02 sec
Start 2: xblat1c
2/98 Test #2: xblat1c .......................... Passed 0.01 sec
Start 3: xblat1d
3/98 Test #3: xblat1d .......................... Passed 0.01 sec
Start 4: xblat1z
4/98 Test #4: xblat1z .......................... Passed 0.02 sec
Start 5: xblat2s
5/98 Test #5: xblat2s .......................... Passed 0.07 sec
Start 6: xblat2c
6/98 Test #6: xblat2c .......................... Passed 0.04 sec
Start 7: xblat2d
7/98 Test #7: xblat2d .......................... Passed 0.03 sec
Start 8: xblat2z
8/98 Test #8: xblat2z .......................... Passed 0.05 sec
Start 9: xblat3s
9/98 Test #9: xblat3s .......................... Passed 0.05 sec
Start 10: xblat3c
10/98 Test #10: xblat3c .......................... Passed 0.07 sec
Start 11: xblat3d
11/98 Test #11: xblat3d .......................... Passed 0.05 sec
Start 12: xblat3z
12/98 Test #12: xblat3z .......................... Passed 0.06 sec
Start 13: xlintsts_stest_in
13/98 Test #13: xlintsts_stest_in ................ Passed 4.59 sec
Start 14: xlintstc_ctest_in
14/98 Test #14: xlintstc_ctest_in ................ Passed 13.10 sec
Start 15: xlintstd_dtest_in
15/98 Test #15: xlintstd_dtest_in ................ Passed 5.26 sec
Start 16: xlintstz_ztest_in
16/98 Test #16: xlintstz_ztest_in ................ Passed 16.50 sec
Start 17: xlintstds_dstest_in
17/98 Test #17: xlintstds_dstest_in .............. Passed 1.22 sec
Start 18: xlintstzc_zctest_in
18/98 Test #18: xlintstzc_zctest_in .............. Passed 3.12 sec
Start 19: xlintstrfs_stest_rfp_in
19/98 Test #19: xlintstrfs_stest_rfp_in .......... Passed 0.43 sec
Start 20: xlintstrfd_dtest_rfp_in
20/98 Test #20: xlintstrfd_dtest_rfp_in .......... Passed 0.45 sec
Start 21: xlintstrfc_ctest_rfp_in
21/98 Test #21: xlintstrfc_ctest_rfp_in .......... Passed 1.36 sec
Start 22: xlintstrfz_ztest_rfp_in
22/98 Test #22: xlintstrfz_ztest_rfp_in .......... Passed 1.62 sec
Start 23: xeigtsts_nep_in
23/98 Test #23: xeigtsts_nep_in .................. Passed 0.15 sec
Start 24: xeigtsts_sep_in
24/98 Test #24: xeigtsts_sep_in .................. Passed 0.79 sec
Start 25: xeigtsts_svd_in
25/98 Test #25: xeigtsts_svd_in .................. Passed 5.11 sec
Start 26: xeigtsts_sec_in
26/98 Test #26: xeigtsts_sec_in .................. Passed 0.58 sec
Start 27: xeigtsts_sed_in
27/98 Test #27: xeigtsts_sed_in .................. Passed 0.24 sec
Start 28: xeigtsts_sgg_in
28/98 Test #28: xeigtsts_sgg_in .................. Passed 0.15 sec
Start 29: xeigtsts_sgd_in
29/98 Test #29: xeigtsts_sgd_in .................. Passed 0.26 sec
Start 30: xeigtsts_ssb_in
30/98 Test #30: xeigtsts_ssb_in .................. Passed 0.02 sec
Start 31: xeigtsts_ssg_in
31/98 Test #31: xeigtsts_ssg_in .................. Passed 0.79 sec
Start 32: xeigtsts_sbal_in
32/98 Test #32: xeigtsts_sbal_in ................. Passed 0.02 sec
Start 33: xeigtsts_sbak_in
33/98 Test #33: xeigtsts_sbak_in ................. Passed 0.01 sec
Start 34: xeigtsts_sgbal_in
34/98 Test #34: xeigtsts_sgbal_in ................ Passed 0.02 sec
Start 35: xeigtsts_sgbak_in
35/98 Test #35: xeigtsts_sgbak_in ................ Passed 0.02 sec
Start 36: xeigtsts_sbb_in
36/98 Test #36: xeigtsts_sbb_in .................. Passed 0.04 sec
Start 37: xeigtsts_glm_in
37/98 Test #37: xeigtsts_glm_in .................. Passed 0.03 sec
Start 38: xeigtsts_gqr_in
38/98 Test #38: xeigtsts_gqr_in .................. Passed 0.04 sec
Start 39: xeigtsts_gsv_in
39/98 Test #39: xeigtsts_gsv_in .................. Passed 0.02 sec
Start 40: xeigtsts_csd_in
40/98 Test #40: xeigtsts_csd_in .................. Passed 0.01 sec
Start 41: xeigtsts_lse_in
41/98 Test #41: xeigtsts_lse_in .................. Passed 0.01 sec
Start 42: xeigtstc_nep_in
42/98 Test #42: xeigtstc_nep_in .................. Passed 0.25 sec
Start 43: xeigtstc_sep_in
43/98 Test #43: xeigtstc_sep_in .................. Passed 1.17 sec
Start 44: xeigtstc_svd_in
44/98 Test #44: xeigtstc_svd_in .................. Passed 8.44 sec
Start 45: xeigtstc_cec_in
45/98 Test #45: xeigtstc_cec_in .................. Passed 0.05 sec
Start 46: xeigtstc_ced_in
46/98 Test #46: xeigtstc_ced_in .................. Passed 0.35 sec
Start 47: xeigtstc_cgg_in
47/98 Test #47: xeigtstc_cgg_in .................. Passed 0.27 sec
Start 48: xeigtstc_cgd_in
48/98 Test #48: xeigtstc_cgd_in .................. Passed 0.29 sec
Start 49: xeigtstc_csb_in
49/98 Test #49: xeigtstc_csb_in .................. Passed 0.03 sec
Start 50: xeigtstc_csg_in
50/98 Test #50: xeigtstc_csg_in .................. Passed 1.14 sec
Start 51: xeigtstc_cbal_in
51/98 Test #51: xeigtstc_cbal_in ................. Passed 0.01 sec
Start 52: xeigtstc_cbak_in
52/98 Test #52: xeigtstc_cbak_in ................. Passed 0.01 sec
Start 53: xeigtstc_cgbal_in
53/98 Test #53: xeigtstc_cgbal_in ................ Passed 0.01 sec
Start 54: xeigtstc_cgbak_in
54/98 Test #54: xeigtstc_cgbak_in ................ Passed 0.01 sec
Start 55: xeigtstc_cbb_in
55/98 Test #55: xeigtstc_cbb_in .................. Passed 0.06 sec
Start 56: xeigtstc_glm_in
56/98 Test #56: xeigtstc_glm_in .................. Passed 0.03 sec
Start 57: xeigtstc_gqr_in
57/98 Test #57: xeigtstc_gqr_in .................. Passed 0.07 sec
Start 58: xeigtstc_gsv_in
58/98 Test #58: xeigtstc_gsv_in .................. Passed 0.04 sec
Start 59: xeigtstc_csd_in
59/98 Test #59: xeigtstc_csd_in .................. Passed 0.02 sec
Start 60: xeigtstc_lse_in
60/98 Test #60: xeigtstc_lse_in .................. Passed 0.02 sec
Start 61: xeigtstd_nep_in
61/98 Test #61: xeigtstd_nep_in .................. Passed 0.17 sec
Start 62: xeigtstd_sep_in
62/98 Test #62: xeigtstd_sep_in .................. Passed 1.16 sec
Start 63: xeigtstd_svd_in
63/98 Test #63: xeigtstd_svd_in .................. Passed 8.03 sec
Start 64: xeigtstd_dec_in
64/98 Test #64: xeigtstd_dec_in .................. Passed 0.58 sec
Start 65: xeigtstd_ded_in
65/98 Test #65: xeigtstd_ded_in .................. Passed 0.29 sec
Start 66: xeigtstd_dgg_in
66/98 Test #66: xeigtstd_dgg_in .................. Passed 0.19 sec
Start 67: xeigtstd_dgd_in
67/98 Test #67: xeigtstd_dgd_in .................. Passed 0.31 sec
Start 68: xeigtstd_dsb_in
68/98 Test #68: xeigtstd_dsb_in .................. Passed 0.02 sec
Start 69: xeigtstd_dsg_in
69/98 Test #69: xeigtstd_dsg_in .................. Passed 1.25 sec
Start 70: xeigtstd_dbal_in
70/98 Test #70: xeigtstd_dbal_in ................. Passed 0.02 sec
Start 71: xeigtstd_dbak_in
71/98 Test #71: xeigtstd_dbak_in ................. Passed 0.01 sec
Start 72: xeigtstd_dgbal_in
72/98 Test #72: xeigtstd_dgbal_in ................ Passed 0.01 sec
Start 73: xeigtstd_dgbak_in
73/98 Test #73: xeigtstd_dgbak_in ................ Passed 0.01 sec
Start 74: xeigtstd_dbb_in
74/98 Test #74: xeigtstd_dbb_in .................. Passed 0.05 sec
Start 75: xeigtstd_glm_in
75/98 Test #75: xeigtstd_glm_in .................. Passed 0.02 sec
Start 76: xeigtstd_gqr_in
76/98 Test #76: xeigtstd_gqr_in .................. Passed 0.04 sec
Start 77: xeigtstd_gsv_in
77/98 Test #77: xeigtstd_gsv_in .................. Passed 0.03 sec
Start 78: xeigtstd_csd_in
78/98 Test #78: xeigtstd_csd_in .................. Passed 0.01 sec
Start 79: xeigtstd_lse_in
79/98 Test #79: xeigtstd_lse_in .................. Passed 0.01 sec
Start 80: xeigtstz_nep_in
80/98 Test #80: xeigtstz_nep_in .................. Passed 0.34 sec
Start 81: xeigtstz_sep_in
81/98 Test #81: xeigtstz_sep_in .................. Passed 1.97 sec
Start 82: xeigtstz_svd_in
82/98 Test #82: xeigtstz_svd_in .................. Passed 12.74 sec
Start 83: xeigtstz_zec_in
83/98 Test #83: xeigtstz_zec_in .................. Passed 0.05 sec
Start 84: xeigtstz_zed_in
84/98 Test #84: xeigtstz_zed_in .................. Passed 0.48 sec
Start 85: xeigtstz_zgg_in
85/98 Test #85: xeigtstz_zgg_in .................. Passed 0.40 sec
Start 86: xeigtstz_zgd_in
86/98 Test #86: xeigtstz_zgd_in .................. Passed 0.38 sec
Start 87: xeigtstz_zsb_in
87/98 Test #87: xeigtstz_zsb_in .................. Passed 0.03 sec
Start 88: xeigtstz_zsg_in
88/98 Test #88: xeigtstz_zsg_in .................. Passed 1.72 sec
Start 89: xeigtstz_zbal_in
89/98 Test #89: xeigtstz_zbal_in ................. Passed 0.02 sec
Start 90: xeigtstz_zbak_in
90/98 Test #90: xeigtstz_zbak_in ................. Passed 0.01 sec
Start 91: xeigtstz_zgbal_in
91/98 Test #91: xeigtstz_zgbal_in ................ Passed 0.02 sec
Start 92: xeigtstz_zgbak_in
92/98 Test #92: xeigtstz_zgbak_in ................ Passed 0.02 sec
Start 93: xeigtstz_zbb_in
93/98 Test #93: xeigtstz_zbb_in .................. Passed 0.07 sec
Start 94: xeigtstz_glm_in
94/98 Test #94: xeigtstz_glm_in .................. Passed 0.03 sec
Start 95: xeigtstz_gqr_in
95/98 Test #95: xeigtstz_gqr_in .................. Passed 0.08 sec
Start 96: xeigtstz_gsv_in
96/98 Test #96: xeigtstz_gsv_in .................. Passed 0.05 sec
Start 97: xeigtstz_csd_in
97/98 Test #97: xeigtstz_csd_in .................. Passed 0.02 sec
Start 98: xeigtstz_lse_in
98/98 Test #98: xeigtstz_lse_in .................. Passed 0.02 sec
100% tests passed, 0 tests failed out of 98
Total Test time (real) = 99.64 sec
-------------- next part --------------
Test project /home/software/lapack/build_quadruple_dir
Start 1: xblat1s
1/98 Test #1: xblat1s .......................... Passed 0.05 sec
Start 2: xblat1c
2/98 Test #2: xblat1c .......................... Passed 0.01 sec
Start 3: xblat1d
3/98 Test #3: xblat1d .......................... Passed 0.03 sec
Start 4: xblat1z
4/98 Test #4: xblat1z .......................... Passed 0.01 sec
Start 5: xblat2s
5/98 Test #5: xblat2s .......................... Passed 0.03 sec
Start 6: xblat2c
6/98 Test #6: xblat2c .......................... Passed 0.04 sec
Start 7: xblat2d
7/98 Test #7: xblat2d .......................... Passed 0.50 sec
Start 8: xblat2z
8/98 Test #8: xblat2z .......................... Passed 1.86 sec
Start 9: xblat3s
9/98 Test #9: xblat3s .......................... Passed 0.04 sec
Start 10: xblat3c
10/98 Test #10: xblat3c .......................... Passed 0.06 sec
Start 11: xblat3d
11/98 Test #11: xblat3d .......................... Passed 1.24 sec
Start 12: xblat3z
12/98 Test #12: xblat3z .......................... Passed 3.96 sec
Start 13: xlintsts_stest_in
13/98 Test #13: xlintsts_stest_in ................ Passed 5.26 sec
Start 14: xlintstc_ctest_in
14/98 Test #14: xlintstc_ctest_in ................ Passed 16.53 sec
Start 15: xlintstd_dtest_in
15/98 Test #15: xlintstd_dtest_in ................ Passed 495.30 sec
Start 16: xlintstz_ztest_in
16/98 Test #16: xlintstz_ztest_in ................***Timeout 1500.01 sec
Start 17: xlintstds_dstest_in
17/98 Test #17: xlintstds_dstest_in .............. Passed 174.01 sec
Start 18: xlintstzc_zctest_in
18/98 Test #18: xlintstzc_zctest_in .............. Passed 551.35 sec
Start 19: xlintstrfs_stest_rfp_in
19/98 Test #19: xlintstrfs_stest_rfp_in .......... Passed 0.48 sec
Start 20: xlintstrfd_dtest_rfp_in
20/98 Test #20: xlintstrfd_dtest_rfp_in .......... Passed 46.00 sec
Start 21: xlintstrfc_ctest_rfp_in
21/98 Test #21: xlintstrfc_ctest_rfp_in .......... Passed 1.64 sec
Start 22: xlintstrfz_ztest_rfp_in
22/98 Test #22: xlintstrfz_ztest_rfp_in .......... Passed 257.58 sec
Start 23: xeigtsts_nep_in
23/98 Test #23: xeigtsts_nep_in .................. Passed 0.21 sec
Start 24: xeigtsts_sep_in
24/98 Test #24: xeigtsts_sep_in .................. Passed 1.43 sec
Start 25: xeigtsts_svd_in
25/98 Test #25: xeigtsts_svd_in .................. Passed 8.01 sec
Start 26: xeigtsts_sec_in
26/98 Test #26: xeigtsts_sec_in .................. Passed 0.57 sec
Start 27: xeigtsts_sed_in
27/98 Test #27: xeigtsts_sed_in .................. Passed 0.29 sec
Start 28: xeigtsts_sgg_in
28/98 Test #28: xeigtsts_sgg_in .................. Passed 0.19 sec
Start 29: xeigtsts_sgd_in
29/98 Test #29: xeigtsts_sgd_in .................. Passed 0.31 sec
Start 30: xeigtsts_ssb_in
30/98 Test #30: xeigtsts_ssb_in .................. Passed 0.02 sec
Start 31: xeigtsts_ssg_in
31/98 Test #31: xeigtsts_ssg_in .................. Passed 1.24 sec
Start 32: xeigtsts_sbal_in
32/98 Test #32: xeigtsts_sbal_in ................. Passed 0.01 sec
Start 33: xeigtsts_sbak_in
33/98 Test #33: xeigtsts_sbak_in ................. Passed 0.01 sec
Start 34: xeigtsts_sgbal_in
34/98 Test #34: xeigtsts_sgbal_in ................ Passed 0.01 sec
Start 35: xeigtsts_sgbak_in
35/98 Test #35: xeigtsts_sgbak_in ................ Passed 0.01 sec
Start 36: xeigtsts_sbb_in
36/98 Test #36: xeigtsts_sbb_in .................. Passed 0.04 sec
Start 37: xeigtsts_glm_in
37/98 Test #37: xeigtsts_glm_in .................. Passed 0.02 sec
Start 38: xeigtsts_gqr_in
38/98 Test #38: xeigtsts_gqr_in .................. Passed 0.04 sec
Start 39: xeigtsts_gsv_in
39/98 Test #39: xeigtsts_gsv_in .................. Passed 0.03 sec
Start 40: xeigtsts_csd_in
40/98 Test #40: xeigtsts_csd_in .................. Passed 0.04 sec
Start 41: xeigtsts_lse_in
41/98 Test #41: xeigtsts_lse_in .................. Passed 0.01 sec
Start 42: xeigtstc_nep_in
42/98 Test #42: xeigtstc_nep_in .................. Passed 0.35 sec
Start 43: xeigtstc_sep_in
43/98 Test #43: xeigtstc_sep_in .................. Passed 1.98 sec
Start 44: xeigtstc_svd_in
44/98 Test #44: xeigtstc_svd_in .................. Passed 12.80 sec
Start 45: xeigtstc_cec_in
45/98 Test #45: xeigtstc_cec_in .................. Passed 0.05 sec
Start 46: xeigtstc_ced_in
46/98 Test #46: xeigtstc_ced_in .................. Passed 0.48 sec
Start 47: xeigtstc_cgg_in
47/98 Test #47: xeigtstc_cgg_in .................. Passed 0.39 sec
Start 48: xeigtstc_cgd_in
48/98 Test #48: xeigtstc_cgd_in .................. Passed 0.38 sec
Start 49: xeigtstc_csb_in
49/98 Test #49: xeigtstc_csb_in .................. Passed 0.03 sec
Start 50: xeigtstc_csg_in
50/98 Test #50: xeigtstc_csg_in .................. Passed 1.70 sec
Start 51: xeigtstc_cbal_in
51/98 Test #51: xeigtstc_cbal_in ................. Passed 0.01 sec
Start 52: xeigtstc_cbak_in
52/98 Test #52: xeigtstc_cbak_in ................. Passed 0.01 sec
Start 53: xeigtstc_cgbal_in
53/98 Test #53: xeigtstc_cgbal_in ................ Passed 0.02 sec
Start 54: xeigtstc_cgbak_in
54/98 Test #54: xeigtstc_cgbak_in ................ Passed 0.01 sec
Start 55: xeigtstc_cbb_in
55/98 Test #55: xeigtstc_cbb_in .................. Passed 0.07 sec
Start 56: xeigtstc_glm_in
56/98 Test #56: xeigtstc_glm_in .................. Passed 0.03 sec
Start 57: xeigtstc_gqr_in
57/98 Test #57: xeigtstc_gqr_in .................. Passed 0.08 sec
Start 58: xeigtstc_gsv_in
58/98 Test #58: xeigtstc_gsv_in .................. Passed 0.05 sec
Start 59: xeigtstc_csd_in
59/98 Test #59: xeigtstc_csd_in .................. Passed 0.02 sec
Start 60: xeigtstc_lse_in
60/98 Test #60: xeigtstc_lse_in .................. Passed 0.02 sec
Start 61: xeigtstd_nep_in
61/98 Test #61: xeigtstd_nep_in .................. Passed 10.37 sec
Start 62: xeigtstd_sep_in
62/98 Test #62: xeigtstd_sep_in .................. Passed 69.46 sec
Start 63: xeigtstd_svd_in
63/98 Test #63: xeigtstd_svd_in .................. Passed 820.82 sec
Start 64: xeigtstd_dec_in
64/98 Test #64: xeigtstd_dec_in .................. Passed 14.64 sec
Start 65: xeigtstd_ded_in
65/98 Test #65: xeigtstd_ded_in .................. Passed 13.62 sec
Start 66: xeigtstd_dgg_in
66/98 Test #66: xeigtstd_dgg_in .................. Passed 15.16 sec
Start 67: xeigtstd_dgd_in
67/98 Test #67: xeigtstd_dgd_in .................. Passed 14.14 sec
Start 68: xeigtstd_dsb_in
68/98 Test #68: xeigtstd_dsb_in .................. Passed 0.56 sec
Start 69: xeigtstd_dsg_in
69/98 Test #69: xeigtstd_dsg_in .................. Passed 66.08 sec
Start 70: xeigtstd_dbal_in
70/98 Test #70: xeigtstd_dbal_in ................. Passed 0.01 sec
Start 71: xeigtstd_dbak_in
71/98 Test #71: xeigtstd_dbak_in ................. Passed 0.01 sec
Start 72: xeigtstd_dgbal_in
72/98 Test #72: xeigtstd_dgbal_in ................ Passed 0.01 sec
Start 73: xeigtstd_dgbak_in
73/98 Test #73: xeigtstd_dgbak_in ................ Passed 0.01 sec
Start 74: xeigtstd_dbb_in
74/98 Test #74: xeigtstd_dbb_in .................. Passed 1.12 sec
Start 75: xeigtstd_glm_in
75/98 Test #75: xeigtstd_glm_in .................. Passed 1.34 sec
Start 76: xeigtstd_gqr_in
76/98 Test #76: xeigtstd_gqr_in .................. Passed 2.71 sec
Start 77: xeigtstd_gsv_in
77/98 Test #77: xeigtstd_gsv_in .................. Passed 2.28 sec
Start 78: xeigtstd_csd_in
78/98 Test #78: xeigtstd_csd_in .................. Passed 0.66 sec
Start 79: xeigtstd_lse_in
79/98 Test #79: xeigtstd_lse_in .................. Passed 0.70 sec
Start 80: xeigtstz_nep_in
80/98 Test #80: xeigtstz_nep_in .................. Passed 37.76 sec
Start 81: xeigtstz_sep_in
81/98 Test #81: xeigtstz_sep_in .................. Passed 163.03 sec
Start 82: xeigtstz_svd_in
82/98 Test #82: xeigtstz_svd_in ..................***Timeout 1500.01 sec
Start 83: xeigtstz_zec_in
83/98 Test #83: xeigtstz_zec_in .................. Passed 1.40 sec
Start 84: xeigtstz_zed_in
84/98 Test #84: xeigtstz_zed_in .................. Passed 47.78 sec
Start 85: xeigtstz_zgg_in
85/98 Test #85: xeigtstz_zgg_in .................. Passed 48.51 sec
Start 86: xeigtstz_zgd_in
86/98 Test #86: xeigtstz_zgd_in .................. Passed 35.53 sec
Start 87: xeigtstz_zsb_in
87/98 Test #87: xeigtstz_zsb_in .................. Passed 2.04 sec
Start 88: xeigtstz_zsg_in
88/98 Test #88: xeigtstz_zsg_in .................. Passed 134.13 sec
Start 89: xeigtstz_zbal_in
89/98 Test #89: xeigtstz_zbal_in ................. Passed 0.02 sec
Start 90: xeigtstz_zbak_in
90/98 Test #90: xeigtstz_zbak_in ................. Passed 0.01 sec
Start 91: xeigtstz_zgbal_in
91/98 Test #91: xeigtstz_zgbal_in ................ Passed 0.02 sec
Start 92: xeigtstz_zgbak_in
92/98 Test #92: xeigtstz_zgbak_in ................ Passed 0.01 sec
Start 93: xeigtstz_zbb_in
93/98 Test #93: xeigtstz_zbb_in .................. Passed 4.23 sec
Start 94: xeigtstz_glm_in
94/98 Test #94: xeigtstz_glm_in .................. Passed 4.15 sec
Start 95: xeigtstz_gqr_in
95/98 Test #95: xeigtstz_gqr_in .................. Passed 9.03 sec
Start 96: xeigtstz_gsv_in
96/98 Test #96: xeigtstz_gsv_in .................. Passed 7.56 sec
Start 97: xeigtstz_csd_in
97/98 Test #97: xeigtstz_csd_in .................. Passed 1.92 sec
Start 98: xeigtstz_lse_in
98/98 Test #98: xeigtstz_lse_in .................. Passed 2.13 sec
98% tests passed, 2 tests failed out of 98
Total Test time (real) = 6120.37 sec
The following tests FAILED:
16 - xlintstz_ztest_in (Timeout)
82 - xeigtstz_svd_in (Timeout)
Errors while running CTest
-------------- next part --------------
implicit none
double precision dlamch
character*1 cmach(10)
integer i
data cmach/"e", "s", "b", "p", "n", "r", "m", "u", "l", "o"/
do i = 1,10
write(*,*), "cmach(i), dlamch(cmach(i)) =",
& cmach(i), dlamch(cmach(i))
enddo
end
-------------- next part --------------
cmach(i), dlamch(cmach(i)) =e 1.11022302462515654E-016
cmach(i), dlamch(cmach(i)) =s 2.22507385850720138E-308
cmach(i), dlamch(cmach(i)) =b 2.0000000000000000
cmach(i), dlamch(cmach(i)) =p 2.22044604925031308E-016
cmach(i), dlamch(cmach(i)) =n 53.000000000000000
cmach(i), dlamch(cmach(i)) =r 1.0000000000000000
cmach(i), dlamch(cmach(i)) =m -1021.0000000000000
cmach(i), dlamch(cmach(i)) =u 2.22507385850720138E-308
cmach(i), dlamch(cmach(i)) =l 1024.0000000000000
cmach(i), dlamch(cmach(i)) =o 1.79769313486231571E+308
-------------- next part --------------
cmach(i), dlamch(cmach(i)) =e 9.62964972193617926527988971292463659E-0035
cmach(i), dlamch(cmach(i)) =s 3.36210314311209350626267781732175260E-4932
cmach(i), dlamch(cmach(i)) =b 2.0000000000000000000000000000000000
cmach(i), dlamch(cmach(i)) =p 1.92592994438723585305597794258492732E-0034
cmach(i), dlamch(cmach(i)) =n 113.00000000000000000000000000000000
cmach(i), dlamch(cmach(i)) =r 1.0000000000000000000000000000000000
cmach(i), dlamch(cmach(i)) =m -16381.000000000000000000000000000000
cmach(i), dlamch(cmach(i)) =u 3.36210314311209350626267781732175260E-4932
cmach(i), dlamch(cmach(i)) =l 16384.000000000000000000000000000000
cmach(i), dlamch(cmach(i)) =o 1.18973149535723176508575932662800702E+4932
|