63 lines
1.2 KiB
Plaintext
63 lines
1.2 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"using Plots\n",
|
|
"include(\"helper.jl\")\n",
|
|
"include(\"p_space.jl\")"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"vertices = [0, 1 - 0.5im, 2, 6]\n",
|
|
"subdivisions = [64, 64, 64]\n",
|
|
"p, w = get_mesh(vertices, subdivisions)\n",
|
|
"\n",
|
|
"scatter(real.(p), imag.(p))"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"# ComplexScaling-FV: Eq. (54)\n",
|
|
"V_pq(p, q) = -10 * g1(1, p, q)\n",
|
|
"\n",
|
|
"H = get_H_matrix(V_pq, p, w)\n",
|
|
"evals = eigen(H).values\n",
|
|
"\n",
|
|
"E_target = 0.258 - 0.164im\n",
|
|
"E = nearest(evals, E_target)\n",
|
|
"\n",
|
|
"print(\"E = $E\")\n",
|
|
"scatter(real.(evals), imag.(evals), xlim = (0,1))"
|
|
]
|
|
}
|
|
],
|
|
"metadata": {
|
|
"kernelspec": {
|
|
"display_name": "Julia 1.9.0",
|
|
"language": "julia",
|
|
"name": "julia-1.9"
|
|
},
|
|
"language_info": {
|
|
"file_extension": ".jl",
|
|
"mimetype": "application/julia",
|
|
"name": "julia",
|
|
"version": "1.9.0"
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 2
|
|
}
|