sub min { eigenstates( any(@_) <= all(@_) ) }

sub max { eigenstates( any(@_) >= all(@_) ) }