From 0f53e4e71990429037d3a7274014ee21d30c12f7 Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Sun, 19 Mar 2023 14:46:00 -0400
Subject: [PATCH] Identified two discrepancies

---
 tester.ipynb | 76 +++++++++++++++++++++++++++++++---------------------
 1 file changed, 45 insertions(+), 31 deletions(-)

diff --git a/tester.ipynb b/tester.ipynb
index 2878b6c..d05e40a 100644
--- a/tester.ipynb
+++ b/tester.ipynb
@@ -6,7 +6,45 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# ./En.run -d 3 -n 3 -e 5 -c eps=0 -c pot=v_gauss,v0=-4,r=2 -N 6 -L 5:14 -c n_imag=1\n",
+    "# ./En.run -d 3 -n 3 -e 5 -c pot=v_gauss,v0=-4,r=2 -N 8 -L 4 -c n_imag=0\n",
+    "# Calculating...\n",
+    "# -> N = 8\n",
+    "# -> L = 4\n",
+    "# -> Spectrum = {-6.07632,-3.81486,-3.71969,-3.71968,-3.38263}...\n",
+    "# Done.\n",
+    "# -> Time used = 12s\n",
+    "\n",
+    "include(\"CPU.jl\")\n",
+    "\n",
+    "T=Float32\n",
+    "\n",
+    "function V_test(r2::T)::T\n",
+    "    return -4*exp(-r2/4)\n",
+    "end\n",
+    "\n",
+    "N=8\n",
+    "L=4.0f0\n",
+    "n_imag=0\n",
+    "\n",
+    "H=HOperator{T}(V_test,3,3,N,L,0.5f0,n_imag)\n",
+    "@time evals,evecs,info=eig(H,5)\n",
+    "print(info.numops,\" operations : \")\n",
+    "println(evals)\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# ./En.run -d 3 -n 3 -e 5 -c pot=v_gauss,v0=-4,r=2 -N 6 -L 14 -c n_imag=0\n",
+    "# Calculating...\n",
+    "# -> N = 6\n",
+    "# -> L = 14\n",
+    "# -> Spectrum = {-8.49538,-3.35492,-3.34356,-3.32830,-3.07909}...\n",
+    "# Done.\n",
+    "# -> Time used = 0.229049s\n",
     "\n",
     "include(\"CPU.jl\")\n",
     "\n",
@@ -17,37 +55,13 @@
     "end\n",
     "\n",
     "N=6\n",
-    "for L::T in 5.0:14.0\n",
-    "    H=HOperator{T}(V_test,3,3,N,L,0.5f0,1)\n",
-    "    @time evals,evecs,info=eig(H,5)\n",
-    "    print(info.numops,\" operations : \")\n",
-    "    println(evals)\n",
-    "end"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "# ./En.run -d 3 -n 2 -e 5 -c eps=0 -c pot=v_gauss,v0=-4,r=2 -N 32 -L 5:14 -c n_imag=0\n",
+    "L=14.0f0\n",
+    "n_imag=0\n",
     "\n",
-    "include(\"CPU.jl\")\n",
-    "\n",
-    "T=Float64\n",
-    "\n",
-    "function V_test(r2::T)::T\n",
-    "    return -4*exp(-r2/4)\n",
-    "end\n",
-    "\n",
-    "N=32\n",
-    "for L::T in 5.0:14.0\n",
-    "    H=HOperator{T}(V_test,3,2,N,L)\n",
-    "    @time evals,evecs,info=eig(H,5)\n",
-    "    print(info.numops,\" operations : \")\n",
-    "    println(evals)\n",
-    "end"
+    "H=HOperator{T}(V_test,3,3,N,L,0.5f0,n_imag)\n",
+    "@time evals,evecs,info=eig(H,5)\n",
+    "print(info.numops,\" operations : \")\n",
+    "println(evals)"
    ]
   }
  ],