Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 93 views

    Numerical Methods Library in FORTRAN: Test Code For

    Uploaded by

    rodwellhead
    User guide to a set of test problems for the fortran routine GLT7 that sets the seven point Gaussian quadrature rule on a standard triangle.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Numerical Methods Library in FORTRAN: Test Code For

    Uploaded by

    rodwellhead
    93 views2 pages
    User guide to a set of test problems for the fortran routine GLT7 that sets the seven point Gaussian quadrature rule on a standard triangle.

    Original Title

    GLT7TESTS.FOR

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    User guide to a set of test problems for the fortran routine GLT7 that sets the seven point Gaussian quadrature rule on a standard triangle.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    93 views2 pages

    Numerical Methods Library in FORTRAN: Test Code For

    Uploaded by

    rodwellhead
    User guide to a set of test problems for the fortran routine GLT7 that sets the seven point Gaussian quadrature rule on a standard triangle.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 2

    www.numerical-methods.

    com

    Numerical Methods Library in FORTRAN


    Module

    SUBROUTINE GLT7TESTS

    Title

    Test code for GLT7.FOR

    Description

    Tests the weights and points , of the Gaussian quadrature rule:


    1

    (, ) ( , ),

    1
    0

    =1

    Which is set by the routine GLT7.FOR.


    The test problems are functions of the form (, ) = on the
    domain of the standard triangle.

    Let
    1

    =
    0

    The values of =

    (1)

    are listed in the following table1.

    K - Reciprocal of I-values
    m

    0
    1
    2
    n
    3
    4
    5

    12

    20

    30

    42

    24

    60 120

    210

    336

    12

    60 180 420

    840 1512

    20

    120 420 1120 2520 5040

    30

    210 840 2520 6300 13860

    42

    336 1512 5040 13860 33264

    Polynomial Test Problems for Quadrature Rules on the Standard Triangle

    www.numerical-methods.com
    Hence
    1

    = 1 .

    Interface

    PROGRAM GLT7TESTS

    Web source of
    code.

    www.numerical-methods.com/fortran/GLT7TESTS.FOR

    Web source of
    this guide

    www.numerical-methods.com/fortan/GLT7TESTS_FOR.htm

    Source of the
    algorithm

    Polynomial Test Problems for Quadrature Rules on the Standard Triangle


    The Laursen and Gellert paper listed in the bibliography states that the
    order precision of the seven point rule is 5. Some experimentation with
    the method on polynomial integrands demosntrates that the seven point
    rule is exact as long as the indices sum to at most five; + 5 .
    In the test problem all non-negative indices in the integrand that sum to
    at most 5 are tested. If the relative error is more that 10-14 an error is
    reported.

    Dependent
    routines

    GLT7.FOR

    Test problems
    or code tested

    Subroutine GLT7 in file GLT7.FOR

    Licence

    This is open source; the software may be used and applied within other
    systems as long as its provenance is appropriately acknowledged.
    See the GNU Licence for more information or contact
    webmaster@numerical-methods.com

    Similar codes
    that may be of
    interest
    Applications

    Bibilography

    Three dimensional boundary element codes at www.boundary-elementmethod.com


    Some criteria for numerically integrated matrices and quadrature
    formulas for triangles" by M.E.Laursen and M.Gellert in International
    Journal for numerical methods in engineering Vol 12 67-76 1978.
    Polynomial Test Problems for Quadrature Rules on the Standard Triangle
    Gaussian Quadrature
    Tutorials on Fortran77

    You might also like