Use scipy's find_minimum in lambertw for MPP #2567
Conversation
|
What do you think of this approach? |
pvlib/singlediode.py
Outdated
| if res.success.all(): | ||
| v_mp = res.x | ||
| p_mp = -1.*res.f_x | ||
| else: | ||
| use_gs = True |
There was a problem hiding this comment.
I am inclined to not include this safety net; if we plan to get rid of the GS function eventually, might as well let any unforeseen issues with the scipy approach bubble to the surface now so that we can address them sooner rather than later.
Or is there reason to expect the scipy function to fail sometimes?
There was a problem hiding this comment.
I can't think of a reason for failure other than bad inputs. Because res.success is a vector, we could pass over to the GS only those index values which failed, for whatever reason. Or skip the success check altogether and just push v_mp forward. I worry a little that will forward unconverged, but still numeric, values.
|
|
||
| # Find the voltage, v_mp, where the power is maximized. | ||
| # Start the golden section search at v_oc * 1.14 | ||
| p_mp, v_mp = _golden_sect_DataFrame(params, 0., v_oc * 1.14, _pwr_optfcn) |
There was a problem hiding this comment.
Just curious... I wonder if anyone knows why 1.14 is the magic number. I see it goes back to the early early days of pvlib-python.
There was a problem hiding this comment.
That was 1.14 * Voc_ref, the STC value of Voc. The factor of 1.14 was meant to push the right endpoint out when temperature is lower than STC.
In pvlib, v_oc reflects temperature and irradiance. I included 1.01 * v_oc in case the algorithm needs to start with a negative value of power at the right endpoint, since it is minimizing negative power. Now that I write this, I don't think that's necessary - it is provable that power has a unique maximum between v=0 and v=v_oc.
optimize.elementwisefunctions #2497docs/sphinx/source/referencefor API changes.docs/sphinx/source/whatsnewfor all changes. Includes link to the GitHub Issue with:issue:`num`or this Pull Request with:pull:`num`. Includes contributor name and/or GitHub username (link with:ghuser:`user`).remote-data) and Milestone are assigned to the Pull Request and linked Issue.Switches to the old golden section if
find_minimumis not available.Since
find_minimumreplaces the golden mean search, if scipy>=1.15, the current tests seem sufficient.