Задача из квантовой механики систем конечной размерности
...
так что бы точность более менее равномерно распределялась по всем собственным значениям
имхо сугубо на своем опыте, а опыт основан на том, что когда-то писал диссер именно по этой теме и размерности моих матриц зашкаливали за миллиард неизвестных, сильно помогает придумать хороший предобуславливатель и решать эту задачу блочным Давидсоном.
Более детально наверное смогу рассказать, когда поподробнее расскажете что именно Вам хочется считать, так как там ньюансов реально много даже начиная от того, что скалярное произведение при больших размерностях просто валится в двойной точности если в лоб, а без скалярного все крыловские и предобусловленные крыловские методы совсем не работают.