To get started, select MATLAB Help or Demos from the Help menu. >> kelley_test 24-Jun-2007 08:18:18 KELLEY_TEST MATLAB version. Test the routines in the KELLEY library. TEST01 Test GMRES on the simple -1,2-1 matrix of order N = 20 Before GMRES, ||X_ERROR|| = 5.357238e+001 Before GMRES, ||A*x-b|| = 2.100000e+001 After GMRES, ||X_ERROR|| = 3.262243e-014 After GMRES, ||A*x-b|| = 7.971989e-015 GMRES took 20 iterations GMRES residual error history: 1 2.100000e+001 2 9.391486e+000 3 5.612486e+000 4 3.834058e+000 5 2.831639e+000 6 2.201398e+000 7 1.774824e+000 8 1.470294e+000 9 1.243933e+000 10 1.070259e+000 11 9.335639e-001 12 8.236878e-001 13 7.337994e-001 14 6.591531e-001 15 5.963599e-001 16 5.429421e-001 17 4.970501e-001 18 4.572787e-001 19 4.225429e-001 20 3.919930e-001 21 0.000000e+000 Before GMRES, ||X_ERROR|| = 5.107753e+001 Before GMRES, ||A*x-b|| = 1.989003e+001 After GMRES, ||X_ERROR|| = 1.146045e-013 After GMRES, ||A*x-b|| = 2.673956e-014 GMRES took 20 iterations GMRES residual error history: 1 1.989003e+001 2 9.044772e+000 3 5.353614e+000 4 3.716610e+000 5 2.745096e+000 6 2.138723e+000 7 1.716104e+000 8 1.426378e+000 9 1.205188e+000 10 1.039259e+000 11 8.886337e-001 12 7.777365e-001 13 6.888578e-001 14 6.195273e-001 15 5.584270e-001 16 5.070788e-001 17 4.632707e-001 18 4.279481e-001 19 3.963195e-001 20 3.722279e-001 21 2.289900e-029 TEST02 Test TFQMR on the simple -1,2-1 matrix of order N = 20 Before TFQMR, ||X_ERROR|| = 5.357238e+001 Before TFQMR, ||A*x-b|| = 2.100000e+001 After TFQMR, ||X_ERROR|| = 8.983893e-015 After TFQMR, ||A*x-b|| = 1.639676e-014 TFQMR took 20 iterations TFQMR residual error history: 1 2.100000e+001 2 5.824352e+000 3 2.757435e+000 4 1.610626e+000 5 1.056625e+000 6 7.466740e-001 7 5.557189e-001 8 4.297364e-001 9 3.422448e-001 10 2.790097e-001 11 2.299639e-001 12 1.872382e-001 13 1.515891e-001 14 1.224511e-001 15 9.877708e-002 16 7.947142e-002 17 6.354826e-002 18 5.014145e-002 19 3.840151e-002 20 2.703330e-002 21 2.425658e-015 Before TFQMR, ||X_ERROR|| = 5.174631e+001 Before TFQMR, ||A*x-b|| = 2.006799e+001 After TFQMR, ||X_ERROR|| = 1.260971e-014 After TFQMR, ||A*x-b|| = 2.198802e-014 TFQMR took 20 iterations TFQMR residual error history: 1 2.006799e+001 2 5.958199e+000 3 2.689475e+000 4 1.449473e+000 5 8.969515e-001 6 6.255194e-001 7 4.687470e-001 8 3.825813e-001 9 3.421511e-001 10 3.110358e-001 11 2.834665e-001 12 2.518393e-001 13 2.232365e-001 14 2.024511e-001 15 1.877948e-001 16 1.787440e-001 17 1.737966e-001 18 1.700736e-001 19 1.614461e-001 20 1.230231e-001 21 1.680772e-015 TEST03 Test BICGSTAB on the simple -1,2-1 matrix of order N = 20 Before BICGSTAB, ||X_ERROR|| = 5.357238e+001 Before BICGSTAB, ||A*x-b|| = 2.100000e+001 After BICGSTAB, ||X_ERROR|| = 1.198114e-008 After BICGSTAB, ||A*x-b|| = 3.981330e-010 BICGSTAB took 20 iterations BICGSTAB residual error history: 1 2.100000e+001 2 6.062178e+000 3 3.299832e+000 4 2.213594e+000 5 1.635007e+000 6 1.281864e+000 7 1.045450e+000 8 8.774544e-001 9 7.524163e-001 10 6.561685e-001 11 5.799848e-001 12 5.183222e-001 13 4.669121e-001 14 4.212527e-001 15 3.753684e-001 16 3.239107e-001 17 2.643018e-001 18 1.977284e-001 19 1.282688e-001 20 6.086974e-002 21 3.981330e-010 Before BICGSTAB, ||X_ERROR|| = 5.154221e+001 Before BICGSTAB, ||A*x-b|| = 2.095361e+001 After BICGSTAB, ||X_ERROR|| = 1.255112e-009 After BICGSTAB, ||A*x-b|| = 4.052968e-011 BICGSTAB took 20 iterations BICGSTAB residual error history: 1 2.095361e+001 2 5.878466e+000 3 3.185794e+000 4 2.132687e+000 5 1.573423e+000 6 1.240938e+000 7 1.008525e+000 8 8.408256e-001 9 7.178353e-001 10 6.265415e-001 11 5.525146e-001 12 4.926597e-001 13 4.433603e-001 14 3.985417e-001 15 3.527851e-001 16 3.007202e-001 17 2.438933e-001 18 1.805172e-001 19 1.164830e-001 20 5.605867e-002 21 4.053062e-011 TEST04 Test GMRESB on the simple -1,2-1 matrix of order N = 20 Before GMRESB, ||X_ERROR|| = 5.357238e+001 Before GMRESB, ||A*x-b|| = 2.100000e+001 After GMRESB, ||X_ERROR|| = 4.727957e-013 After GMRESB, ||A*x-b|| = 2.792404e-014 GMRESB took 20 iterations GMRESB residual error history: 1 2.100000e+001 2 9.391486e+000 3 5.612486e+000 4 3.834058e+000 5 2.831639e+000 6 2.201398e+000 7 1.774824e+000 8 1.470294e+000 9 1.243933e+000 10 1.070259e+000 11 9.335639e-001 12 8.236878e-001 13 7.337994e-001 14 6.591531e-001 15 5.963599e-001 16 5.429421e-001 17 4.970501e-001 18 4.572787e-001 19 4.225429e-001 20 3.919930e-001 21 2.792404e-014 Before GMRESB, ||X_ERROR|| = 5.137995e+001 Before GMRESB, ||A*x-b|| = 2.053184e+001 After GMRESB, ||X_ERROR|| = 7.289694e-014 After GMRESB, ||A*x-b|| = 4.603344e-014 GMRESB took 20 iterations GMRESB residual error history: 1 2.053184e+001 2 9.176457e+000 3 5.435724e+000 4 3.750841e+000 5 2.776325e+000 6 2.145081e+000 7 1.725387e+000 8 1.436599e+000 9 1.218714e+000 10 1.045945e+000 11 9.008410e-001 12 7.864122e-001 13 6.968358e-001 14 6.270186e-001 15 5.683769e-001 16 5.170964e-001 17 4.730314e-001 18 4.364870e-001 19 4.040206e-001 20 3.752478e-001 21 3.341020e-014 TEST05 Test PCGSOL on the simple -1,2-1 matrix of order N = 20 Before PCGSOL, ||X_ERROR|| = 5.357238e+001 Before PCGSOL, ||A*x-b|| = 2.100000e+001 After PCGSOL, ||X_ERROR|| = 1.253127e-014 After PCGSOL, ||A*x-b|| = 3.249826e-014 PCGSOL took 20 iterations PCGSOL residual error history: 1 2.100000e+001 2 1.050000e+001 3 7.000000e+000 4 5.250000e+000 5 4.200000e+000 6 3.500000e+000 7 3.000000e+000 8 2.625000e+000 9 2.333333e+000 10 2.100000e+000 11 1.909091e+000 12 1.750000e+000 13 1.615385e+000 14 1.500000e+000 15 1.400000e+000 16 1.312500e+000 17 1.235294e+000 18 1.166667e+000 19 1.105263e+000 20 1.050000e+000 21 2.529851e-015 Before PCGSOL, ||X_ERROR|| = 5.176373e+001 Before PCGSOL, ||A*x-b|| = 1.950752e+001 After PCGSOL, ||X_ERROR|| = 2.302529e-014 After PCGSOL, ||A*x-b|| = 3.157877e-014 PCGSOL took 20 iterations PCGSOL residual error history: 1 1.950752e+001 2 1.009476e+001 3 6.734063e+000 4 5.242227e+000 5 4.036716e+000 6 3.453586e+000 7 2.897504e+000 8 2.504013e+000 9 2.251101e+000 10 2.006504e+000 11 1.718821e+000 12 1.602310e+000 13 1.519175e+000 14 1.466792e+000 15 1.310363e+000 16 1.275856e+000 17 1.142483e+000 18 1.059585e+000 19 9.538062e-001 20 4.583517e-001 21 2.378361e-015 TEST06 Test NSOL on the Chandrasekhar function. 1.0000 0.0669 0.1480 2.0000 0.0012 0.0182 3.0000 0.0000 0.0003 0 0.4524 0.4524 1.0000 0.0669 0.1480 2.0000 0.0012 0.0182 3.0000 0.0000 0.0003 Error flag IERR = 0 Residual norm: 0.000002 it_hist = 0.4524 0.0669 0.0012 0.0000 KELLEY_TEST Normal end of execution. 24-Jun-2007 08:18:27 >>