Testing for v1.0.2 of SIAMFANLEquations.jl

Author: ctkelley

SIAMFANL.ipynb

@@ -678,7 +678,7 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.0"
+   "version": "1.10.2"
   }
  },
  "nbformat": 4,

SIAMFANLCh1.ipynb

@@ -2004,7 +2004,7 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.0"
+   "version": "1.10.2"
   }
  },
  "nbformat": 4,

SIAMFANLCh1s.ipynb

@@ -745,10 +745,10 @@
        "  ────────────────────────────────────────────────────────────────────────────\n",
        "\n",
        "\u001b[1m  Example for nsoli\u001b[22m\n",
-       "\u001b[1m  –––––––––––––––––––\u001b[22m\n",
+       "\u001b[1m  –––––––––––––––––\u001b[22m\n",
        "\n",
        "\u001b[1m  Simple 2D problem.\u001b[22m\n",
-       "\u001b[1m  --------------------\u001b[22m\n",
+       "\u001b[1m  ------------------\u001b[22m\n",
        "\n",
        "  You should get the same results as for nsol.jl because GMRES will solve the\n",
        "  equation for the step exactly in two iterations. Finite difference Jacobians\n",
@@ -896,7 +896,7 @@
        " 1.55516e-04  1.84641e-03  1.84641e-03\n",
        " 1.33169e-06  1.82363e-04  1.82363e-04\n",
        " 1.13965e-08  2.34291e-06  2.34291e-06\n",
-       " 9.75259e-11  2.42602e-11  2.42602e-11"
+       " 9.75259e-11  2.42596e-11  2.42596e-11"
       ]
      },
      "execution_count": 5,
@@ -945,11 +945,11 @@
      "output_type": "stream",
      "text": [
       "Shamanskii, n=5\n",
-      "  94.777 ms (74 allocations: 739.66 KiB)\n",
+      "  95.012 ms (75 allocations: 740.95 KiB)\n",
       "Newton-GMRES, fixed eta\n",
-      "  1.185 ms (303 allocations: 513.23 KiB)\n",
+      "  1.040 ms (319 allocations: 513.47 KiB)\n",
       "Newton-GMRES, Eisenstat-Walker\n",
-      "  1.187 ms (303 allocations: 513.23 KiB)\n"
+      "  1.041 ms (319 allocations: 513.47 KiB)\n"
      ]
     }
    ],
@@ -983,9 +983,9 @@
      "output_type": "stream",
      "text": [
       "Newton-GMRES, fixed eta\n",
-      "  1.228 ms (370 allocations: 509.06 KiB)\n",
+      "  1.058 ms (385 allocations: 508.58 KiB)\n",
       "Newton-GMRES, Eisenstat-Walker\n",
-      "  1.229 ms (370 allocations: 509.06 KiB)\n"
+      "  1.056 ms (385 allocations: 508.58 KiB)\n"
      ]
     }
    ],
@@ -1060,7 +1060,7 @@
      "output_type": "stream",
      "text": [
       "nsol, sham=5\n",
-      "  3.571 ms (416 allocations: 6.29 MiB)\n"
+      "  3.631 ms (421 allocations: 6.30 MiB)\n"
      ]
     }
    ],
@@ -1099,7 +1099,7 @@
      "output_type": "stream",
      "text": [
       "nsoli, not preconditioned\n",
-      "  4.380 ms (3899 allocations: 986.98 KiB)\n"
+      "  3.645 ms (4259 allocations: 1004.08 KiB)\n"
      ]
     }
    ],
@@ -1132,7 +1132,7 @@
      "text": [
       "nsoli, preconditioned, \n",
       "    Eisenstat-Walker forcing term\n",
-      "  1.377 ms (772 allocations: 653.78 KiB)\n"
+      "  1.334 ms (792 allocations: 654.31 KiB)\n"
      ]
     }
    ],
@@ -1162,7 +1162,7 @@
      "output_type": "stream",
      "text": [
       "nsoli, preconditioned, fixed eta\n",
-      "  1.729 ms (975 allocations: 927.78 KiB)\n"
+      "  1.684 ms (1007 allocations: 928.69 KiB)\n"
      ]
     }
    ],
