Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

#

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

Jump to bottom

Visualization #48

Merged
merged 3 commits into from
Nov 5, 2022
Merged

Visualization #48

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
d14d296
add 2D irreducible plot
fsxbhyy Nov 4, 2022
2b8ece4
merge
fsxbhyy Nov 4, 2022
5f3e41d
Merge commit 'fdd8' into visualization
fsxbhyy Nov 4, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
21 changes: 20 additions & 1 deletion toolbox/plotlyjs_wignerseitz.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -364,7 +364,7 @@ function plot(c::Cell{2}, layout::Layout=Layout();
return PlotlyJS.plot(ts, layout; config=config)
end

function plot(clist::Array{Cell{2}}, idx_center::Int, layout::Layout=Layout();
function plot(clist::Array{Cell{2}}, idx_center::Int, layout::Layout=Layout(); ibz =nothing,
config::PlotConfig=PlotConfig(responsive=true, displaylogo=false))

layout = merge(DEFAULT_PLOTLY_LAYOUT_2D, layout)
Expand All @@ -380,6 +380,25 @@ function plot(clist::Array{Cell{2}}, idx_center::Int, layout::Layout=Layout();
max_x, max_y = maximum(v -> abs(v[1]), basis(c_center)), maximum(v -> abs(v[2]), basis(c_center))
get!(layout[:xaxis], :range, [-max_x - scale / 15, max_x + scale / 15])
get!(layout[:yaxis], :range, [-max_y - scale / 15, max_y + scale / 15])


if !isnothing(ibz)
setting(ibz) !== CARTESIAN && (ibz = cartesianize(ibz))
# merge_coplanar!(ibz) # Coplanar triangles have to be merged after calling mesh3d.
# BZ
merge_coplanar!(ibz)
tbz = Vector{GenericTrace{Dict{Symbol,Any}}}(undef, length(ibz))
for (i, poly) in enumerate(ibz)
tbz[i] = PlotlyJS.scatter(
x=push!(getindex.(poly, 1), poly[1][1]),
y=push!(getindex.(poly, 2), poly[1][2]);
mode="lines", hovertext="Cell", hoverinfo="text+x+y",
line=attr(color=BZ_COL[], width=3 ), fill = "toself", fillcolor = "rgb(50, 200, 200)"
)
end
push!(tbz_list, tbz)
end

for (ci, c) in enumerate(clist)
if (!isnothing(c))
# Cell boundaries
Expand Down
26 changes: 15 additions & 11 deletions toolbox/test_extend.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,23 +106,22 @@ end

# Wigner-Seitz cells visualization
# 2D
#Rs = [[1.0, 0.0], [-0.5, √3/2]]
#N, DIM = 4, 2
#origin = [0.0 0.0]
#lattice =Matrix([2 0; 1 sqrt(3)]')
#msize = (3,3)
lattice =Matrix([2 0; 1 sqrt(3)]')
#br = BZMeshes.Brillouin(lattice = lattice)

# latvec = [1 0; 0 1]'
# 3D

lattice = [[0 0.5 0.5]; [0.5 0 0.5]; [0.5 0.5 0.0]]
#lattice = [[0 0.5 0.5]; [0.5 0 0.5]; [0.5 0.5 0.0]]
#atoms = [1,1]
#positions = [ones(3) / 8, -ones(3) / 8]

#lattice = [[1.0 0.0 0.0]; [0.0 1.0 0.0]; [0.0 0.0 1.0]]
atoms = [1]
positions = [zeros(3)]
#positions = [zeros(3)]
positions = [zeros(2)]

atom_pos = hvcat(size(positions, 1), positions...)
ibzformat = "convex hull"
Expand All @@ -132,13 +131,15 @@ convention = "ordinary"

br = BZMeshes.Brillouin(lattice=lattice, atoms=atoms, positions=positions)
#br = BZMeshes.Brillouin(lattice=lattice)
msize = (4, 4, 4)

bzmesh = UniformBZMesh(br=br, size=msize, shift=[true, true, true])
msize = (3,3)
bzmesh = UniformBZMesh(br=br, size=msize, shift=[false, false])

#msize = (4, 4, 4)
#bzmesh = UniformBZMesh(br=br, size=msize, shift=[true, true, true])
meshmap = MeshMap(bzmesh)

recip_lattice = lattice_vector(bzmesh)
#latvec = mapslices(x -> [x],recip_lattice , dims=1)[:]
latvec = mapslices(x -> [x], recip_lattice, dims=1)[:]

#generate irreducible brillouin zone
Expand Down Expand Up @@ -184,8 +185,11 @@ end

P = plot(clist, idx_center, ibz=c)

addtraces!(P, scatter3d(x=[r[1] for r in fullmesh], y=[r[2] for r in fullmesh], z=[r[3] for r in fullmesh], mode="markers", marker=attr(size=3)))
addtraces!(P, scatter3d(x=[r[1] for r in reducedmesh], y=[r[2] for r in reducedmesh], z=[r[3] for r in reducedmesh], mode="markers", marker=attr(size=3)))
addtraces!(P, scatter(x=[r[1] for r in fullmesh], y=[r[2] for r in fullmesh], mode="markers", marker=attr(size=5)))
addtraces!(P, scatter(x=[r[1] for r in reducedmesh], y=[r[2] for r in reducedmesh], mode="markers", marker=attr(size=5)))

#addtraces!(P, scatter3d(x=[r[1] for r in fullmesh], y=[r[2] for r in fullmesh], z=[r[3] for r in fullmesh], mode="markers", marker=attr(size=3)))
#addtraces!(P, scatter3d(x=[r[1] for r in reducedmesh], y=[r[2] for r in reducedmesh], z=[r[3] for r in reducedmesh], mode="markers", marker=attr(size=3)))

n = length(P.plot.data) # number of traces, the last two are the full and reduced meshes
allinvis, vis1, vis2 = [true for i in 1:n], [true for i in 1:n], [true for i in 1:n]
Expand Down Expand Up @@ -236,5 +240,5 @@ d = Dict(
relayout!(P, d)

display(P)
# readline()
readline()