Velocity_Interpolation_wErrorCov not working for single turbine

[ufechner7]

The following Matlab test fails:

clear all
addpath('./WindField/Velocity_Interpolation_wErrorCov/');

% Simple linear data: speed = 5 at t=0, speed = 15 at t=10
times = [0.0, 10.0];
speeds = [5.0, 15.0];
data = [times(:), speeds(:)];  % Column-wise concatenation

% Simple 2x2 identity Cholesky: noise still generated but uncorrelated
chol = chol(eye(2));  % chol returns upper triangular by default in MATLAB

% Define a WindVelMatrix function equivalent
wind = WindVelMatrix(data, chol);

% Test 5: Single turbine (scalar iT input)
iT_scalar = 1;
vel_scalar = getWindSpeedT(wind, iT_scalar, 5.0);

assert(length(vel_scalar) == 1);

[MarcusBecker-GitHub]

The test uses a single turbine while defining a 2x2 matrix for the noise, which should be of dimension (nT x nT).
The function readCovMatrix generates this matrix based on the provided input.

If I modify the code as follows it works:

clear all
addpath('./WindField/Velocity_Interpolation_wErrorCov/');

% Simple linear data: speed = 5 at t=0, speed = 15 at t=10
times = [0.0, 10.0];
speeds = [5.0, 15.0];
data = [times(:), speeds(:)];  % Column-wise concatenation

% Simple 2x2 identity Cholesky: noise still generated but uncorrelated
WindVel.CholSig = chol(eye(1));  % chol returns upper triangular by default in MATLAB

% Define a WindVelMatrix function equivalent
WindVel.Data = data;

% Test 5: Single turbine (scalar iT input)
iT_scalar = 1;
vel_scalar = getWindSpeedT(WindVel, iT_scalar, 5.0);

assert(length(vel_scalar) == 1);

Changes made:

• chol(eye(1)) instead of chol(eye(2))
• implemented WindVel struct