Hi Dominik,

Thanks.

I generated the plot for Bi calculation as attached. But there are very high dense KPOINTS. I thought I change optimize it by set the parameter iterator=range as following:

result_0 = z2pack.surface.run(

    system=system,

    surface=lambda s, t: [0, s / 2, t],

    iterator=range(10,24,2),

    save_file='./results/res_0.p',

    load=True

But it generate same result.

Or do I need to set s limit by pos_tol?

Regards

Wei Li