Exemplo n.º 1
0
def test_dWdx_absolute_square_on_grid():
    W = Superpotential(3)
    x = np.array([[[1, 1, 1], [1, 1, 1], [0, 0, 0]],
                  [[1, 1, 1], [1, 1, 1], [0, 0, 0]]])
    print("expected: [[24.02876087+0.j 24.02876087+0.j 24.02876087+0.j]\
 [24.02876087+0.j 24.02876087+0.j 24.02876087+0.j]\
 [ 0.        +0.j  0.        +0.j  0.        +0.j]]")
    print("result: ", W.dWdx_absolute_square_on_grid(x))
 def get_potential_energy_density(self):
     W = Superpotential(self.N)
     ped = (1 / 4) * W.dWdx_absolute_square_on_grid(self.x)
     ped = np.real(ped)  #it is real anyway
     return ped