Numpy Matmul erhöht die Kluft durch Nullwarnung ⇐ Python

[Anonymous]

Ich implementiere Gaußsche Prozessregression und arbeite derzeit daran, Stichproben aus dem GP zu erhalten. Dafür verwende ich Gleichung 2.19 aus Rasmussen/Williams:
, wobei das POWORGE die PoteniOrg. Die Operation von matmul erzeugt eine Kluft durch Null Warnung. Ich versuche zu verstehen, warum eine Matmul -Operation überhaupt eine Abteilung durchführen würde. Es scheint auch, dass der Fehler nur bei der Vorhersage von 500 INIS+ Punkten auftritt, was mich zu der Annahme veranlasst, dass es eine Speicherverbindung gibt. Kann jemand erklären, warum dieser Fehler auftritt und wie ich es verhindern kann?

Code: Select all

import numpy as np
from numpy.random import RandomState

np.seterr(all='raise') # to error out

SEED = RandomState(42)
JITTER = 1e-6
observed_x = np.array([0.0, 1.0, 3.0])
observed_y = np.array([1., -1.0, 3.0])

noise = SEED.randn(10, 1000)

lengthscale = 1.
k_ = lambda xi, xj: np.exp(- 1 * (xi[:, None] - xj) ** 2 / (2 * lengthscale ** 2))
K = (
k_(observed_x, observed_x)
)
K += JITTER * np.eye(len(observed_x))

x_star = np.linspace(-1, 4, 1000)
K_star = k_(x_star, observed_x)
L = np.linalg.cholesky(K)
alpha = np.linalg.solve(L.T, np.linalg.solve(L, observed_y))

mu_pred = K_star @ alpha

Umgebung info/pyproject.toml

Code: Select all

[project]
name = "gp"
version = "0.1.0"
readme = "README.md"
requires-python = ">=3.10"
dependencies = [
"ipykernel>=6.29.5",
"ipympl>=0.9.7",
"jupyter>=1.1.1",
"matplotlib>=3.10.1",
"numpy>=2.2.4",
"scikit-learn>=1.6.1",
"seaborn>=0.13.2",
"notebook==6.5.6",
]