search

sigma-py / quadpy

:triangular_ruler: Numerical integration (quadrature, cubature) in Python
757 stars 82 forks source link

Return points and values for integrate_adaptive #467

Closed zezhong-zhang closed 1 year ago

zezhong-zhang commented 1 year ago

Just wondering if we could return the points and values evaluated for c1.integrate_adaptive? This would be really helpful for many applications in scientific computing.

If there is no plan for such feature in near future, is it possible to suggest best practice implementation strategy so that I can try myself. I noticed that in the code, it uses a while np.any(~is_good) loop to check the interval adaptively and sum the good results. But this also throws away the points and values evaluated in each loop.

zezhong-zhang commented 1 year ago

Temporarily solved by append a list for the points and values inside of the function to be evaluated by 'integrate_adaptive'. Close the issue for now.