public void estimate(final Access1D<?> x, final Access1D<?> y) { final int tmpRowDim = (int) Math.min(x.count(), y.count()); final int tmpColDim = this.size(); final PhysicalStore<ComplexNumber> tmpBody = ComplexDenseStore.FACTORY.makeZero(tmpRowDim, tmpColDim); final PhysicalStore<ComplexNumber> tmpRHS = ComplexDenseStore.FACTORY.makeZero(tmpRowDim, 1); for (int i = 0; i < tmpRowDim; i++) { ComplexNumber tmpX = ComplexNumber.ONE; final ComplexNumber tmpXfactor = ComplexNumber.valueOf((Number) x.get(i)); final ComplexNumber tmpY = ComplexNumber.valueOf((Number) y.get(i)); for (int j = 0; j < tmpColDim; j++) { tmpBody.set(i, j, tmpX); tmpX = tmpX.multiply(tmpXfactor); } tmpRHS.set(i, 0, tmpY); } final QR<ComplexNumber> tmpQR = QR.makeComplex(); tmpQR.decompose(tmpBody); this.set(tmpQR.solve(tmpRHS)); }
public static void assertEquals( final String message, final Access1D<?> expected, final Access1D<?> actual, final NumberContext context) { for (long i = 0; i < expected.count(); i++) { TestUtils.assertEquals(message + " @ " + i, expected.get(i), actual.get(i), context); } }
public void setValues(final Access1D<?> aValue) { for (int p = 0; p < myProcesses.length; p++) { myProcesses[p].setValue(aValue.doubleValue(p)); } }
public void set(final Access1D<?> someCoefficient) { final int tmpLimit = (int) Math.min(this.size(), someCoefficient.count()); for (int p = 0; p < tmpLimit; p++) { this.set(p, ComplexNumber.valueOf((Number) someCoefficient.get(p))); } }