@Test(expectedExceptions = IllegalArgumentException.class) public void testOutOfRange3() { RANGE_LIMITS.transformGradient(-3); }
@Test(expectedExceptions = IllegalArgumentException.class) public void testOutOfRange4() { UPPER_LIMIT.transformGradient(1.01); }
@Test(expectedExceptions = IllegalArgumentException.class) public void testOutOfRange4() { RANGE_LIMITS.transformGradient(1.01); }
@Test(expectedExceptions = IllegalArgumentException.class) public void testOutOfRange3() { LOWER_LIMIT.transformGradient(-3); }
/** * Calculates the Jacobian of the transform from function parameters to fitting parameters - * the i,j element will be the partial derivative of i^th fitting parameter with respect. * to the j^th function parameter * @param functionParameters The function parameters * @return matrix of partial derivative of fitting parameter with respect to function parameters */ // TODO not tested @Override public DoubleMatrix jacobian(DoubleArray functionParameters) { ArgChecker.notNull(functionParameters, "function parameters"); ArgChecker.isTrue(functionParameters.size() == _nMP, "functionParameters wrong dimension"); double[][] jac = new double[_nFP][_nMP]; for (int i = 0, j = 0; i < _nMP; i++) { if (_freeParameters[i]) { jac[j][i] = _transforms[i].transformGradient(functionParameters.get(i)); j++; } } return DoubleMatrix.copyOf(jac); }
protected void assertGradientRoundTrip(final ParameterLimitsTransform transform, final double modelParam) { final double g = transform.transformGradient(modelParam); final double fp = transform.transform(modelParam); final double gInv = transform.inverseTransformGradient(fp); assertEquals(g, 1.0 / gInv, 1e-8); }
protected void assertGradient(final ParameterLimitsTransform transform, final double modelParam) { final double eps = 1e-5; final double g = transform.transformGradient(modelParam); double fdg; try { final double down = transform.transform(modelParam - eps); final double up = transform.transform(modelParam + eps); fdg = (up - down) / 2 / eps; } catch (final IllegalArgumentException e) { final double fp = transform.transform(modelParam); try { final double up = transform.transform(modelParam + eps); fdg = (up - fp) / eps; } catch (final IllegalArgumentException e2) { final double down = transform.transform(modelParam - eps); fdg = (fp - down) / eps; } } assertEquals(g, fdg, 1e-6); }