Java Code Examples for org.apache.commons.math3.analysis.differentiation.DerivativeStructure#add()

public DerivativeStructure value(DerivativeStructure[] variables) {
    DerivativeStructure[] y = new DerivativeStructure[factors.getRowDimension()];
    for (int i = 0; i < y.length; ++i) {
        y[i] = variables[0].getField().getZero();
        for (int j = 0; j < factors.getColumnDimension(); ++j) {
            y[i] = y[i].add(variables[j].multiply(factors.getEntry(i, j)));
        }
    }

    DerivativeStructure sum = variables[0].getField().getZero();
    for (int i = 0; i < y.length; ++i) {
        DerivativeStructure ri = y[i].subtract(target[i]);
        sum = sum.add(ri.multiply(ri));
    }
    return sum;
}
 

public DerivativeStructure value(DerivativeStructure[] variables) {
    DerivativeStructure[] y = new DerivativeStructure[factors.getRowDimension()];
    for (int i = 0; i < y.length; ++i) {
        y[i] = variables[0].getField().getZero();
        for (int j = 0; j < factors.getColumnDimension(); ++j) {
            y[i] = y[i].add(variables[j].multiply(factors.getEntry(i, j)));
        }
    }

    DerivativeStructure sum = variables[0].getField().getZero();
    for (int i = 0; i < y.length; ++i) {
        DerivativeStructure ri = y[i].subtract(target[i]);
        sum = sum.add(ri.multiply(ri));
    }
    return sum;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure x1 = variables[0];
    DerivativeStructure x2 = variables[1];
    DerivativeStructure x3 = variables[2];
    DerivativeStructure tmp1 = variables[0].getField().getZero();
    if (x1.getValue() == 0) {
        tmp1 = tmp1.add((x2.getValue() >= 0) ? 0.25 : -0.25);
    } else {
        tmp1 = x2.divide(x1).atan().divide(twoPi);
        if (x1.getValue() < 0) {
            tmp1 = tmp1.add(0.5);
        }
    }
    DerivativeStructure tmp2 = x1.multiply(x1).add(x2.multiply(x2)).sqrt();
    return new DerivativeStructure[] {
        x3.subtract(tmp1.multiply(10)).multiply(10),
        tmp2.subtract(1).multiply(10),
        x3
    };
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure x1 = variables[0];
    DerivativeStructure x2 = variables[1];
    DerivativeStructure x3 = variables[2];
    DerivativeStructure tmp1 = variables[0].getField().getZero();
    if (x1.getValue() == 0) {
        tmp1 = tmp1.add((x2.getValue() >= 0) ? 0.25 : -0.25);
    } else {
        tmp1 = x2.divide(x1).atan().divide(twoPi);
        if (x1.getValue() < 0) {
            tmp1 = tmp1.add(0.5);
        }
    }
    DerivativeStructure tmp2 = x1.multiply(x1).add(x2.multiply(x2)).sqrt();
    return new DerivativeStructure[] {
        x3.subtract(tmp1.multiply(10)).multiply(10),
        tmp2.subtract(1).multiply(10),
        x3
    };
}
 

/** {@inheritDoc} */
public DerivativeStructure[] value(DerivativeStructure[] point) {

    // extract parameters
    final double[] parameters = new double[point.length];
    for (int k = 0; k < point.length; ++k) {
        parameters[k] = point[k].getValue();
    }

    // compute the residuals
    final DerivativeStructure[] values = new DerivativeStructure[observations.size()];
    int i = 0;
    for (WeightedObservedPoint observed : observations) {

        // build the DerivativeStructure by adding first the value as a constant
        // and then adding derivatives
        DerivativeStructure vi = new DerivativeStructure(point.length, 1, f.value(observed.getX(), parameters));
        for (int k = 0; k < point.length; ++k) {
            vi = vi.add(new DerivativeStructure(point.length, 1, k, 0.0));
        }

        values[i++] = vi;

    }

    return values;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    for (int i = 0; i < (m - 2); ++i) {
        double div = (i + 1) / 29.0;
        DerivativeStructure s1 = variables[0].getField().getZero();
        DerivativeStructure dx = variables[0].getField().getOne();
        for (int j = 1; j < n; ++j) {
            s1 = s1.add(dx.multiply(j).multiply(variables[j]));
            dx = dx.multiply(div);
        }
        DerivativeStructure s2 = variables[0].getField().getZero();
        dx = variables[0].getField().getOne();
        for (int j = 0; j < n; ++j) {
            s2 = s2.add(dx.multiply(variables[j]));
            dx = dx.multiply(div);
        }
        f[i] = s1.subtract(s2.multiply(s2)).subtract(1);
    }

    DerivativeStructure x1 = variables[0];
    DerivativeStructure x2 = variables[1];
    f[m - 2] = x1;
    f[m - 1] = x2.subtract(x1.multiply(x1)).subtract(1);

    return f;

}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
  DerivativeStructure sum = variables[0].getField().getZero();
  for (int i = 0; i < n; ++i) {
    sum = sum.add(variables[i]);
  }
  DerivativeStructure t  = sum.multiply(2.0 / m).add(1);
  DerivativeStructure[] f = new DerivativeStructure[m];
  for (int i = 0; i < n; ++i) {
    f[i] = variables[i].subtract(t);
  }
  Arrays.fill(f, n, m, t.negate());
  return f;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    DerivativeStructure sum = variables[0].getField().getZero();
  for (int i = 1; i < (n - 1); ++i) {
      sum = sum.add(variables[i].multiply(i + 1));
  }
  for (int i = 0; i < (m - 1); ++i) {
    f[i] = sum.multiply(i).subtract(1);
  }
  f[m - 1] = variables[0].getField().getOne().negate();
  return f;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure x1 = variables[0];
    DerivativeStructure x2 = variables[1];
    DerivativeStructure x3 = variables[2];
    DerivativeStructure x4 = variables[3];
  return new DerivativeStructure[] {
    x1.add(x2.multiply(10)),
    x3.subtract(x4).multiply(sqrt5),
    x2.subtract(x3.multiply(2)).multiply(x2.subtract(x3.multiply(2))),
    x1.subtract(x4).multiply(x1.subtract(x4)).multiply(sqrt10)
  };
}
 

