You are FAST!  That's amazing. Works exactly as advertised.

Re qng:

"The function returns the final approximation, result, an estimate of
the absolute error, abserr and the number of function evaluations
used, neval."

https://www.gnu.org/software/gsl/manual/html_node/QNG-non_002dadaptive-Gauss_002dKronrod-integration.html

This is great.

Bill.

On 25 October 2015 at 22:43, Kurt Pagani <[email protected]> wrote:
>
> As you commanded (files attached) :)
>
> What I do not understand yet: what's the meaning of the other two
> values?  Did you already look into the gsll sources?
>
> *(gsll:integration-qng (lambda (x) (exp (- (* x x)))) 0.0 1.0)
>
> 0.7468241328124271
> 8.291413475940725e-15
> 21
>
>
> kfp@helix:~/Development/GSL$ fricas -nox
> viewman not present, disabling graphics
> Checking for foreign routines
> AXIOM="/usr/local/lib/fricas/target/x86_64-unknown-linux"
> spad-lib="/usr/local/lib/fricas/target/x86_64-unknown-linux/lib/libspad.so"
> foreign routines found
> openServer result 0
>                        FriCAS Computer Algebra System
>                          Version: FriCAS 2014-12-18
>                    Timestamp: Son Okt 25 19:16:29 CET 2015
> -----------------------------------------------------------------------------
>    Issue )copyright to view copyright notices.
>    Issue )summary for a summary of useful system commands.
>    Issue )quit to leave FriCAS and return to shell.
> -----------------------------------------------------------------------------
>
>
> (1) -> )co gsl
>    Compiling FriCAS source code from file
>       /home/kfp/Development/GSL/gsl.spad using old system compiler.
> Value = T
> To load "gsll":
>   Load 1 ASDF system:
>     gsll
> ; Loading "gsll"
> .............
> Value = ("gsll")
> Value = T
>    GSL abbreviates package gsl
> ------------------------------------------------------------------------
>    initializing NRLIB GSL for gsl
>    compiling into NRLIB GSL
>    compiling exported gslIntegrationQng : (DoubleFloat ->
> DoubleFloat,DoubleFloat,DoubleFloat) -> DoubleFloat
> Time: 0.01 SEC.
>
> (time taken in buildFunctor:  0)
>
> ;;;     ***       |gsl| REDEFINED
>
> ;;;     ***       |gsl| REDEFINED
> Time: 0.00 SEC.
>
>
>    Cumulative Statistics for Constructor gsl
>       Time: 0.01 seconds
>
>    finalizing NRLIB GSL
>    Processing gsl for Browser database:
> --------constructor---------
> --------(gslIntegrationQng ((DoubleFloat) (Mapping (DoubleFloat)
> (DoubleFloat)) (DoubleFloat) (DoubleFloat)))---------
> ; compiling file "/home/kfp/Development/GSL/GSL.NRLIB/GSL.lsp" (written
> 26 OCT 2015 03:36:40 AM):
>
> ; file: /home/kfp/Development/GSL/GSL.NRLIB/GSL.lsp
> ; in: SDEFUN |GSL;gslIntegrationQng;M3Df;1|
> ;     (BOOT::SDEFUN BOOT::|GSL;gslIntegrationQng;M3Df;1|
> ;                   ((BOOT::|f| BOOT::|Mapping| (BOOT::|DoubleFloat|)
> ;                     (BOOT::|DoubleFloat|))
> ;                    (BOOT::|a| BOOT::|DoubleFloat|)
> ;                    (BOOT::|b| BOOT::|DoubleFloat|) (BOOT::$
> BOOT::|DoubleFloat|))
> ;                   (GSLL:INTEGRATION-QNG (BOOT::|mkLispFunction1|
> BOOT::|f|)
> ;                                         BOOT::|a| BOOT::|b|))
> ; --> DEFUN PROGN SB-IMPL::%DEFUN SB-IMPL::%DEFUN SB-INT:NAMED-LAMBDA
> ; ==>
> ;   #'(SB-INT:NAMED-LAMBDA BOOT::|GSL;gslIntegrationQng;M3Df;1|
> ;         (BOOT::|f| BOOT::|a| BOOT::|b| BOOT::$)
> ;       (BLOCK BOOT::|GSL;gslIntegrationQng;M3Df;1|
> ;         (GSLL:INTEGRATION-QNG (BOOT::|mkLispFunction1| BOOT::|f|)
> BOOT::|a|
> ;                               BOOT::|b|)))
> ;
> ; caught STYLE-WARNING:
> ;   The variable $ is defined but never used.
> ;
> ; compilation unit finished
> ;   caught 1 STYLE-WARNING condition
>
> ; /home/kfp/Development/GSL/GSL.NRLIB/GSL.fasl written
> ; compilation finished in 0:00:00.054
> ------------------------------------------------------------------------
>    gsl is now explicitly exposed in frame frame1
>    gsl will be automatically loaded when needed from
>       /home/kfp/Development/GSL/GSL.NRLIB/GSL
>
> (1) -> DF ==> DoubleFloat
> f:=(x:DF):DF+->exp(-x^2)
> g(x:DF):DF == exp(-x^2)
> gslIntegrationQng(f,0.0::DF,1.0::DF)
>                                                                    Type:
> Void
> (2) ->
>    (2)  theMap(*1;anonymousFunction;0;frame1;internal)
>                                            Type: (DoubleFloat ->
> DoubleFloat)
> (3) ->    Function declaration g : DoubleFloat -> DoubleFloat has been added
>       to workspace.
>                                                                    Type:
> Void
> (4) ->
>    (4)  0.7468241328124271
>                                                             Type:
> DoubleFloat
> (5) -> gslIntegrationQng(g,0.0::DF,1.0::DF)
>    Compiling function g with type DoubleFloat -> DoubleFloat
>
>    (5)  0.7468241328124271
>                                                             Type:
> DoubleFloat
> (6) ->
>

-- 
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 post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/fricas-devel.
For more options, visit https://groups.google.com/d/optout.

Reply via email to