@@ -1193,7 +1193,7 @@
      "text": [
       "nsoli, preconditioned, \n",
       "    Eisenstat-Walker forcing term\n",
-      "  1.364 ms (792 allocations: 741.77 KiB)\n"
+      "  1.329 ms (816 allocations: 742.42 KiB)\n"
      ]
     }
    ],
@@ -1225,7 +1225,7 @@
      "text": [
       "nsoli, left preconditioned, \n",
       "    Eisenstat-Walker forcing term\n",
-      "  1.555 ms (896 allocations: 747.52 KiB)\n"
+      "  1.396 ms (990 allocations: 751.45 KiB)\n"
      ]
     }
    ],
@@ -1258,7 +1258,7 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "  2.438 ms (1634 allocations: 2.05 MiB)\n"
+      "  2.317 ms (1806 allocations: 2.06 MiB)\n"
      ]
     }
    ],
@@ -1437,8 +1437,7 @@
        "\n",
        "julia> function PreCondBeam(v, x, bdata)\n",
        "          J = bdata.D2\n",
-       "          ptv = J\u000b",
-       "\n",
+       "          ptv = J\\v\n",
        "       end\n",
        "PreCondBeam (generic function with 1 method)\n",
        "\n",
@@ -1626,8 +1625,7 @@
        "\n",
        "julia> function PreCondBeam(v, x, bdata)\n",
        "          J = bdata.D2\n",
-       "          ptv = J\u000b",
-       "\n",
+       "          ptv = J\\v\n",
        "       end\n",
        "PreCondBeam (generic function with 1 method)\n",
        "\n",
@@ -1863,10 +1861,10 @@
        "  iteration + 1. So, for scalar equations, it's a row vector.\n",
        "\n",
        "\u001b[1m  Example for ptcsol\u001b[22m\n",
-       "\u001b[1m  ––––––––––––––––––––\u001b[22m\n",
+       "\u001b[1m  ––––––––––––––––––\u001b[22m\n",
        "\n",
        "\u001b[1m  The buckling beam problem.\u001b[22m\n",
-       "\u001b[1m  ----------------------------\u001b[22m\n",
+       "\u001b[1m  --------------------------\u001b[22m\n",
        "\n",
        "  You'll need to use TestProblems for this to work. The preconditioner is a\n",
        "  solver for the high order term.\n",
@@ -1875,8 +1873,7 @@
        "\u001b[36m  \u001b[39m\n",
        "\u001b[36m  julia> function PreCondBeam(v, x, bdata)\u001b[39m\n",
        "\u001b[36m            J = bdata.D2\u001b[39m\n",
-       "\u001b[36m            ptv = J\u000b",
-       "\u001b[39m\n",
+       "\u001b[36m            ptv = J\\v\u001b[39m\n",
        "\u001b[36m         end\u001b[39m\n",
        "\u001b[36m  PreCondBeam (generic function with 1 method)\u001b[39m\n",
        "\u001b[36m  \u001b[39m\n",
@@ -1986,7 +1983,7 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "  575.041 μs (571 allocations: 2.75 MiB)\n"
+      "  517.959 μs (490 allocations: 2.47 MiB)\n"
      ]
     }
    ],
@@ -2036,7 +2033,7 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "  1.847 ms (1891 allocations: 6.01 MiB)\n"
+      "  1.632 ms (1945 allocations: 6.11 MiB)\n"
      ]
     }
    ],
@@ -2125,7 +2122,7 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.0"
+   "version": "1.10.2"
   }
  },
  "nbformat": 4,

SIAMFANLCh2.ipynb

@@ -46,7 +46,7 @@ fpsize=fsize(printlabel)
 figure(1)
 subplot(1,2,1)
 ~printlabel || title("Figure 1.6 from print book")
-semilogy(pival,abs.(ptcerr),"k--",nival,abs.(newterr),"k-")
+semilogy(pival,abs.(ptcerr).+1.e-18,"k--",nival,abs.(newterr),"k-")
 PTClabel=L"$\Psi$TC"
 legend((PTClabel,"Newton"),fontsize=fpsize)
 xlabel("Iterations",fontsize=fpsize)