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 TestPolarized { @FunctionalInterface private static interface CrossSectionZEnergyThetaPhiWrapper { public double execute(int Z, double energy, double theta, double phi); } static Stream<Arguments> goodCrossSectionZEnergyThetaPhiValuesProvider() { return Stream.of( arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Rayl, Math.PI / 4.0, 0.17394690792051704), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Rayl, 0.0, 0.5788126901827545), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Compt, Math.PI / 4.0, 0.005489497545806118) ); } @ParameterizedTest(name="test_cross_sections_z_energy_theta_phi_good_values {index} -> {0} {1}") @MethodSource("goodCrossSectionZEnergyThetaPhiValuesProvider") public void test_cross_sections_z_energy_theta_phi_good_values(CrossSectionZEnergyThetaPhiWrapper wrapper, double theta, double expected) { double cs = wrapper.execute(26, 10.0, theta, Math.PI / 4.0); assertEquals(cs, expected, 1E-4); } static Stream<Arguments> badCrossSectionZEnergyThetaPhiValuesProvider() { return Stream.of( arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Rayl, 0, 10.0, Math.PI / 4.0, Xraylib.Z_OUT_OF_RANGE), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Rayl, Xraylib.ZMAX + 1, 10.0, Math.PI / 4.0, Xraylib.Z_OUT_OF_RANGE), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Rayl, 26, 0.0, Math.PI / 4.0, Xraylib.NEGATIVE_ENERGY), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Compt, 0, 10.0, Math.PI / 4.0, Xraylib.Z_OUT_OF_RANGE), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Compt, Xraylib.ZMAX + 1, 10.0, Math.PI / 4.0, Xraylib.Z_OUT_OF_RANGE), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Compt, 26, 0.0, Math.PI / 4.0, Xraylib.NEGATIVE_ENERGY), arguments((CrossSectionZEnergyThetaPhiWrapper) Xraylib::DCSP_Compt, 26, 10.0, 0.0, Xraylib.NEGATIVE_Q) ); } @ParameterizedTest(name="test_cross_sections_z_energy_theta_phi_bad_values {index} -> {0} {1}") @MethodSource("badCrossSectionZEnergyThetaPhiValuesProvider") public void test_cross_sections_z_energy_theta_phi_bad_values(CrossSectionZEnergyThetaPhiWrapper wrapper, int Z, double energy, double theta, String message) { IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> { double cs = wrapper.execute(Z, energy, theta, Math.PI / 4.0); }); assertEquals(exc.getMessage(), message); } @Test public void test_dcsp_kn_good() { double cs = Xraylib.DCSP_KN(10.0, Math.PI / 4, Math.PI / 4); assertEquals(cs, 0.05888029282784654, 1E-6); } @Test public void test_dcsp_kn_bad() { IllegalArgumentException exc = assertThrows(IllegalArgumentException.class, () -> { double cs = Xraylib.DCSP_KN(0.0, Math.PI / 4, Math.PI / 4); }); assertEquals(exc.getMessage(), Xraylib.NEGATIVE_ENERGY); } @Test public void test_dcsp_thoms() { double cs = Xraylib.DCSP_Thoms(Math.PI / 4, Math.PI / 4); assertEquals(cs, 0.05955590775, 1E-6); } }