Current state of the simulation

This commit is contained in:
dario
2024-04-04 15:22:50 +02:00
committed by dario
parent 828af75fe3
commit dedd445470
4 changed files with 49 additions and 15 deletions

View File

@@ -13,15 +13,25 @@ class GaussianNoise(Transform):
self.__mu = mu
self.__sigma = sigma
def __call__(self, t: float, x: ArrayLike) -> ArrayLike:
assert np.shape(self.__mu) == np.shape(x), "Mu and x have to match in shape."
def __call__(self, _: float, x: ArrayLike) -> ArrayLike:
if np.isscalar(x):
noise = np.random.normal(0, 1)
x += self.__sigma * noise + self.__mu
else:
dim = len(x)
if np.isscalar(self.__sigma):
sigma = np.identity(dim) * self.__sigma
else:
sigma = self.__sigma
if np.isscalar(self.__mu):
mu = np.ones(dim)
else:
mu = self.__mu
noise = np.random.normal(0, 1, np.shape(x))
x += self.__sigma @ noise + self.__mu
x += sigma @ noise + mu
return x