コード例 #1
0
 @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);
   }
 }
コード例 #2
0
 @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);
   }
 }
コード例 #3
0
ファイル: LocalizationUtils.java プロジェクト: jacobke/imglib
 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);
   }
 }
コード例 #4
0
 @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);
   }
 }