Scribd
Upload
129 views2 pages

Numerical Methods Library in FORTRAN: Test Code For

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.

Uploaded by

rodwellhead
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
129 views2 pages

Numerical Methods Library in FORTRAN: Test Code For

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.

Uploaded by

rodwellhead
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 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

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy