Complex Nuc1(d3vec G,double norm) Complex Nuc2(d3vec G,double norm)Return values:
- Nuc1: norm*exp(-i*(G.x*r1x+G.y*r1y+G.z*r1z)), where (r1x,r1y,r1z)=(0,0,0)
- Nuc2: norm*exp(-i*(G.x*r2x+G.y*r2y+G.z*r2z)), where (r2x,r2y,r2z)=(1.5,0,0)
Input:
- G: G vector
- norm: input parameter (usually will be set to Z/Vol, where Vol is the volume of the cell and Z is the charge on the nucleus.
Provide the functions of the above two prototypes to be able to initialize J(rhonuc) for nuclei at two different locations.