// Qu' is 0 except for one value at position i, so it's really (Qui')*Yi
float[] QuiYi = Yi.clone();
for (int i = 0; i < QuiYi.length; i++) {
QuiYi[i] *= targetQui;
}
newXu = YTYsolver.solveFToD(QuiYi);
}
}
// Similarly for Y vs X
double[] newYi = null;