@Test
public final void testCylinder() {
try {
PhysicalParameters parameters = new PhysicalParameters(25., TemperatureType.C);
double z0 = parameters.calcZ0(BaseRadius);
double waveNumber = parameters.calcWaveNumber(BaseFrequency);
TerminationCalculator term = new IdealOpenEndCalculator();
StateVector sv = term.calcStateVector(null, waveNumber, parameters);
TransferMatrix tm = Tube.calcCylinderMatrix(waveNumber, BaseLength, BaseRadius, parameters);
assertEquals("Determinant incorrect", 1.0, tm.determinant().getReal(), 0.0001);
assertEquals("Determinant.imag incorrect", 0.0, tm.determinant().getImaginary(), 0.0001);
Complex zLoad = tm.multiply(sv).Impedance().divide(z0);
assertEquals("Re(Z) incorrect", 0.03712, zLoad.getReal(), 0.00001);
assertEquals("Im(Z) incorrect", -0.48647, zLoad.getImaginary(), 0.00001);
BoreSectionCalculator boreCalc = new DefaultBoreSectionCalculator();
BoreSection bore = new BoreSection(BaseLength, BaseRadius, BaseRadius);
TransferMatrix tm2 = boreCalc.calcTransferMatrix(bore, waveNumber, parameters);
assertEquals("Determinant 2 incorrect", 1.0, tm2.determinant().getReal(), 0.0001);
assertEquals("Determinant.imag 2 incorrect", 0.0, tm2.determinant().getImaginary(), 0.0001);
Complex zLoad2 = tm2.multiply(sv).Impedance().divide(z0);
assertEquals("Re(Z2) incorrect", 0.03712, zLoad2.getReal(), 0.00001);
assertEquals("Im(Z2) incorrect", -0.48647, zLoad2.getImaginary(), 0.00001);
} catch (Exception e) {
fail(e.getMessage());
}
}
