use of uk.ac.sussex.gdsc.test.junit5.SeededTest in project GDSC-SMLM by aherbert.
the class ConvolutionTest method canComputeDoubleScaledConvolution.
@SeededTest
void canComputeDoubleScaledConvolution(RandomSeed seed) {
final UniformRandomProvider random = RngUtils.create(seed.getSeed());
final TDoubleArrayList list = new TDoubleArrayList();
int size = 10;
for (int i = 0; i < sizeLoops / 2; i++) {
double sd = 0.5;
for (int j = 0; j < sdLoops; j++) {
final double[] data1 = randomData(random, size);
final double[] data2 = randomData(random, size);
final double[] kernel = createKernel(sd);
for (int scale = 2; scale < 5; scale++) {
final double[] e1 = convolve(kernel, data1, list, scale);
final double[] e2 = convolve(kernel, data2, list, scale);
final double[][] o = Convolution.convolve(kernel, data1, data2, scale);
final double[][] o2 = new double[2][o[0].length];
Convolution.convolve(kernel, data1, data2, scale, new DoubleConvolutionValueProcedure() {
int index = 0;
@Override
public boolean execute(double value1, double value2) {
o2[0][index] = value1;
o2[1][index] = value2;
index++;
return true;
}
});
Assertions.assertArrayEquals(e1, o[0]);
Assertions.assertArrayEquals(e1, o2[0]);
Assertions.assertArrayEquals(e2, o[1]);
Assertions.assertArrayEquals(e2, o2[1]);
}
sd *= 2;
}
size *= 2;
}
}
use of uk.ac.sussex.gdsc.test.junit5.SeededTest in project GDSC-SMLM by aherbert.
the class ConvolutionTest method canComputeDoubleConvolution.
@SeededTest
void canComputeDoubleConvolution(RandomSeed seed) {
final UniformRandomProvider random = RngUtils.create(seed.getSeed());
int size = 10;
for (int i = 0; i < sizeLoops; i++) {
double sd = 0.5;
for (int j = 0; j < sdLoops; j++) {
final double[] data1 = randomData(random, size);
final double[] data2 = randomData(random, size);
final double[] kernel = createKernel(sd);
double[] e1;
double[] e2;
double[][] r1;
for (int fft = 0; fft < 2; fft++) {
if (fft == 1) {
e1 = Convolution.convolveFft(kernel, data1);
e2 = Convolution.convolveFft(kernel, data2);
r1 = Convolution.convolveFft(kernel, data1, data2);
} else {
e1 = Convolution.convolve(kernel, data1);
e2 = Convolution.convolve(kernel, data2);
r1 = Convolution.convolve(kernel, data1, data2);
}
Assertions.assertEquals(r1.length, 2);
Assertions.assertEquals(e1.length, r1[0].length);
Assertions.assertEquals(e2.length, r1[1].length);
for (int k = 0; k < e1.length; k++) {
// Exact match
Assertions.assertEquals(e1[k], r1[0][k]);
Assertions.assertEquals(e2[k], r1[1][k]);
}
}
sd *= 2;
}
size *= 2;
}
}
use of uk.ac.sussex.gdsc.test.junit5.SeededTest in project GDSC-SMLM by aherbert.
the class ConvolutionTest method canComputeDoubleScaledConvolutionWithEarlyExit.
@SeededTest
void canComputeDoubleScaledConvolutionWithEarlyExit(RandomSeed seed) {
final UniformRandomProvider random = RngUtils.create(seed.getSeed());
int size = 10;
final int sizeLoops = 4;
final int sLoops = 2;
for (int i = 0; i < sizeLoops; i++) {
double sd = 0.5;
for (int j = 0; j < sLoops; j++) {
final double[] data1 = randomData(random, size);
final double[] data2 = randomData(random, size);
final double[] kernel = createKernel(sd);
for (int scale = 2; scale < 5; scale++) {
final double[][] e = Convolution.convolve(kernel, data1, data2, scale);
final double[][] o = new double[2][e[0].length];
final int limit = data1.length;
Convolution.convolve(kernel, data1, data2, scale, new DoubleConvolutionValueProcedure() {
int index = 0;
@Override
public boolean execute(double value1, double value2) {
o[0][index] = value1;
o[1][index] = value1;
index++;
return index < limit;
}
});
int index = 0;
for (; index < limit; index++) {
Assertions.assertEquals(e[0][index], o[0][index]);
Assertions.assertEquals(e[0][index], o[1][index]);
}
while (index < o.length) {
Assertions.assertEquals(0, o[0][index]);
Assertions.assertEquals(0, o[1][index]);
index++;
}
}
sd *= 2;
}
size *= 2;
}
}
use of uk.ac.sussex.gdsc.test.junit5.SeededTest in project GDSC-SMLM by aherbert.
the class ConvolutionTest method canComputeScaledConvolutionWithEarlyExit.
@SeededTest
void canComputeScaledConvolutionWithEarlyExit(RandomSeed seed) {
final UniformRandomProvider random = RngUtils.create(seed.getSeed());
int size = 10;
final int sizeLoops = 4;
final int sLoops = 2;
for (int i = 0; i < sizeLoops; i++) {
double sd = 0.5;
for (int j = 0; j < sLoops; j++) {
final double[] data = randomData(random, size);
final double[] kernel = createKernel(sd);
for (int scale = 2; scale < 5; scale++) {
final double[] e = Convolution.convolve(kernel, data, scale);
final double[] o = new double[e.length];
final int limit = data.length;
Convolution.convolve(kernel, data, scale, new ConvolutionValueProcedure() {
int index = 0;
@Override
public boolean execute(double value) {
o[index++] = value;
return index < limit;
}
});
int index = 0;
for (; index < limit; index++) {
Assertions.assertEquals(e[index], o[index]);
}
while (index < o.length) {
Assertions.assertEquals(0, o[index++]);
}
}
sd *= 2;
}
size *= 2;
}
}
use of uk.ac.sussex.gdsc.test.junit5.SeededTest in project GDSC-SMLM by aherbert.
the class ConvolutionTest method doSpeedTest.
@SpeedTag
@SeededTest
void doSpeedTest(RandomSeed seed) {
Assumptions.assumeTrue(logger.isLoggable(Level.INFO));
Assumptions.assumeTrue(TestSettings.allow(TestComplexity.HIGH));
final UniformRandomProvider rg = RngUtils.create(seed.getSeed());
int size = 10;
for (int i = 0; i < sizeLoops; i++) {
double sd = 0.5;
for (int j = 0; j < sdLoops; j++) {
speedTest(rg, size, sd);
sd *= 2;
}
size *= 2;
}
}
Aggregations