-
Notifications
You must be signed in to change notification settings - Fork 244
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
[RomApp] Fix Bug in ECM #11732
[RomApp] Fix Bug in ECM #11732
Conversation
just a question: not enforcing that constraint implies giving back with the capability of exactly integrating the areas, is that correct? |
The intention of the constraint is to avoid the trivial solution of all weights equal to zero. A by-product of the contraint is that the sum of the weights equal the number of elements. If the contraint is not enforced, the sum of the weights can be whatever value the greedy algorithm comes up with, which can be very large or very small according to the specific problem. Imposing this condition should not have a negative effect besides requiring a single element extra, with respect to not imposing it. However, in practice we are observing that enforcing this constraint sometimes prevents the algorithm from converging. We are looking into the reasons why this is happening. In the meantime, this PR provides an easy fix by allowing the user to choose not to impose the constraint if convergence is difficult for the specific case. |
Can you please take a look why this is not passing the CI? |
@@ -17,7 +17,7 @@ class EmpiricalCubatureMethod(): | |||
def __init__( | |||
self, | |||
ECM_tolerance = 0, | |||
Filter_tolerance = 1e-16, | |||
Filter_tolerance = 0, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Why is the tolerance here set to 0 by default?
I mean, do you think it will most likely need to use all elements/conditions to converge?
Everything else LGTM. |
Could we approve this then? |
no objection on my side. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM
📝 Description
This PR is intended to fix two problems in the ECM algorithm:
🆕 Changelog