@Test
public final void testGradient() {
double[] pos = new double[2];
// dG/dA
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in dG/dA value.",
val[i] / params[2],
g.grad(pos, params, 2),
TOLERANCE);
}
// dG/dx0
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in dG/dx0 value.",
dgdx0_val[i],
g.grad(pos, params, 0),
TOLERANCE);
}
// dG/dy0
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in dG/dy0 value.",
dgdy0_val[i],
g.grad(pos, params, 1),
TOLERANCE);
}
// dG/db
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in dG/db value.",
dgdb_val[i],
g.grad(pos, params, 3),
1e2 * TOLERANCE);
}
}
@Test
public final void testVal() {
double[] pos = new double[2];
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in gaussian value.",
val[i],
g.val(pos, params),
TOLERANCE);
}
}
public static final <T extends RealType<T>> void addGaussianSpotToImage(
RandomAccessibleInterval<T> img, double[] params) {
IterableInterval<T> iterImg = Views.iterable(img);
Cursor<T> lc = iterImg.localizingCursor();
int nDims = img.numDimensions();
double[] position = new double[nDims];
double val;
T var = iterImg.firstElement().createVariable();
while (lc.hasNext()) {
lc.fwd();
lc.localize(position);
val = gaussian.val(position, params);
var.setReal(val);
lc.get().add(var);
}
}
@Test
public final void testHessian() {
double[] pos = new double[2];
// d2G / dA2
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dA² value.",
0,
g.hessian(pos, params, 2, 2),
TOLERANCE);
}
// d2G / dAdx0
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dAdx0 value, ",
d2gdAdx0_val[i],
g.hessian(pos, params, 2, 0),
TOLERANCE);
}
// d2G / dx0dA paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dx0dA value, ",
d2gdAdx0_val[i],
g.hessian(pos, params, 0, 2),
TOLERANCE);
}
// d2G / dAdy0
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dAdy0 value, ",
d2gdAdy0_val[i],
g.hessian(pos, params, 2, 1),
TOLERANCE);
}
// d2G / dy0dA paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dy0dA value, ",
d2gdAdy0_val[i],
g.hessian(pos, params, 1, 2),
TOLERANCE);
}
// d2G / dAdb
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dAdb value, ",
d2gdAdb_val[i],
g.hessian(pos, params, 2, 3),
TOLERANCE);
}
// d2G / dbdA paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dbdA value, ",
d2gdAdb_val[i],
g.hessian(pos, params, 3, 2),
TOLERANCE);
}
// d2G / dx0dy0
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dx0dy0 value, ",
d2gdx0dy0_val[i],
g.hessian(pos, params, 0, 1),
TOLERANCE);
}
// d2G / dy0dx0 paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dy0dx0 value, ",
d2gdx0dy0_val[i],
g.hessian(pos, params, 1, 0),
TOLERANCE);
}
// d2G / dx0db
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dx0db value, ",
d2gdx0db_val[i],
g.hessian(pos, params, 0, 3),
10 * TOLERANCE);
}
// d2G / dbdx0 paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dbdx0 value, ",
d2gdx0db_val[i],
g.hessian(pos, params, 3, 0),
10 * TOLERANCE);
}
// d2G / dy0db
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dy0db value, ",
d2gdy0db_val[i],
g.hessian(pos, params, 1, 3),
10 * TOLERANCE);
}
// d2G / dbdy0 paranoid I tell you
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dbdy0 value, ",
d2gdy0db_val[i],
g.hessian(pos, params, 3, 1),
10 * TOLERANCE);
}
// d2G / db2
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/db² value, ",
d2gdb2_val[i],
g.hessian(pos, params, 3, 3),
100 * TOLERANCE);
}
// d2G / dx02
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dx0² value, ",
d2gdx02_val[i],
g.hessian(pos, params, 0, 0),
TOLERANCE);
}
// d2G / dy02
for (int i = 0; i < X.length; i++) {
pos[0] = X[i];
pos[1] = Y[i];
assertEquals(
"Bad accuracy for x=" + pos[0] + ", y=" + pos[1] + " in d²G/dy0² value, ",
d2gdy02_val[i],
g.hessian(pos, params, 1, 1),
TOLERANCE);
}
}