def __ep(self, x):
d = np.linalg.norm(x - self.coordinates, axis=1)
v = pe_harmonic(d, 0.0, 0, self.__P, self.__B) * self.__A
return v
def __eval_ep_harmonic(self, x):
z, d = get_angle(x, center=self.coordinates, ref_frame=self.coordinate_system)
v = pe_harmonic(d, z, self.__l, self.__P, self.__B) * self.__A
return v
r = np.arange(1.0, 5.0, 0.1)
B = 1.0
l = 1.0
P = 0.0
sg1 = short_generalized(r, B, l, P)
sg2 = long_generalized(r, B, l, P)
# CASE 1
r = 2.25
P = 1.0
l = 1
B = 2.24
t = np.pi / 3
V = pe_harmonic(r, t, l, P, B)
print("r={:3.2f} t={:3.2f} P=1, l=1, B={:3.2f} -> V={:18.12f}".format(
r, t, B, V))
r = 2.25
P = 0.0
l = 2
B = 2.24
t = np.pi / 3
V = pe_harmonic(r, t, l, P, B)
print("r={:3.2f} t={:3.2f} P=0, l=2, B={:3.2f} -> V={:18.12f}".format(
r, t, B, V))
r = 2.25
P = 1