diff --git a/.coverage b/.coverage new file mode 100644 index 0000000000000000000000000000000000000000..a9320ef354670acb84953dfaf7fc2b2f67903333 Binary files /dev/null and b/.coverage differ diff --git a/coverage.xml b/coverage.xml new file mode 100644 index 0000000000000000000000000000000000000000..43946851ee4b397af5acbec809d92a6001249ec3 --- /dev/null +++ b/coverage.xml @@ -0,0 +1,640 @@ +<?xml version="1.0" ?> +<coverage version="7.6.11" timestamp="1739271601600" lines-valid="591" lines-covered="560" line-rate="0.9475" branches-covered="0" branches-valid="0" branch-rate="0" complexity="0"> + <!-- Generated by coverage.py: https://coverage.readthedocs.io/en/7.6.11 --> + <!-- Based on https://raw.githubusercontent.com/cobertura/web/master/htdocs/xml/coverage-04.dtd --> + <sources> + <source>/Users/janhabscheid/Documents/git/fxdgm/fxdgm</source> + </sources> + <packages> + <package name="." line-rate="0.9475" branch-rate="0" complexity="0"> + <classes> + <class name="ElectrolyticDiode.py" filename="ElectrolyticDiode.py" complexity="0" line-rate="0.9384" branch-rate="0"> + <methods/> + <lines> + <line number="10" hits="1"/> + <line number="11" hits="1"/> + <line number="12" hits="1"/> + <line number="13" hits="1"/> + <line number="14" hits="1"/> + <line number="15" hits="1"/> + <line number="16" hits="1"/> + <line number="17" hits="1"/> + <line number="18" hits="1"/> + <line number="19" hits="1"/> + <line number="21" hits="1"/> + <line number="82" hits="1"/> + <line number="83" hits="0"/> + <line number="84" hits="1"/> + <line number="85" hits="1"/> + <line number="86" hits="1"/> + <line number="87" hits="1"/> + <line number="88" hits="0"/> + <line number="89" hits="0"/> + <line number="90" hits="0"/> + <line number="93" hits="1"/> + <line number="94" hits="1"/> + <line number="95" hits="0"/> + <line number="97" hits="1"/> + <line number="98" hits="0"/> + <line number="100" hits="1"/> + <line number="101" hits="1"/> + <line number="103" hits="1"/> + <line number="104" hits="1"/> + <line number="106" hits="1"/> + <line number="107" hits="1"/> + <line number="108" hits="1"/> + <line number="109" hits="1"/> + <line number="111" hits="1"/> + <line number="112" hits="1"/> + <line number="113" hits="1"/> + <line number="114" hits="1"/> + <line number="116" hits="1"/> + <line number="117" hits="1"/> + <line number="118" hits="1"/> + <line number="119" hits="1"/> + <line number="121" hits="1"/> + <line number="122" hits="1"/> + <line number="123" hits="1"/> + <line number="124" hits="1"/> + <line number="127" hits="1"/> + <line number="130" hits="1"/> + <line number="135" hits="1"/> + <line number="136" hits="1"/> + <line number="139" hits="1"/> + <line number="140" hits="1"/> + <line number="141" hits="1"/> + <line number="144" hits="1"/> + <line number="145" hits="1"/> + <line number="146" hits="1"/> + <line number="147" hits="1"/> + <line number="149" hits="1"/> + <line number="150" hits="1"/> + <line number="151" hits="1"/> + <line number="152" hits="1"/> + <line number="153" hits="1"/> + <line number="154" hits="1"/> + <line number="155" hits="1"/> + <line number="156" hits="1"/> + <line number="157" hits="1"/> + <line number="158" hits="1"/> + <line number="159" hits="1"/> + <line number="160" hits="1"/> + <line number="162" hits="1"/> + <line number="163" hits="1"/> + <line number="164" hits="1"/> + <line number="165" hits="1"/> + <line number="166" hits="1"/> + <line number="167" hits="1"/> + <line number="168" hits="1"/> + <line number="169" hits="1"/> + <line number="170" hits="1"/> + <line number="171" hits="1"/> + <line number="172" hits="1"/> + <line number="173" hits="1"/> + <line number="176" hits="1"/> + <line number="177" hits="1"/> + <line number="178" hits="1"/> + <line number="179" hits="1"/> + <line number="181" hits="1"/> + <line number="182" hits="1"/> + <line number="183" hits="1"/> + <line number="184" hits="1"/> + <line number="186" hits="1"/> + <line number="187" hits="1"/> + <line number="188" hits="1"/> + <line number="189" hits="1"/> + <line number="192" hits="1"/> + <line number="212" hits="1"/> + <line number="246" hits="1"/> + <line number="247" hits="1"/> + <line number="250" hits="1"/> + <line number="251" hits="1"/> + <line number="253" hits="1"/> + <line number="254" hits="1"/> + <line number="262" hits="1"/> + <line number="274" hits="1"/> + <line number="275" hits="1"/> + <line number="276" hits="1"/> + <line number="277" hits="1"/> + <line number="278" hits="1"/> + <line number="279" hits="1"/> + <line number="280" hits="1"/> + <line number="281" hits="1"/> + <line number="282" hits="1"/> + <line number="283" hits="1"/> + <line number="284" hits="1"/> + <line number="285" hits="1"/> + <line number="286" hits="1"/> + <line number="287" hits="1"/> + <line number="295" hits="1"/> + <line number="298" hits="1"/> + <line number="299" hits="1"/> + <line number="300" hits="1"/> + <line number="301" hits="1"/> + <line number="303" hits="1"/> + <line number="304" hits="1"/> + <line number="305" hits="1"/> + <line number="306" hits="1"/> + <line number="309" hits="1"/> + <line number="311" hits="1"/> + <line number="312" hits="1"/> + <line number="313" hits="1"/> + <line number="314" hits="1"/> + <line number="315" hits="1"/> + <line number="316" hits="1"/> + <line number="319" hits="1"/> + <line number="320" hits="1"/> + <line number="321" hits="1"/> + <line number="322" hits="1"/> + <line number="325" hits="1"/> + <line number="326" hits="1"/> + <line number="327" hits="1"/> + <line number="328" hits="1"/> + <line number="330" hits="1"/> + <line number="331" hits="1"/> + <line number="333" hits="1"/> + <line number="334" hits="1"/> + <line number="335" hits="0"/> + <line number="336" hits="0"/> + <line number="338" hits="0"/> + </lines> + </class> + <class name="Eq02.py" filename="Eq02.py" complexity="0" line-rate="0.9604" branch-rate="0"> + <methods/> + <lines> + <line number="8" hits="1"/> + <line number="9" hits="1"/> + <line number="10" hits="1"/> + <line number="11" hits="1"/> + <line number="12" hits="1"/> + <line number="13" hits="1"/> + <line number="14" hits="1"/> + <line number="15" hits="1"/> + <line number="16" hits="1"/> + <line number="19" hits="1"/> + <line number="93" hits="1"/> + <line number="94" hits="0"/> + <line number="96" hits="1"/> + <line number="97" hits="1"/> + <line number="100" hits="1"/> + <line number="101" hits="1"/> + <line number="103" hits="1"/> + <line number="104" hits="1"/> + <line number="107" hits="1"/> + <line number="108" hits="1"/> + <line number="110" hits="1"/> + <line number="113" hits="1"/> + <line number="116" hits="1"/> + <line number="117" hits="1"/> + <line number="120" hits="1"/> + <line number="121" hits="1"/> + <line number="122" hits="1"/> + <line number="125" hits="1"/> + <line number="126" hits="1"/> + <line number="129" hits="1"/> + <line number="130" hits="1"/> + <line number="131" hits="1"/> + <line number="132" hits="1"/> + <line number="133" hits="1"/> + <line number="134" hits="1"/> + <line number="137" hits="1"/> + <line number="138" hits="1"/> + <line number="139" hits="1"/> + <line number="140" hits="1"/> + <line number="141" hits="1"/> + <line number="142" hits="1"/> + <line number="146" hits="1"/> + <line number="147" hits="1"/> + <line number="148" hits="1"/> + <line number="149" hits="1"/> + <line number="151" hits="1"/> + <line number="152" hits="1"/> + <line number="156" hits="1"/> + <line number="158" hits="1"/> + <line number="159" hits="1"/> + <line number="160" hits="1"/> + <line number="161" hits="1"/> + <line number="162" hits="1"/> + <line number="163" hits="1"/> + <line number="164" hits="1"/> + <line number="166" hits="1"/> + <line number="167" hits="1"/> + <line number="168" hits="1"/> + <line number="169" hits="1"/> + <line number="170" hits="1"/> + <line number="171" hits="1"/> + <line number="172" hits="1"/> + <line number="176" hits="1"/> + <line number="178" hits="1"/> + <line number="179" hits="1"/> + <line number="182" hits="1"/> + <line number="183" hits="1"/> + <line number="186" hits="1"/> + <line number="187" hits="1"/> + <line number="189" hits="1"/> + <line number="196" hits="1"/> + <line number="199" hits="1"/> + <line number="202" hits="1"/> + <line number="203" hits="1"/> + <line number="204" hits="1"/> + <line number="205" hits="1"/> + <line number="206" hits="1"/> + <line number="208" hits="0"/> + <line number="209" hits="0"/> + <line number="211" hits="0"/> + <line number="212" hits="1"/> + <line number="213" hits="1"/> + <line number="216" hits="1"/> + <line number="217" hits="1"/> + <line number="218" hits="1"/> + <line number="219" hits="1"/> + <line number="222" hits="1"/> + <line number="225" hits="1"/> + <line number="226" hits="1"/> + <line number="227" hits="1"/> + <line number="228" hits="1"/> + <line number="231" hits="1"/> + <line number="232" hits="1"/> + <line number="234" hits="1"/> + <line number="235" hits="1"/> + <line number="237" hits="1"/> + <line number="238" hits="1"/> + <line number="239" hits="1"/> + <line number="240" hits="1"/> + <line number="241" hits="1"/> + <line number="243" hits="1"/> + </lines> + </class> + <class name="Eq04.py" filename="Eq04.py" complexity="0" line-rate="0.9558" branch-rate="0"> + <methods/> + <lines> + <line number="8" hits="1"/> + <line number="9" hits="1"/> + <line number="10" hits="1"/> + <line number="11" hits="1"/> + <line number="12" hits="1"/> + <line number="13" hits="1"/> + <line number="14" hits="1"/> + <line number="15" hits="1"/> + <line number="16" hits="1"/> + <line number="19" hits="1"/> + <line number="89" hits="1"/> + <line number="90" hits="1"/> + <line number="93" hits="1"/> + <line number="94" hits="1"/> + <line number="96" hits="1"/> + <line number="97" hits="1"/> + <line number="100" hits="1"/> + <line number="101" hits="1"/> + <line number="103" hits="1"/> + <line number="106" hits="1"/> + <line number="109" hits="1"/> + <line number="110" hits="1"/> + <line number="113" hits="1"/> + <line number="114" hits="1"/> + <line number="115" hits="1"/> + <line number="118" hits="1"/> + <line number="119" hits="1"/> + <line number="120" hits="1"/> + <line number="121" hits="1"/> + <line number="124" hits="1"/> + <line number="125" hits="1"/> + <line number="126" hits="1"/> + <line number="127" hits="1"/> + <line number="128" hits="1"/> + <line number="129" hits="1"/> + <line number="130" hits="1"/> + <line number="131" hits="1"/> + <line number="132" hits="1"/> + <line number="133" hits="1"/> + <line number="136" hits="1"/> + <line number="137" hits="1"/> + <line number="138" hits="1"/> + <line number="139" hits="1"/> + <line number="140" hits="1"/> + <line number="141" hits="1"/> + <line number="142" hits="1"/> + <line number="143" hits="1"/> + <line number="144" hits="1"/> + <line number="145" hits="1"/> + <line number="149" hits="1"/> + <line number="150" hits="1"/> + <line number="151" hits="1"/> + <line number="152" hits="1"/> + <line number="154" hits="1"/> + <line number="155" hits="1"/> + <line number="157" hits="1"/> + <line number="158" hits="1"/> + <line number="160" hits="1"/> + <line number="161" hits="1"/> + <line number="164" hits="1"/> + <line number="169" hits="1"/> + <line number="171" hits="1"/> + <line number="172" hits="1"/> + <line number="175" hits="1"/> + <line number="176" hits="1"/> + <line number="179" hits="1"/> + <line number="180" hits="1"/> + <line number="184" hits="1"/> + <line number="194" hits="1"/> + <line number="195" hits="1"/> + <line number="201" hits="1"/> + <line number="202" hits="1"/> + <line number="205" hits="1"/> + <line number="206" hits="1"/> + <line number="209" hits="1"/> + <line number="210" hits="1"/> + <line number="212" hits="1"/> + <line number="213" hits="1"/> + <line number="215" hits="1"/> + <line number="225" hits="1"/> + <line number="228" hits="1"/> + <line number="229" hits="1"/> + <line number="230" hits="1"/> + <line number="231" hits="1"/> + <line number="233" hits="1"/> + <line number="234" hits="1"/> + <line number="235" hits="1"/> + <line number="236" hits="1"/> + <line number="239" hits="1"/> + <line number="242" hits="1"/> + <line number="243" hits="1"/> + <line number="244" hits="1"/> + <line number="245" hits="1"/> + <line number="246" hits="1"/> + <line number="248" hits="0"/> + <line number="249" hits="0"/> + <line number="251" hits="0"/> + <line number="252" hits="1"/> + <line number="253" hits="1"/> + <line number="256" hits="1"/> + <line number="257" hits="1"/> + <line number="258" hits="1"/> + <line number="259" hits="1"/> + <line number="262" hits="1"/> + <line number="265" hits="1"/> + <line number="266" hits="1"/> + <line number="267" hits="1"/> + <line number="268" hits="1"/> + <line number="269" hits="1"/> + <line number="270" hits="1"/> + <line number="272" hits="1"/> + <line number="273" hits="0"/> + <line number="274" hits="0"/> + </lines> + </class> + <class name="EqN.py" filename="EqN.py" complexity="0" line-rate="0.9364" branch-rate="0"> + <methods/> + <lines> + <line number="8" hits="1"/> + <line number="9" hits="1"/> + <line number="10" hits="1"/> + <line number="11" hits="1"/> + <line number="12" hits="1"/> + <line number="13" hits="1"/> + <line number="14" hits="1"/> + <line number="15" hits="1"/> + <line number="16" hits="1"/> + <line number="18" hits="1"/> + <line number="82" hits="1"/> + <line number="83" hits="0"/> + <line number="84" hits="1"/> + <line number="85" hits="0"/> + <line number="87" hits="1"/> + <line number="88" hits="1"/> + <line number="91" hits="1"/> + <line number="92" hits="1"/> + <line number="94" hits="1"/> + <line number="95" hits="1"/> + <line number="98" hits="1"/> + <line number="99" hits="1"/> + <line number="101" hits="1"/> + <line number="104" hits="1"/> + <line number="107" hits="1"/> + <line number="108" hits="1"/> + <line number="109" hits="1"/> + <line number="110" hits="1"/> + <line number="113" hits="1"/> + <line number="114" hits="1"/> + <line number="115" hits="1"/> + <line number="117" hits="1"/> + <line number="118" hits="1"/> + <line number="119" hits="1"/> + <line number="120" hits="1"/> + <line number="123" hits="1"/> + <line number="124" hits="1"/> + <line number="127" hits="1"/> + <line number="128" hits="1"/> + <line number="129" hits="1"/> + <line number="130" hits="1"/> + <line number="131" hits="1"/> + <line number="132" hits="1"/> + <line number="135" hits="1"/> + <line number="136" hits="1"/> + <line number="137" hits="1"/> + <line number="138" hits="1"/> + <line number="139" hits="1"/> + <line number="140" hits="1"/> + <line number="143" hits="1"/> + <line number="144" hits="1"/> + <line number="145" hits="1"/> + <line number="146" hits="1"/> + <line number="148" hits="1"/> + <line number="149" hits="1"/> + <line number="152" hits="1"/> + <line number="155" hits="1"/> + <line number="156" hits="1"/> + <line number="157" hits="1"/> + <line number="158" hits="1"/> + <line number="159" hits="1"/> + <line number="160" hits="1"/> + <line number="161" hits="1"/> + <line number="162" hits="1"/> + <line number="165" hits="1"/> + <line number="167" hits="1"/> + <line number="168" hits="1"/> + <line number="169" hits="1"/> + <line number="170" hits="1"/> + <line number="171" hits="1"/> + <line number="174" hits="1"/> + <line number="175" hits="1"/> + <line number="176" hits="1"/> + <line number="177" hits="1"/> + <line number="180" hits="1"/> + <line number="187" hits="1"/> + <line number="188" hits="1"/> + <line number="191" hits="1"/> + <line number="192" hits="0"/> + <line number="194" hits="0"/> + <line number="195" hits="1"/> + <line number="198" hits="1"/> + <line number="199" hits="1"/> + <line number="202" hits="1"/> + <line number="203" hits="1"/> + <line number="204" hits="1"/> + <line number="205" hits="1"/> + <line number="208" hits="1"/> + <line number="211" hits="1"/> + <line number="212" hits="1"/> + <line number="213" hits="1"/> + <line number="214" hits="1"/> + <line number="215" hits="1"/> + <line number="217" hits="0"/> + <line number="218" hits="0"/> + <line number="220" hits="0"/> + <line number="221" hits="1"/> + <line number="222" hits="1"/> + <line number="225" hits="1"/> + <line number="226" hits="1"/> + <line number="227" hits="1"/> + <line number="228" hits="1"/> + <line number="231" hits="1"/> + <line number="232" hits="1"/> + <line number="233" hits="1"/> + <line number="234" hits="1"/> + <line number="235" hits="1"/> + <line number="236" hits="1"/> + <line number="237" hits="1"/> + <line number="238" hits="1"/> + </lines> + </class> + <class name="Helper_DoubleLayerCapacity.py" filename="Helper_DoubleLayerCapacity.py" complexity="0" line-rate="0.977" branch-rate="0"> + <methods/> + <lines> + <line number="4" hits="1"/> + <line number="5" hits="1"/> + <line number="9" hits="1"/> + <line number="23" hits="1"/> + <line number="25" hits="1"/> + <line number="41" hits="1"/> + <line number="43" hits="1"/> + <line number="59" hits="1"/> + <line number="61" hits="1"/> + <line number="77" hits="1"/> + <line number="78" hits="1"/> + <line number="79" hits="1"/> + <line number="80" hits="1"/> + <line number="82" hits="1"/> + <line number="108" hits="1"/> + <line number="109" hits="1"/> + <line number="110" hits="1"/> + <line number="112" hits="1"/> + <line number="144" hits="1"/> + <line number="145" hits="1"/> + <line number="146" hits="1"/> + <line number="147" hits="1"/> + <line number="148" hits="1"/> + <line number="150" hits="1"/> + <line number="188" hits="1"/> + <line number="189" hits="0"/> + <line number="190" hits="1"/> + <line number="191" hits="1"/> + <line number="192" hits="1"/> + <line number="194" hits="1"/> + <line number="195" hits="1"/> + <line number="196" hits="1"/> + <line number="199" hits="1"/> + <line number="200" hits="1"/> + <line number="202" hits="1"/> + <line number="203" hits="1"/> + <line number="204" hits="1"/> + <line number="206" hits="1"/> + <line number="207" hits="1"/> + <line number="208" hits="1"/> + <line number="210" hits="1"/> + <line number="212" hits="1"/> + <line number="214" hits="1"/> + <line number="215" hits="1"/> + <line number="217" hits="1"/> + <line number="218" hits="1"/> + <line number="221" hits="1"/> + <line number="265" hits="1"/> + <line number="266" hits="1"/> + <line number="267" hits="1"/> + <line number="268" hits="1"/> + <line number="269" hits="1"/> + <line number="273" hits="1"/> + <line number="311" hits="1"/> + <line number="312" hits="0"/> + <line number="313" hits="1"/> + <line number="314" hits="1"/> + <line number="315" hits="1"/> + <line number="317" hits="1"/> + <line number="318" hits="1"/> + <line number="319" hits="1"/> + <line number="321" hits="1"/> + <line number="323" hits="1"/> + <line number="324" hits="1"/> + <line number="326" hits="1"/> + <line number="328" hits="1"/> + <line number="329" hits="1"/> + <line number="331" hits="1"/> + <line number="333" hits="1"/> + <line number="334" hits="1"/> + <line number="335" hits="1"/> + <line number="337" hits="1"/> + <line number="338" hits="1"/> + <line number="340" hits="1"/> + <line number="341" hits="1"/> + <line number="343" hits="1"/> + <line number="344" hits="1"/> + <line number="346" hits="1"/> + <line number="348" hits="1"/> + <line number="350" hits="1"/> + <line number="355" hits="1"/> + <line number="403" hits="1"/> + <line number="404" hits="1"/> + <line number="406" hits="1"/> + <line number="407" hits="1"/> + <line number="408" hits="1"/> + <line number="409" hits="1"/> + </lines> + </class> + <class name="RefinedMesh1D.py" filename="RefinedMesh1D.py" complexity="0" line-rate="0.8519" branch-rate="0"> + <methods/> + <lines> + <line number="8" hits="1"/> + <line number="9" hits="1"/> + <line number="10" hits="1"/> + <line number="11" hits="1"/> + <line number="12" hits="1"/> + <line number="16" hits="1"/> + <line number="32" hits="1"/> + <line number="33" hits="1"/> + <line number="34" hits="1"/> + <line number="35" hits="1"/> + <line number="36" hits="1"/> + <line number="37" hits="1"/> + <line number="38" hits="0"/> + <line number="39" hits="0"/> + <line number="40" hits="0"/> + <line number="41" hits="0"/> + <line number="42" hits="1"/> + <line number="43" hits="1"/> + <line number="44" hits="1"/> + <line number="45" hits="1"/> + <line number="46" hits="1"/> + <line number="47" hits="1"/> + <line number="48" hits="1"/> + <line number="49" hits="1"/> + <line number="50" hits="1"/> + <line number="51" hits="1"/> + <line number="52" hits="1"/> + </lines> + </class> + <class name="__init__.py" filename="__init__.py" complexity="0" line-rate="1" branch-rate="0"> + <methods/> + <lines> + <line number="1" hits="1"/> + <line number="3" hits="1"/> + <line number="4" hits="1"/> + <line number="5" hits="1"/> + <line number="6" hits="1"/> + <line number="7" hits="1"/> + <line number="8" hits="1"/> + </lines> + </class> + </classes> + </package> + </packages> +</coverage> diff --git a/fxdgm/ElectrolyticDiode.py b/fxdgm/ElectrolyticDiode.py index 836b48c19f548142f22d74f01577dde207448d40..5016d424344596edacd7af66c3d2ac068a7876a2 100644 --- a/fxdgm/ElectrolyticDiode.py +++ b/fxdgm/ElectrolyticDiode.py @@ -337,7 +337,7 @@ def ElectrolyticDiode(Bias_type:str, phi_bias:float, g_phi:float, z_A:float, z_C else: raise ValueError('Invalid return_type') -if __name__ == '__main__': +if __name__ == '__main__': # pragma: no cover # dont cover main in coverage phi_bias = 10#10 Bias_type = 'ForwardBias' # 'ForwardBias', 'NoBias', 'BackwardBias' g_phi = 350#5 diff --git a/fxdgm/Eq02.py b/fxdgm/Eq02.py index bf40e4165e2fe706d5ec02d6081faacabb06e358..15f96f44e393e36fe193625f4ab8d91b3af0f1e4 100644 --- a/fxdgm/Eq02.py +++ b/fxdgm/Eq02.py @@ -242,47 +242,47 @@ def solve_System_2eq(phi_left:float, phi_right:float, p_right:float, z_A:float, return y_A_vals, y_C_vals, phi_vals, p_vals, x_vals -if __name__ == '__main__': - # Define the parameters - phi_left = 5.0 - phi_right = 0.0 - p_right = 0.0 - y_A_R = 1/3 - y_C_R = 1/3 - z_A = -1.0 - z_C = 1.0 - K = 'incompressible' - Lambda2 = 8.553e-6 - a2 = 7.5412e-4 - number_cells = 1024 - relax_param = .1 - rtol = 1e-4 - max_iter = 500 +# if __name__ == '__main__': +# # Define the parameters +# phi_left = 5.0 +# phi_right = 0.0 +# p_right = 0.0 +# y_A_R = 1/3 +# y_C_R = 1/3 +# z_A = -1.0 +# z_C = 1.0 +# K = 'incompressible' +# Lambda2 = 8.553e-6 +# a2 = 7.5412e-4 +# number_cells = 1024 +# relax_param = .1 +# rtol = 1e-4 +# max_iter = 500 - # Solve the system - y_A, y_C, phi, p, x = solve_System_2eq(phi_left, phi_right, p_right, z_A, z_C, y_A_R, y_C_R, K, Lambda2, a2, number_cells, relax_param=relax_param, x0=0, x1=1, refinement_style='uniform', return_type='Vector', max_iter=max_iter, rtol=rtol) +# # Solve the system +# y_A, y_C, phi, p, x = solve_System_2eq(phi_left, phi_right, p_right, z_A, z_C, y_A_R, y_C_R, K, Lambda2, a2, number_cells, relax_param=relax_param, x0=0, x1=1, refinement_style='uniform', return_type='Vector', max_iter=max_iter, rtol=rtol) - # Plot the solution - plt.plot(x, phi) - plt.xlim(0,0.05) - plt.grid() - plt.xlabel('x [-]') - plt.ylabel('$\\varphi$ [-]') - plt.show() +# # Plot the solution +# plt.plot(x, phi) +# plt.xlim(0,0.05) +# plt.grid() +# plt.xlabel('x [-]') +# plt.ylabel('$\\varphi$ [-]') +# plt.show() - plt.plot(x, y_A, '--', color='tab:blue', label='$y_A$') - plt.plot(x, y_C, '-', color='tab:blue', label='$y_C$') - plt.plot(x, 1 - y_A - y_C, ':', color='tab:blue', label='$y_S$') - plt.xlim(0,0.05) - plt.legend() - plt.grid() - plt.xlabel('x [-]') - plt.ylabel('$y_\\alpha$ [-]') - plt.show() +# plt.plot(x, y_A, '--', color='tab:blue', label='$y_A$') +# plt.plot(x, y_C, '-', color='tab:blue', label='$y_C$') +# plt.plot(x, 1 - y_A - y_C, ':', color='tab:blue', label='$y_S$') +# plt.xlim(0,0.05) +# plt.legend() +# plt.grid() +# plt.xlabel('x [-]') +# plt.ylabel('$y_\\alpha$ [-]') +# plt.show() - plt.plot(x, p) - plt.xlim(0,0.05) - plt.grid() - plt.xlabel('x [-]') - plt.ylabel('$p$ [-]') - plt.show() +# plt.plot(x, p) +# plt.xlim(0,0.05) +# plt.grid() +# plt.xlabel('x [-]') +# plt.ylabel('$p$ [-]') +# plt.show() diff --git a/fxdgm/Eq04.py b/fxdgm/Eq04.py index a370fffc787d8e08fd2e20a613ce08fe14630b94..95cbdf92e67930a15631bc3030e077952a7d9eca 100644 --- a/fxdgm/Eq04.py +++ b/fxdgm/Eq04.py @@ -274,7 +274,7 @@ def solve_System_4eq(phi_left:float, phi_right:float, p_right:float, z_A:float, return y_A, y_C, phi, p, msh -if __name__ == '__main__': +if __name__ == '__main__': # pragma: no cover # dont cover main in coverage # Define the parameters phi_left = 10.0 phi_right = 0.0 diff --git a/fxdgm/EqN.py b/fxdgm/EqN.py index b24a07734fbd85973ca6ec48f979286c5aa725c0..5c847d92dd7439f830c14618eb573ff57f6efb8a 100644 --- a/fxdgm/EqN.py +++ b/fxdgm/EqN.py @@ -238,7 +238,7 @@ def solve_System_Neq(phi_left:float, phi_right:float, p_right:float, z_alpha:lis return y, phi, p, x -if __name__ == '__main__': +if __name__ == '__main__': # pragma: no cover # dont cover main in coverage # Define the parameters phi_left = 8.0 phi_right = 0.0