@Reptorian
Looks like we make the same kind of diagonals on quantum_torus_small.png @ N=30. They diff to zero. I like your diagonalizer the best. I think we can check off the diag() implementation on the todo. Got to work a bit on @David_Tschumperle -meigen to pull the eigen_vectors from the math expression environment through to the command line. Or use the math expression parser with -exp{math_lib} and pull from meig(A,m,n) via the math expression’s store().
Don’t know when I can pick that up; it’s up for grabs. Hopefully can pick it up toward the end of the week. Stay tuned… Almost home…
gosgood@bertha ~/git_repositories/SchrodingerWellPython $ gmic osg_hamiltonian.cimg rep_hamiltonian.cimg -sub[-2,-1][gmic]-0./ Start G'MIC interpreter.
[gmic]-0./ Input file 'osg_hamiltonian.cimg' at position 0 (1 image 900x900x1x1).
[gmic]-1./ Input file 'rep_hamiltonian.cimg' at position 1 (1 image 900x900x1x1).
[gmic]-2./ Subtract images [0,1].
[gmic]-1./ Display image [0] = 'osg_hamiltonian.cimg'.
[0] = 'osg_hamiltonian.cimg':
size = (900,900,1,1) [3164 Kio of floats].
data = (0,0,0,0,0,0,0,0,0,0,0,0,(...),0,0,0,0,0,0,0,0,0,0,0,0).
min = 0, max = 0, mean = 0, std = 0, coords_min = (0,0,0,0), coords_max = (0,0,0,0).
[gmic]-1./ End G'MIC interpreter.