“…The GL energy functional introduced above was extensively studied for constant or regular magnetic fields in smooth domains (see, e.g., [CDR,CR2,CR3,CR4,FK1] for the 2D case and [AG,FK2,FKP,FMP,P] for the 3D one) and in domains with corners at the boundary [Cor,CG1,CG2,Gia,HK]. More recently, 2D GL models with piecewise-constant magnetic fields were also considered in [Ass1,Ass2,AK,.…”