A few changes are needed within the classdef for eigenstate calculations
- Check treatment of sparsity in lines 72-75 of setup.m
- Check setup of potential matrices in lines 68-89 of setup.m
- Set cutoff to default value, e.g. the machine epsilon in the classdef eigen.m
- Check calculation of density of Hamiltonian matrix in diag.m
- Change from option 'sm' (smallest magitude) to option 'sr' (smallest real) when using sparse eigensolver eigs withinin diag.m.
- Another option: compute eigenvalues near a numeric sigma value that is nearly equal to an eigenvalue
Diff:
Diff:
Apparently, Octave has problems with the sigma value. For example Demos\ExcitonPhonon\TISE\Cyclic\N=03 is crashing