I am trying to generate a numerical initial condition for a PDE that looks like this: it is a Bravais lattice, with a gaussian at each site (in 2D). This is not difficult to do by hand for a square lattice with two atoms per unit cell, for example: using Plots const N = 256 dx = 1/143 dy = dx x =…

Bravais lattice with a gaussian at every site

rafael.guerra January 23, 2022, 2:24pm 3

FWIW, the packages indicated in this thread handle Bravais lattices, and there is also the PWDFT.jl package which deals with them.

1 Like