wip

example/compositemesh.jl

@@ -14,7 +14,7 @@ br = BZMeshes.Cell(lattice=lattice, atoms=atoms, positions=positions)
 
 _dispersion(k) = -sum(cos.(k)) + 0.1
 # _dispersion(k) = dot(k, k) - π^2 / 4
-T = 0.1
+T = 0.001
 
 function dispersion(k)
     # constrain to 1st bz
@@ -39,24 +39,25 @@ function int_mesh(mesh)
     return integrate(data, mesh)
 end
 
-N = 16
+N = 8
 bound = [-π, π]
 theta = SimpleGrid.Uniform(bound, N; isperiodic=true)
 # bzmesh = BZMeshes.PolarMesh(dispersion=dispersion, anglemesh=theta, cell=br,
 #     kmax=π, Nloggrid=5, Nbasegrid=4, minterval=0.001)
 pm = BZMeshes.CompositePolarMesh(dispersion=dispersion,
     anglemesh=theta, cell=br, basegridtype=:uniform,
-    kmax=π * sqrt(2.1), Nloggrid=12, Nbasegrid=2, minterval=0.1T, N=4)
+    kmax=π * sqrt(2.4), Nloggrid=8, Nbasegrid=3, minterval=0.1T, N=4)
 
 pm2 = BZMeshes.CompositePolarMesh(dispersion=dispersion,
     anglemesh=theta, cell=br, basegridtype=:uniform,
-    kmax=π * sqrt(2.1), Nloggrid=8, Nbasegrid=4, minterval=0.5T, N=6)
+    kmax=π * sqrt(2.4), Nloggrid=10, Nbasegrid=6, minterval=0.01T, N=6)
 
 
 println("N=$(length(pm)), ", int_mesh(pm))
 println("N=$(length(pm2)), ", int_mesh(pm2))
 
-um = BZMeshes.UniformBZMesh(cell=br, size=(1000, 1000))
-um2 = BZMeshes.UniformBZMesh(cell=br, size=(2000, 2000))
+Nuni = 1000
+um = BZMeshes.UniformBZMesh(cell=br, size=(Nuni, Nuni))
+um2 = BZMeshes.UniformBZMesh(cell=br, size=(2Nuni, 2Nuni))
 println("N=$(length(um)), ", int_mesh(um))
 println("N=$(length(um2)), ", int_mesh(um2))