Quasiperiodic TTG: wave functions, conductivity and superconductivity

This is the code to calculate the wave functions, conductivity and superconductivity in quasiperiodic twisted trilayer graphene. The single-particle wave function is calculated with the Chebyshev filter and shake method The normal state conductivity is calculated with kernel polynomial method (KPM) expansion of the Kubo formula The superconductivity order parameter is calculated with KPM-based self-consistent BdG, while the superfluid stiffness is calculated with double KPM method. The code is written in fortran and the makefile is alongside with the code. The default compiler we use is the intel fortran compiler and we take advantange of the MKL library functions for (small) matrix diagonalizations. The code can be compiled with other compilers, e.g. gfortran, with minimal modificiton by replacing the MKL library functions.