def test_impedance_array(self):
high_freq = 10**8 # Hz
low_freq = 0.01 # Hz
decades = 10
R = 100 # ohm
C = 10E-6 # F
f_range = circuits.freq_gen(high_freq, low_freq, decades)
circuit = circuits.cir_RC_parallel(f_range[1], R=R, C=C)
impedance = impedance_array(circuit)
assert isinstance(decades, int), 'the number of decades should\
be an integer'
assert high_freq >= low_freq, 'the low frequency should be smaller\
than the high frequency limit value. Check again.'
np.testing.assert_almost_equal(len(f_range[1]),
impedance[0].shape[0],
decimal=18,
err_msg='the impedance is\
not correclty computed. The number of points in the array is not correct.')
assert len(f_range[1]) == impedance[0].shape[0], 'the frequency and\
the impedance respons do not match in length. Check again.'
assert len(impedance) == 5, 'the impedance array inputted is not the\
right dimensions. The number of columns exceed the expected value (5)'
assert isinstance(impedance[0][1], complex), 'the first column of the\
impedance response should be populated by complex numberes'
assert circuit[0].real == impedance[1][0], 'the complex impedance is\
not separated into its real and imaginary parts correctly.'
assert circuit[0].imag == impedance[2][0], 'the complex impedance is\
def test_rc_parallel(self):
high_freq = 10**8 # Hz
low_freq = 0.01 # Hz
decades = 7
assert isinstance(decades, int),\
'The number of decades should be an integer'
assert high_freq >= low_freq,\
'The low frequency should be smaller than the high\
frequency limit value. Check again.'
f_range = circuits.freq_gen(high_freq, low_freq, decades)
Resistance = 10
Capacitance = 10**-6
assert np.positive(Resistance), 'The input resistance\
is invalid'
assert np.positive(Capacitance), 'The input capacitance\
is invalid'
assert isinstance(Capacitance, float), 'the capacitance should\
be a float, not an integer'
assert Capacitance <= 1, 'the capacitance value is probably too high.'
response = circuits.cir_RC_parallel(f_range[1], Resistance,
Capacitance)
assert len(response) == len(f_range[1]), 'The returned response\
is not valid'
for item in response:
assert isinstance(item, complex), 'The returned response\
def test_RC_parallel(self):
response = circuits.cir_RC_parallel(f_range[1],
R=Resistance,
C=Capacitance)
assert np.positive(Resistance), \
'The input resistance is invalid'
assert np.positive(Capacitance), \
'The input capacitance is invalid'
assert isinstance(Capacitance, float), \
'the capacitance should be a float, not an integer'
assert len(response) == len(f_range[1]), \
'The returned response is not valid'
for item in response:
assert isinstance(item, complex),\
'The returned response includes invalid impedance'
def test_to_dataframe(self):
high_freq = 10**8 # Hz
low_freq = 0.01 # Hz
decades = 10
R = 100 # ohm
C = 10E-6 # F
f_range = circuits.freq_gen(high_freq, low_freq, decades)
circuit = circuits.cir_RC_parallel(f_range[1], R=R, C=C)
impedance_arr = impedance_array(circuit)
dataframe = to_dataframe(f_range, impedance_arr)
assert dataframe.columns[0] == 'freq [Hz]', 'the first column should\
contain the frequency respose.'
assert dataframe.columns[1] == 'angular_freq [1/s]', 'the second\
column should contain the angular frequency respose.'
assert dataframe.columns[3] == 'Re_Z [ohm]', 'the fourth column should\
contain the real impedance part.'
assert dataframe.columns[4] == 'Im_Z [ohm]', 'the fifth column\
should contain the imaginary impedance part.'
assert dataframe.columns[5] == '|Z| [ohm]', 'the sixth column\
should contain the magnitude of the impedance.'
assert dataframe.columns[6] == 'phase_angle [rad]', 'the seventh\
column should contain the phase angle of the impedance.'
assert dataframe['freq [Hz]'][0] == f_range[0][0], 'the first column\
should contain the frequency respose. Check again'
assert dataframe['angular_freq [1/s]'][0] == f_range[1][0],\
'the second column should contain the angular frequency respose.\
Check again'
assert dataframe['Re_Z [ohm]'][0] == circuit[0].real, 'the third\
column should contain the real impedance respose. Check again'
assert dataframe['Im_Z [ohm]'][0] == circuit[0].imag, 'the fourth\
column should contain the real impedance respose. Check again'
assert isinstance(dataframe, pd.DataFrame), 'The output should be\
def test_to_dataframe(self):
high_freq = 10**8 # Hz
low_freq = 0.01 # Hz
decades = 7
R = 100 # ohm
C = 10E-6 # F
n_points = np.round(decades * np.log10(int(high_freq)) -
np.log10(low_freq))
f_range = circuits.freq_gen(high_freq, low_freq, decades=7)
circuit = circuits.cir_RC_parallel(f_range[1], R, C)
impedance_arr = impedance_array(circuit)
dataframe = to_dataframe(f_range, impedance_arr)
assert isinstance(decades, int),\
'the number of decades should be an integer'
assert high_freq >= low_freq,\
'the low frequency should be smaller than the high\
frequency limit value. Check again.'
np.testing.assert_almost_equal(len(f_range[1]),
impedance_arr[0].shape[0],
decimal=18,
err_msg='the impedance is not correclty\
computed. The number of points in the array is not correct.')
assert len(f_range[1]) == impedance_arr[0].shape[0],\
'the frequency and the impedance respons do not match in length.\
Check again.'
assert len(impedance_arr) == 5, 'the impedance array inputted is not\
the right dimensions. The number of columns exceed the expected value (5)'
assert isinstance(impedance_arr[0][1], complex), 'the first column of\
the impedance response should be populated by complex numberes'
assert dataframe.columns[0] == 'freq [Hz]', 'the first column should\
contain the frequency respose.'
assert dataframe.columns[1] == 'angular_freq [1/s]', 'the second\
column should contain the angular frequency respose.'
assert dataframe.columns[2] == 'Re_Z [ohm]', 'the third column should\
contain the real impedance part.'
assert dataframe.columns[3] == 'Im_Z [ohm]', 'the fourth column\
should contain the imaginary impedance part.'
assert dataframe.columns[4] == '|Z| [ohm]', 'the fifth column \
should contain the magnitude of the impedance.'
assert dataframe.columns[5] == 'phase_angle [rad]', 'the sixth column\
should contain the phase angle of the impedance.'
assert dataframe['freq [Hz]'][0] == f_range[0][0], 'the first column\
should contain the frequency respose. Check again'
assert dataframe['angular_freq [1/s]'][0] == f_range[1][0],\
'the second column should contain the angular frequency respose.\
Check again'
assert dataframe['Re_Z [ohm]'][0] == circuit[0].real, 'the third\
column should contain the real impedance respose. Check again'
assert dataframe['Im_Z [ohm]'][0] == circuit[0].imag, 'the fourth\
column should contain the real impedance respose. Check again'
assert isinstance(dataframe, pd.DataFrame), 'The output should be\
def test_RC_simuation(self):
high_freq = 10**8 # Hz
low_freq = 0.01 # Hz
decades = 7
R = 100 # ohm
C = 10E-6 # F
n_points = np.round(decades * np.log10(int(high_freq)) -
np.log10(low_freq))
f_range = circuits.freq_gen(high_freq, low_freq, decades=7)
# Define the RC parallel simulation
circuit_configuration_p = 'parallel'
circuit_parallel = circuits.cir_RC_parallel(f_range[1], R, C)
impedance_data_p = impedance_array(circuit_parallel)
impedance_data_p_df = to_dataframe(f_range, impedance_data_p)
# Define the RC series simulation
circuit_configuration_s = 'series'
circuit_series = circuits.cir_RC_series(f_range[1], R, C)
impedance_data_s = impedance_array(circuit_series)
impedance_data_s_df = to_dataframe(f_range, impedance_data_s)
assert isinstance(decades, int), 'the number of decades should be\
an integer'
assert high_freq >= low_freq, 'the low frequency should be smaller\
than the high frequency limit value. Check again.'
np.testing.assert_almost_equal(len(f_range[1]),
impedance_data_p[0].shape[0],
decimal=18,
err_msg='the impedance\
is not correclty computed. The number of points in the array is not correct.')
np.testing.assert_almost_equal(len(f_range[1]),
impedance_data_s[0].shape[0],
decimal=18,
err_msg='the impedance\
is not correclty computed. The number of points in the array is not correct.')
assert isinstance(impedance_data_p_df, pd.DataFrame), \
'The output should be a pandas.DataFrame'
assert isinstance(impedance_data_s_df, pd.DataFrame), \
'The output should be a pandas.DataFrame'
assert len(f_range[1]) == impedance_data_p[0].shape[0], 'the frequency\
and the impedance respons do not match in length. Check again.'
assert isinstance(circuit_configuration_p, str), 'the circuit\
configuration should be a string.'
assert len(impedance_data_p) == 5, 'the impedance array inputted\
is not the right dimensions. The number of columns exceed the\
expected value (5)'
for item in impedance_data_p[0][:]:
assert isinstance(item, complex), 'the first\
column of the impedance response should be populated by complex numberes'
assert circuit_parallel[0].real == impedance_data_p[1][0], 'the\
complex impedance is not separated into its real and imaginary\
parts correctly.'
assert circuit_parallel[0].imag == impedance_data_p[2][0], 'the\
complex impedance is not separated into its real and imaginary parts\
correctly.'
assert isinstance(circuit_configuration_s, str), 'the circuit\
configuration should be a string.'
assert len(f_range[1]) == impedance_data_s[0].shape[0], 'the frequency\
and the impedance respons do not match in length. Check again.'
assert len(impedance_data_s) == 5, 'the impedance array inputted is\
not the right dimensions. The number of columns exceed the expected value (5)'
for item in impedance_data_s[0][:]:
assert isinstance(item, complex), 'the first\
column of the impedance response should be populated by complex numberes'
assert circuit_series[0].real == impedance_data_s[1][0], 'the complex\
impedance is not separated into its real and imaginary parts correctly.'
assert circuit_series[0].imag == impedance_data_s[2][0], 'the complex\
impedance is not separated into its real and imaginary parts correctly.'
assert isinstance(C, float), 'the capacitance should be a float,\
not an integer'
assert C <= 1, 'the capacitance value is probably too high.'