1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
import static org.junit.jupiter.api.Assertions.assertEquals;
import static org.junit.jupiter.api.Assertions.assertArrayEquals;
import static org.junit.jupiter.api.Assertions.assertThrows;
import java.util.stream.Stream;
import org.junit.jupiter.api.Test;
import org.junit.jupiter.api.DisplayName;
import org.junit.jupiter.params.ParameterizedTest;
import org.junit.jupiter.params.provider.ValueSource;
import org.junit.jupiter.params.provider.MethodSource;
import org.junit.jupiter.params.provider.NullAndEmptySource;
import org.junit.jupiter.params.provider.Arguments;
import static org.junit.jupiter.params.provider.Arguments.arguments;
import com.github.tschoonj.xraylib.Xraylib;
public class TestCrossSections {
@FunctionalInterface
private static interface CrossSectionWrapper {
public double execute(int Z, double energy);
}
static Stream<Arguments> goodValuesProvider() {
return Stream.of(
arguments((CrossSectionWrapper) Xraylib::CS_Photo, 11.451033638148562),
arguments((CrossSectionWrapper) Xraylib::CS_Compt, 0.11785269096475783),
arguments((CrossSectionWrapper) Xraylib::CS_Rayl, 0.39841164641058013),
arguments((CrossSectionWrapper) Xraylib::CS_Total, 11.451033638148562 + 0.11785269096475783 + 0.39841164641058013),
arguments((CrossSectionWrapper) Xraylib::CS_Energy, 11.420221747941419)
);
}
@ParameterizedTest(name="test_good_values {index} -> {0} {1}")
@MethodSource("goodValuesProvider")
public void test_good_values(CrossSectionWrapper wrapper, double expected) {
double cs = wrapper.execute(10, 10.0);
assertEquals(cs, expected, 1E-4);
}
static Stream<Arguments> badValuesProvider() {
return Stream.of(
arguments((CrossSectionWrapper) Xraylib::CS_Photo, 0.09, 1001.0, 0.1, 999.0),
arguments((CrossSectionWrapper) Xraylib::CS_Compt, 0.09, 801.0, 0.1, 800.0),
arguments((CrossSectionWrapper) Xraylib::CS_Rayl, 0.09, 801.0, 0.1, 800.0),
arguments((CrossSectionWrapper) Xraylib::CS_Total, 0.09, 801.0, 0.1, 800.0),
arguments((CrossSectionWrapper) Xraylib::CS_Energy, 0.9, 20001.0, 1.0, 20000.0)
);
}
@ParameterizedTest(name="test_bad_values {index} -> {0} {1} {2}")
@MethodSource("badValuesProvider")
public void test_bad_values(CrossSectionWrapper wrapper, double bad_min, double bad_max, double good_min, double good_max) {
IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(-1, 10.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(Xraylib.ZMAX, 10.0);
});
assertEquals(exc.getMessage(), Xraylib.Z_OUT_OF_RANGE);
exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(26, 0.0);
});
assertEquals(exc.getMessage(), Xraylib.NEGATIVE_ENERGY);
exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(26, -1.0);
});
assertEquals(exc.getMessage(), Xraylib.NEGATIVE_ENERGY);
exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(26, bad_min);
});
assertEquals(exc.getMessage(), Xraylib.SPLINT_X_TOO_LOW);
exc = assertThrows(IllegalArgumentException.class, () -> {
double cs = wrapper.execute(26, bad_max);
});
assertEquals(exc.getMessage(), Xraylib.SPLINT_X_TOO_HIGH);
wrapper.execute(26, good_min);
wrapper.execute(26, good_max);
}
}
|
