Hello Waldek and Ric Thanks so much for the function and comments you sent. The integration works pretty fast for the example function you showed.
For my own use-case, I expect to need very high precision (eg. the Mathematica code I'm trying to port uses 300 digits for the same model). Is DoubleFloat the correct Type to use? Given that the target is 300 digits precision, I tried to run the adaptive integration with lower and lower errors (< 1.0e-16). It ended up taking a much longer time (didn't complete in 3mins). Do any of you have ideas on how to overcome this. Attached are the current versions of the 'beamshapes.input' and the 'a_gaussianquad.input' which holds the integration functions. best, Thejasvi On Wed, Apr 14, 2021 at 12:34 AM Waldek Hebisch <[email protected]> wrote: > On Tue, Apr 13, 2021 at 09:11:55PM +0200, [email protected] wrote: > > @Waldek: I mean, you said that I = \int_a^b f(x)dx = ad_gauss(f, a, b, > 1.0e-12, 25), but looking at the use of h (third parameter) in the code I > would say that I = ad_gauss(f, a, b-a, 1.0e-12, 25). Ric > > Right, thanks for correction. > > -- > Waldek Hebisch > > -- > You received this message because you are subscribed to a topic in the > Google Groups "FriCAS - computer algebra system" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/fricas-devel/WLoTCKr-I3M/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/fricas-devel/20210413223424.GD16648%40math.uni.wroc.pl > . > -- You received this message because you are subscribed to the Google Groups "FriCAS - computer algebra system" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/fricas-devel/CABR_hHdtrthGWemf0-h-%2BoHTOUvCyJCO8e6vt3pEMtwb7sCepA%40mail.gmail.com.
a_gaussianquad.input
Description: Binary data
beamshapes.input
Description: Binary data
