�
�F�V( � @ sT d d l Z d d l Z d d l m Z d d l m Z d d l Z d d l Z d d l m Z e
d d � e j j d � d d � D� � d k r� e
d � � n d d
d d g Z e d e j d d g d d g g � d e j d d g d d g g � d e j d d g d d g g � d e j d d g d d g g � � Z e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _ e j e j e j � e _! e j e j e j � e _" e j e j e j � e _# e j e j e j � e _$ e j e j e j � e _% d d d d d d d d d d d d d d d d � Z& d d � Z' d d d d d d d
� Z( d S) � N)�copy)�SimpleNamespace)�
wraparoundc c s | ] } t | � Vq d S)N)�int)�.0�i� r �../code/functions.py� <genexpr> s r
�.� � � z numpy >= (1, 8, 0) is required�spectrum�hamiltonian_array�h_k�pauli�s0g �?g �sx�syy �?�szc sj t j d f d t j d f g } | d k r: t � } n | j j � | | | g } d d � | D� } � f d d � t | � D� } | \ } } } t | | | | | d � \ } } t | � d k r� t j j | � � n t | � d k r� t
d � � ngt | � d
k r�| d k rb| d d d
j � k rQd j | d d � } qb| d d } n | d k rwd } n t
j | d d
� f d | | g �} i } t | t j � r�t | � | d <n | d k r�n
| | d <t |
t j � r t |
� | d <n t |
t � r%|
| d <n | d k r=t | � n t d � } |
d k rat |
� n t d � }
t | � r�| j | | � � } n | | |
f d i | d 6� St | � d k rZi � | d k r�| d d d
j � k r�| � d <n |
d k r-| d
d d
j � k r-| � d <n | d k r~| d d d
j � k rmd j | d d � } q~| d d } n | d k r�| d
d d
j � k r�d j | d
d � } q�| d
d } n | d k r�d } n | d k r�| � d <n | d k r@t j t | d d
� t | d d
� g d � } n |
d k r�t j t | d d
� t | d d
� g d � }
n | d k r�d } n i | d |
d | d | d
|
d
| d
f d 6| | g d 6| g d 6� t
j � � � f d d � t � j d � D� � } t | � r?| j | | � � } n | d i i d d 6d 6� St
d � � d S)!ah Function that plots system spectrum for varying parameters or momenta.
Parameters:
-----------
sys : kwant.Builder object
The un-finalized (in)finite system.
p : SimpleNamespace object
A container used to store Hamiltonian parameters. The parameters that
are sequences are used as plot axes.
k_x, k_y, k_z : floats or sequences of floats
Real space momenta at which the Hamiltonian has to be evaluated.
If the system dimensionality is low, extra momenta are ignored.
title : function
Function that takes p as argument and generates a string.
xdim, ydim, zdim : holoviews.Dimension or string
The labels of the axes. Default to best guess, and extra ones
are ignored.
xticks, yticks zticks : list
Lists of axes xticks, extra ones are ignored.
xlims, ylims, zlims : tuple
Upper and lower plot limit of the axes. If None the upper and lower
limits of the ticks are used. Extra ones are ignored.
Returns:
--------
plot : holoviews.Path object
Plot of varying parameter vs. spectrum.
z$-\pi$r �$0$z$\pi$Nc S s> g | ]4 } | d k r4 t j t j t j d � n | � q S)N�e )�np�linspace�pi)r r r r r �
<listcomp>K s z spectrum.<locals>.<listcomp>c s. g | ]$ \ } } | � k r$ | n d � q S)r r )r �jr )�dimensionalityr r r M s Tr
� zA 0D plot requestedz k_x k_y k_zz${}$z$E$�kdims�xticks�yticks�plot�Path�zticks�extents�vdimsc sD g | ]: } t j � d d � d d � | f � � d � � � q S)Nr# )�hv�Surface)r r )�energies�kwargs�styler r r � s �� �fig_size�OverlayzCannot make 4D plots yet.)r r )r
r )NN�����) r r r �symmetry�num_directions� enumerater �len�linalg�eigvalsh�
ValueError�split�formatr( r$ �
isinstance� collections�Iterable�listr �slice�callable�relabel�round�min�maxr/ �range�shape)�sys�p�k_x�k_y�k_z�title�xdim�ydim�zdimr! r" r% ZxlimsZylimsZzlimsZpi_ticks�k� hamiltonians� variablesr# �ticksr )r r* r+ r, r r ' s� "
! )
$$ &
&
7 7
c C s t | | | � d S)zNFunction that returns the Hamiltonian of a kwant 1D system as a momentum.
r )r )rF rG Zmomentumr r r r � s Fc
s� | d k r t � } n y �
j j j d � Wn t k
rI d � Yn X�
j j � t | � � � d k r� �
j � �
d d � � n� t �
j j � d k r� � f d d � � nT t j
�
j j � j } | j t j
j | j j | � � � � � � f d d � � t �
� j � �
t � } x< � j j � D]+ \ } } t | t j � r5| | | <q5q5Wxr d � f d � f d
� f g D]O \ } } | | k r�t d j | � � � n t | t j � r�| | | <q�q�W| s%�
j � g � � � � g � d d
�d }
| r |
g f S|
Sn � � � � � �
f d d � � t t | j � � � \ � } � � f d d � t j | � D� }
t |
d j � } t j
|
� j d d � | D� | � }
| r�|
t t � | � � f S|
Sd S)a4 Evaluate the Hamiltonian of a system over a grid of parameters.
Parameters:
-----------
sys : kwant.Builder object
The un-finalized kwant system whose Hamiltonian is calculated.
p : SimpleNamespace object
A container of Hamiltonian parameters. The parameters that are
sequences are used to loop over.
k_x, k_y, k_z : floats or sequences of floats
Momenta at which the Hamiltonian has to be evaluated. If the system
only has 1 translation symmetry, only `k_x` is used, and interpreted as
lattice momentum. Otherwise the momenta are in reciprocal space.
return_grid : bool
Whether to also return the names of the variables used for expansion,
and their values.
Returns:
--------
hamiltonians : numpy.ndarray
An array with the Hamiltonians. The first n-2 dimensions correspond to
the expanded variables.
parameters : list of tuples
Names and ranges of values that were used in evaluation of the
Hamiltonians.
Examples:
---------
>>> hamiltonian_array(sys, SimpleNamespace(t=1, mu=np.linspace(-2, 2)),
... k_x=np.linspace(-np.pi, np.pi))
>>> hamiltonian_array(sys_2d, p, np.linspace(-np.pi, np.pi),
... np.linspace(-np.pi, np.pi))
Nr
r c S s g S)Nr )rO r r r �momentum_to_lattice� s z.hamiltonian_array.<locals>.momentum_to_latticec s0 t | � d � � r% t d � � n | d g S)Nz0Dispersion is 1D, but more momenta are provided.r )�anyr7 )rO )r r r rS � s c sc t j j � | d � � � d d � \ } } t j t | � d k � rY t d � � n t | � S)Nr gH�����z>z>Requested momentum doesn't correspond to any lattice momentum.)r r5 �lstsqrT �abs� RuntimeErrorr= )rO � residuals)�A�space_dimensionalityr r rS � s /rH rI rJ zUOne of the system parameters is {}, which is reserved for momentum. Please rename it.�sparseF.c sl � j j | � | j d � � | j d � � | j d � � g } � | � } � j d � g | d d � S)NrH rI rJ �argsr[ F)�__dict__�update�get�hamiltonian_submatrix)�valuesrO )rH rI rJ rS �parsrF r r � hamiltonian s
z&hamiltonian_array.<locals>.hamiltonianc s+ g | ]! } � t t � | � � � � q Sr )�dict�zip)r �value)rc �namesr r r s z%hamiltonian_array.<locals>.<listcomp>c S s g | ] } t | � � q Sr )r4 )r rf r r r r s r0 )N.) r r1 �periodsrE �AttributeErrorr2 r � finalizedr4 r �array�T�dotr5 �invr rd r] �itemsr: r; r< rW r9 r` re �sorted� itertools�productr= �reshape)
rF rG rH rI rJ Z return_grid�BZchanging�keyrf rP ra �sizer ) rY r rc rH rI rJ rS rg rb rZ rF r r � sX #
' .
! )r
r r y � �g �))rq r; r �typesr �numpyr � holoviewsr( r �tuple� __version__r8 rW �__all__rk r �kronr Zs0s0r Zs0sxr Zs0syr Zs0szZsxs0ZsxsxZsxsyZsxszZsys0ZsysxZsysyZsyszZszs0ZszsxZszsyZszszr r r r r r r �<module> sD 5'!!$� 