public void testCosine() { float w = 0.02f*2.0f*FLT_PI; int n1 = 1001; float shift = 10.0f; float sigma = 8.0f*shift; float[] f = cos(rampfloat(w*shift,w,n1)); float[] g = cos(rampfloat(0.0f,w,n1)); float[] u = new float[n1]; float[] c = new float[n1]; float[] d = new float[n1]; int min1 = -2*(int)shift; int max1 = 2*(int)shift; LocalShiftFinder lsf = new LocalShiftFinder(shift); lsf.find1(min1,max1,f,g,u,c,d); for (int i1=n1/4; i1<3*n1/4; ++i1) { assertEquals(shift,u[i1],0.02f); assertEquals(1.0f,c[i1],0.02f); assertEquals(1.0f,d[i1],0.02f); } /* edu.mines.jtk.mosaic.SimplePlot.asPoints(f); edu.mines.jtk.mosaic.SimplePlot.asPoints(g); edu.mines.jtk.mosaic.SimplePlot.asPoints(u); edu.mines.jtk.mosaic.SimplePlot.asPoints(c); edu.mines.jtk.mosaic.SimplePlot.asPoints(d); */ } }