Hi Barış,
you can for instance make the lead_onsite different from onsite and then print syst[lat(i,j)] for different regions to check. 

Le sam. 6 juil. 2019 à 22:43, Barış Özgüler <ozguler@wisc.edu> a écrit :
Hi Ousmane,

How can I make sure that there is no overlap between the lead and normal region? Lead_shape() returns -L <= x <= L and we want the lead to satisfy x <= 0. If we don't use color_site(site), then we would get the attached plot. The black region not only represents the normal region but also represents the lead. However, the lead may also cover the region of x > 0 as the function Lead_shape() suggests. Is there another way to check what part of the shape is scattering region, what part is lead?


A. Barış Özgüler

PhD Candidate in Physics


University of Wisconsin–Madison
Madison, Wisconsin 53706, USA