Skip to content

Instantly share code, notes, and snippets.

@mfherbst
mfherbst / lobpcg_test.jl
Created May 22, 2019 07:49
A test for the LOBPCG eigensolver from IterativeSolvers.jl
using IterativeSolvers
function run_test(;N=1000)
A = [7.44333, 7.44333, 5.77377, 6.35301, 6.35301, 6.18264, 7.44333, 5.77377,
6.35301, 5.77377, 3.35461, 3.18425, 5.26268, 6.35301, 3.18425, 2.26429,
3.59312, 7.17075, 5.26268, 3.59312, 4.17236, 7.00039, 7.17075, 7.00039,
6.35301, 6.18264, 6.35301, 3.18425, 2.26429, 3.59312, 7.17075, 6.18264,
2.26429, 0.594723, 1.17396, 4.00199, 3.59312, 1.17396, 1.0036, 3.08203,
7.40926, 7.17075, 4.00199, 3.08203, 4.41086, 7.9885, 7.40926, 7.9885,
5.26268, 3.59312, 4.17236, 7.00039, 3.59312, 1.17396, 1.0036, 3.08203,
@mfherbst
mfherbst / lobpcg_test2.jl
Last active June 3, 2019 09:36
lobpcg_test2.jl
using IterativeSolvers
function run_test(;N=1000, nev=10)
A = [4.71217, 4.59973, 4.44981, 3.96257, 4.59973, 4.44981, 4.71217, 3.58777, 3.10053,
3.73769, 4.93705, 2.95061, 2.08857, 2.35093, 3.73769, 3.43785, 2.20101, 2.08857,
3.10053, 3.43785, 2.95061, 3.58777, 4.93705, 3.58777, 2.72573, 2.98809, 4.37485,
4.93705, 2.57581, 1.33897, 1.22653, 2.23849, 4.37485, 2.68825, 1.07661, 0.589372,
1.22653, 2.98809, 3.92509, 1.93865, 1.07661, 1.33897, 2.72573, 3.92509, 2.68825,
2.57581, 3.58777, 4.93705, 4.71217, 3.47533, 3.36289, 4.37485, 3.32541, 1.71377,
1.22653, 1.86369, 3.62525, 3.06305, 1.07661, 0.214573, 0.476932, 1.86369, 4.37485,
Ecut = 3.0
Ecut = 4.0
Ecut = 5.0
^C
signal (2): Unterbrechung
in expression starting at /thread_reproduce.jl:13
do_futex_wait.constprop.0 at /lib/x86_64-linux-gnu/libpthread.so.0 (unknown line)
__new_sem_wait_slow.constprop.0 at /lib/x86_64-linux-gnu/libpthread.so.0 (unknown line)
lock_planner_mutex at /home/mfherbst/.julia/artifacts/1f7802d864a42a8958dafb40427722f8a7abbe07/lib/libfftw3.so (unknown line)
fftw_destroy_plan at /home/mfherbst/.julia/artifacts/1f7802d864a42a8958dafb40427722f8a7abbe07/lib/libfftw3.so (unknown line)