public DerivativeStructure getRadius(DerivativeStructure cx, DerivativeStructure cy) {
    DerivativeStructure r = cx.getField().getZero();
    for (Vector2D point : points) {
        r = r.add(distance(point, cx, cy));
    }
    return r.divide(points.size());
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure x1 = variables[0];
    DerivativeStructure x2 = variables[1];
    DerivativeStructure x3 = variables[2];
    DerivativeStructure x4 = variables[3];
  return new DerivativeStructure[] {
    x1.add(x2.multiply(10)),
    x3.subtract(x4).multiply(sqrt5),
    x2.subtract(x3.multiply(2)).multiply(x2.subtract(x3.multiply(2))),
    x1.subtract(x4).multiply(x1.subtract(x4)).multiply(sqrt10)
  };
}
 

/** {@inheritDoc} */
public DerivativeStructure[] value(DerivativeStructure[] point) {

    // extract parameters
    final double[] parameters = new double[point.length];
    for (int k = 0; k < point.length; ++k) {
        parameters[k] = point[k].getValue();
    }

    // compute the residuals
    final DerivativeStructure[] values = new DerivativeStructure[observations.size()];
    int i = 0;
    for (WeightedObservedPoint observed : observations) {

        // build the DerivativeStructure by adding first the value as a constant
        // and then adding derivatives
        DerivativeStructure vi = new DerivativeStructure(point.length, 1, f.value(observed.getX(), parameters));
        for (int k = 0; k < point.length; ++k) {
            vi = vi.add(new DerivativeStructure(point.length, 1, k, 0.0));
        }

        values[i++] = vi;

    }

    return values;
}
 

public DerivativeStructure value(DerivativeStructure[] variables)  {
    DerivativeStructure radius = getRadius(variables[0], variables[1]);

    DerivativeStructure sum = variables[0].getField().getZero();
    for (Vector2D point : points) {
        DerivativeStructure di = distance(point, variables[0], variables[1]).subtract(radius);
        sum = sum.add(di.multiply(di));
    }

    return sum;
}
 

/** {@inheritDoc} */
public DerivativeStructure[] value(DerivativeStructure[] point) {

    // extract parameters
    final double[] parameters = new double[point.length];
    for (int k = 0; k < point.length; ++k) {
        parameters[k] = point[k].getValue();
    }

    // compute the residuals
    final DerivativeStructure[] values = new DerivativeStructure[observations.size()];
    int i = 0;
    for (WeightedObservedPoint observed : observations) {

        // build the DerivativeStructure by adding first the value as a constant
        // and then adding derivatives
        DerivativeStructure vi = new DerivativeStructure(point.length, 1, f.value(observed.getX(), parameters));
        for (int k = 0; k < point.length; ++k) {
            vi = vi.add(new DerivativeStructure(point.length, 1, k, 0.0));
        }

        values[i++] = vi;

    }

    return values;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    DerivativeStructure sum = variables[0].getField().getZero();
  for (int i = 1; i < (n - 1); ++i) {
      sum = sum.add(variables[i].multiply(i + 1));
  }
  for (int i = 0; i < (m - 1); ++i) {
    f[i] = sum.multiply(i).subtract(1);
  }
  f[m - 1] = variables[0].getField().getOne().negate();
  return f;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    DerivativeStructure sum = variables[0].getField().getZero();
    for (int i = 0; i < n; ++i) {
        sum = sum.add(variables[i].multiply(i + 1));
    }
    for (int i = 0; i < m; ++i) {
        f[i] = sum.multiply(i + 1).subtract(1);
    }
    return f;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    DerivativeStructure sum = variables[0].getField().getZero();
    for (int i = 0; i < n; ++i) {
        sum = sum.add(variables[i].multiply(i + 1));
    }
    for (int i = 0; i < m; ++i) {
        f[i] = sum.multiply(i + 1).subtract(1);
    }
    return f;
}
 

public DerivativeStructure getRadius(DerivativeStructure cx, DerivativeStructure cy) {
    DerivativeStructure r = cx.getField().getZero();
    for (Vector2D point : points) {
        r = r.add(distance(point, cx, cy));
    }
    return r.divide(points.size());
}
 

public DerivativeStructure value(DerivativeStructure[] variables)  {
    DerivativeStructure radius = getRadius(variables[0], variables[1]);

    DerivativeStructure sum = variables[0].getField().getZero();
    for (Vector2D point : points) {
        DerivativeStructure di = distance(point, variables[0], variables[1]).subtract(radius);
        sum = sum.add(di.multiply(di));
    }

    return sum;
}
 

@Override
public DerivativeStructure[] value(DerivativeStructure[] variables) {
    DerivativeStructure[] f = new DerivativeStructure[m];
    DerivativeStructure sum = variables[0].getField().getZero();
    for (int i = 0; i < n; ++i) {
        sum = sum.add(variables[i].multiply(i + 1));
    }
    for (int i = 0; i < m; ++i) {
        f[i] = sum.multiply(i + 1).subtract(1);
    }
    return f;
}