06-Feb-2019 08:15:27 kelley_test MATLAB version. Test kelley. KELLEY_TEST01 Test GMRES on the simple -1,2-1 matrix of order N = 20 Before GMRES, ||X_ERROR|| = 5.357238e+01 Before GMRES, ||A*x-b|| = 2.100000e+01 After GMRES, ||X_ERROR|| = 5.240347e-14 After GMRES, ||A*x-b|| = 1.414652e-14 GMRES took 20 iterations GMRES residual error history: 1 2.100000e+01 2 9.391486e+00 3 5.612486e+00 4 3.834058e+00 5 2.831639e+00 6 2.201398e+00 7 1.774824e+00 8 1.470294e+00 9 1.243933e+00 10 1.070259e+00 11 9.335639e-01 12 8.236878e-01 13 7.337994e-01 14 6.591531e-01 15 5.963599e-01 16 5.429421e-01 17 4.970501e-01 18 4.572787e-01 19 4.225429e-01 20 3.919930e-01 21 0.000000e+00 Before GMRES, ||X_ERROR|| = 5.120581e+01 Before GMRES, ||A*x-b|| = 2.000876e+01 After GMRES, ||X_ERROR|| = 4.772978e-14 After GMRES, ||A*x-b|| = 2.135210e-14 GMRES took 20 iterations GMRES residual error history: 1 2.000876e+01 2 8.792554e+00 3 5.351533e+00 4 3.669813e+00 5 2.704665e+00 6 2.122435e+00 7 1.708333e+00 8 1.414728e+00 9 1.204133e+00 10 1.030822e+00 11 8.928321e-01 12 7.798649e-01 13 6.940372e-01 14 6.236210e-01 15 5.636029e-01 16 5.136067e-01 17 4.700804e-01 18 4.323902e-01 19 4.020984e-01 20 3.741155e-01 21 7.634256e-29 KELLEY_TEST02 Test TFQMR on the simple -1,2-1 matrix of order N = 20 Before TFQMR, ||X_ERROR|| = 5.357238e+01 Before TFQMR, ||A*x-b|| = 2.100000e+01 After TFQMR, ||X_ERROR|| = 9.691427e-15 After TFQMR, ||A*x-b|| = 1.527757e-14 TFQMR took 20 iterations TFQMR residual error history: 1 2.100000e+01 2 5.824352e+00 3 2.757435e+00 4 1.610626e+00 5 1.056625e+00 6 7.466740e-01 7 5.557189e-01 8 4.297364e-01 9 3.422448e-01 10 2.790097e-01 11 2.299639e-01 12 1.872382e-01 13 1.515891e-01 14 1.224511e-01 15 9.877708e-02 16 7.947142e-02 17 6.354826e-02 18 5.014145e-02 19 3.840151e-02 20 2.703330e-02 21 2.425658e-15 Before TFQMR, ||X_ERROR|| = 5.215078e+01 Before TFQMR, ||A*x-b|| = 2.207620e+01 After TFQMR, ||X_ERROR|| = 1.291872e-14 After TFQMR, ||A*x-b|| = 2.485412e-14 TFQMR took 20 iterations TFQMR residual error history: 1 2.207620e+01 2 5.510127e+00 3 2.581347e+00 4 1.571795e+00 5 1.020879e+00 6 7.086899e-01 7 5.127420e-01 8 3.811200e-01 9 2.989733e-01 10 2.376329e-01 11 1.884997e-01 12 1.489086e-01 13 1.213892e-01 14 9.953533e-02 15 8.188868e-02 16 6.735469e-02 17 5.414254e-02 18 4.436310e-02 19 3.566820e-02 20 2.549748e-02 21 6.969579e-15 KELLEY_TEST03 Test BICGSTAB on the simple -1,2-1 matrix of order N = 20 Before BICGSTAB, ||X_ERROR|| = 5.357238e+01 Before BICGSTAB, ||A*x-b|| = 2.100000e+01 After BICGSTAB, ||X_ERROR|| = 6.558253e-09 After BICGSTAB, ||A*x-b|| = 2.183363e-10 BICGSTAB took 20 iterations BICGSTAB residual error history: 1 2.100000e+01 2 6.062178e+00 3 3.299832e+00 4 2.213594e+00 5 1.635007e+00 6 1.281864e+00 7 1.045450e+00 8 8.774544e-01 9 7.524163e-01 10 6.561685e-01 11 5.799848e-01 12 5.183222e-01 13 4.669121e-01 14 4.212527e-01 15 3.753684e-01 16 3.239107e-01 17 2.643018e-01 18 1.977284e-01 19 1.282688e-01 20 6.086974e-02 21 2.183366e-10 Before BICGSTAB, ||X_ERROR|| = 5.150157e+01 Before BICGSTAB, ||A*x-b|| = 2.011416e+01 After BICGSTAB, ||X_ERROR|| = 3.319787e-10 After BICGSTAB, ||A*x-b|| = 8.410791e-12 BICGSTAB took 20 iterations BICGSTAB residual error history: 1 2.011416e+01 2 5.841420e+00 3 3.202741e+00 4 2.133794e+00 5 1.597559e+00 6 1.248405e+00 7 1.008487e+00 8 8.303602e-01 9 7.088704e-01 10 6.142410e-01 11 5.443437e-01 12 4.889824e-01 13 4.383012e-01 14 3.935513e-01 15 3.423691e-01 16 2.860595e-01 17 2.232568e-01 18 1.575175e-01 19 9.333852e-02 20 3.698558e-02 21 8.410801e-12 KELLEY_TEST04 Test GMRESB on the simple -1,2-1 matrix of order N = 20 Before GMRESB, ||X_ERROR|| = 5.357238e+01 Before GMRESB, ||A*x-b|| = 2.100000e+01 After GMRESB, ||X_ERROR|| = 1.671256e-13 After GMRESB, ||A*x-b|| = 4.925150e-14 GMRESB took 20 iterations GMRESB residual error history: 1 2.100000e+01 2 9.391486e+00 3 5.612486e+00 4 3.834058e+00 5 2.831639e+00 6 2.201398e+00 7 1.774824e+00 8 1.470294e+00 9 1.243933e+00 10 1.070259e+00 11 9.335639e-01 12 8.236878e-01 13 7.337994e-01 14 6.591531e-01 15 5.963599e-01 16 5.429421e-01 17 4.970501e-01 18 4.572787e-01 19 4.225429e-01 20 3.919930e-01 21 4.975583e-14 Before GMRESB, ||X_ERROR|| = 5.177254e+01 Before GMRESB, ||A*x-b|| = 2.001571e+01 After GMRESB, ||X_ERROR|| = 2.207147e-13 After GMRESB, ||A*x-b|| = 4.890277e-14 GMRESB took 20 iterations GMRESB residual error history: 1 2.001571e+01 2 8.952209e+00 3 5.379749e+00 4 3.697509e+00 5 2.770780e+00 6 2.167318e+00 7 1.747509e+00 8 1.438216e+00 9 1.212309e+00 10 1.041721e+00 11 9.023415e-01 12 7.936172e-01 13 7.069887e-01 14 6.322409e-01 15 5.705450e-01 16 5.179862e-01 17 4.758391e-01 18 4.392913e-01 19 4.069644e-01 20 3.782890e-01 21 6.352708e-14 KELLEY_TEST05 Test PCGSOL on the simple -1,2-1 matrix of order N = 20 Before PCGSOL, ||X_ERROR|| = 5.357238e+01 Before PCGSOL, ||A*x-b|| = 2.100000e+01 After PCGSOL, ||X_ERROR|| = 1.177936e-14 After PCGSOL, ||A*x-b|| = 1.822559e-14 PCGSOL took 20 iterations PCGSOL residual error history: 1 2.100000e+01 2 1.050000e+01 3 7.000000e+00 4 5.250000e+00 5 4.200000e+00 6 3.500000e+00 7 3.000000e+00 8 2.625000e+00 9 2.333333e+00 10 2.100000e+00 11 1.909091e+00 12 1.750000e+00 13 1.615385e+00 14 1.500000e+00 15 1.400000e+00 16 1.312500e+00 17 1.235294e+00 18 1.166667e+00 19 1.105263e+00 20 1.050000e+00 21 5.752886e-15 Before PCGSOL, ||X_ERROR|| = 5.149830e+01 Before PCGSOL, ||A*x-b|| = 2.160416e+01 After PCGSOL, ||X_ERROR|| = 4.398725e-14 After PCGSOL, ||A*x-b|| = 1.590683e-14 PCGSOL took 20 iterations PCGSOL residual error history: 1 2.160416e+01 2 1.025293e+01 3 6.919156e+00 4 5.269702e+00 5 4.128192e+00 6 3.382017e+00 7 2.963860e+00 8 2.594522e+00 9 2.240644e+00 10 2.031090e+00 11 1.791799e+00 12 1.682394e+00 13 1.508757e+00 14 1.395133e+00 15 1.351295e+00 16 1.259587e+00 17 1.172750e+00 18 1.113847e+00 19 1.081189e+00 20 1.018333e+00 21 5.970440e-15 KELLEY_TEST06 Test NSOL on the Chandrasekhar function. 1.000000 0.066935 0.147959 2.0000000 0.0012205 0.0182334 3.00000000000 0.00000034967 0.00028650715 0.00000 0.45239 0.45239 1.00000 0.06693 0.14796 2.00000 0.00122 0.01823 3.00000 0.00000 0.00029 Error flag IERR = 0 Residual norm: 0.000002 it_hist = 0.45238815323 0.06693491011 0.00122045336 0.00000034967 kelley_test Normal end of execution. 06-Feb-2019 08:15:29