@@ -590,7 +590,7 @@ def _fit_sandia_field_a(pac, pdc, vdc, pac0, vdc0):
590590 'Parameter fitting requires statsmodels' )
591591 # select data. Avoid very low power, clipping and DC voltage far from
592592 # nominal
593- u = (pac > 0.05 * pac0 ) & (pac < pac0 ) & (np .abs (vdc - vdc0 )/ vdc0 < 0.05 )
593+ u = (pac > 0.05 * pac0 ) & (pac < pac0 ) & (np .abs (vdc - vdc0 )/ vdc0 < 0.05 )
594594 Y = pac [u ]
595595 X = np .array ([pdc [u ]** 2 , pdc [u ]]).T
596596 X = sm .add_constant (X )
@@ -616,10 +616,10 @@ def _f2(params, pac0, vdc0, pdc0, ps0, C0, vdc, pdc, pac):
616616 p ['Vdco' ] = vdc0 # DC V
617617 p ['Pdco' ] = pdc0 # DC power
618618 p ['Pso' ] = ps0 # DC power, small
619- p ['C0' ] = C0 # unitless, tiny
620- p ['C1' ] = params [0 ] # 1/V, tiny
619+ p ['C0' ] = C0 # unitless, tiny
620+ p ['C1' ] = params [0 ] # 1/V, tiny
621621 p ['C2' ] = params [1 ] # 1/V, tiny
622- p ['C3' ] = params [2 ] # 1/V tiny
622+ p ['C3' ] = params [2 ] # 1/V tiny
623623 diff = (_sandia_eff (vdc , pdc , p ) - pac )
624624 return np .sqrt (np .dot (diff , diff ))
625625
@@ -653,7 +653,7 @@ def _fit_sandia_field_b(resid, pac, pdc, vdc, pac0, vdc0, pdc0, ps0, C0):
653653 Parameters for the Sandia inverter model including C1, C2 and C3.
654654 '''
655655 # select data. Avoid very low power and clipping
656- u = (pac > 0.05 * pac0 ) & (pac < pac0 )
656+ u = (pac > 0.05 * pac0 ) & (pac < pac0 )
657657
658658 # initial guess
659659 x0 = np .array ([0. , 0. , 0. ])
0 commit comments