{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "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",
    "\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=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",
    "\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"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 1.8.5",
   "language": "julia",
   "name": "julia-1.8"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "1.8.5"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}