NEA-0930 FELPO.

FELPO, 2-D Minimization of Quadratic Functionals by Finite Elements Method

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROBLEM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHOR, MATERIAL, CATEGORIES

[ top ]

1. NAME OR DESIGNATION OF PROGRAM: FELPO.

[ top ]

2. COMPUTERS

To submit a request, click below on the link of the version you wish to order. Only liaison officers are authorised to submit online requests. Rules for requesters are available here.

Program name	Package id	Status	Status date
FELPO	NEA-0930/01	Tested	12-JUN-1985

Machines used:

Package ID	Orig. computer	Test computer
NEA-0930/01	DEC VAX 11/780	DEC VAX 11/780

[ top ]

3. DESCRIPTION OF PROBLEM OR FUNCTION

FELPO is a finite element program for the solution of minimization problems of quadratic functionals on two-dimensional domains, as they arise in the numerical solution of elliptical partial differential equations of second order.

[ top ]

4. METHOD OF SOLUTION

Consider a domain G(X,Y) with the boundary DG,  where a function U(X,Y) has to be calculated, which minimzes the quadratic functional over G with the integrand:

F(U) = F1*UXX2 + 2*F2*UX*UY + F3*UYY2 + F4*U2 + 2*F5*U

The Fi are functions of X and Y, which have to satisfy the conditions:

     F1,F3 > 0; F2 >=0   and   F1*F3 - F2*F2 > 0

On the boundary DGE, essential boundary conditions

                   U=U(DGE)

have to be prescribed. On the rest of the boundary DG-DGE, the solution fulfills the natural boundary conditions

   (F1*UX+F2*UY)*NX+(F2*UX+F3*UY)*NY = 0

automatically. NX and NY are the direction-cosines of the outer normal of the boundary.

The solution is equivalent to the solution of the partial differential equation

   DX(F1*UX+F2*UY) + DY(F2*UX+F3*UY) = F4*U + F5

under the boundary conditions described above.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM:

[ top ]

6. TYPICAL RUNNING TIME:

[ top ]

7. UNUSUAL FEATURES OF THE PROGRAM:

[ top ]

8. RELATED AND AUXILIARY PROGRAMS:

[ top ]

9. STATUS

Package ID	Status date	Status
NEA-0930/01	12-JUN-1985	Tested at NEADB

[ top ]

10. REFERENCES:

[ top ]

11. MACHINE REQUIREMENTS

Less than 180K bytes of main storage are required for execution on VAX-11/780.

[ top ]

12. PROGRAMMING LANGUAGE(S) USED

Package ID	Computer language
NEA-0930/01	FORTRAN-IV

[ top ]

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:
VMS V3.4 (VAX-11/780).

[ top ]

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

[ top ]

15. NAME AND ESTABLISHMENT OF AUTHOR

     Jens Poepplau
     Scientific Control Systems GmbH (SCS)
     Oehleckerring 40
     D-2000 HAMBURG 62
     Federal Republic of Germany

[ top ]

16. MATERIAL AVAILABLE

NEA-0930/01

File name	File description	Records
NEA0930_01.003	INFORMATION FILE	51
NEA0930_01.004	FELPO SOURCE PROGRAM	1936
NEA0930_01.005	FKOFF ROUTINE FOR TEST CASE 1	32
NEA0930_01.006	FKOFF ROUTINE FOR TEST CASE 2	31
NEA0930_01.007	FKOFF ROUTINE FOR TEST CASE 3	41
NEA0930_01.008	FKOFF ROUTINE FOR TEST CASE 4	35
NEA0930_01.009	COMMAND FILE TO RUN TEST CASES	43
NEA0930_01.010	FELPO TEST CASE 1.1 INPUT DATA	48
NEA0930_01.011	FELPO TEST CASE 1.1 PRINTED OUTPUT	150
NEA0930_01.012	FELPO TEST CASE 1.2 INPUT DATA	48
NEA0930_01.013	FELPO TEST CASE 1.2 PRINTED OUTPUT	150
NEA0930_01.014	FELPO TEST CASE 1.3 INPUT DATA	27
NEA0930_01.015	FELPO TEST CASE 1.3 PRINTED OUTPUT	105
NEA0930_01.016	FELPO TEST CASE 2.1 INPUT DATA	87
NEA0930_01.017	FELPO TEST CASE 2.1 PRINTED OUTPUT	218
NEA0930_01.018	FELPO TEST CASE 3.1 INPUT DATA	294
NEA0930_01.019	FELPO TEST CASE 3.1 PRINTED OUTPUT	644
NEA0930_01.020	FELPO TEST CASE 3.2 INPUT DATA	294
NEA0930_01.021	FELPO TEST CASE 3.2 PRINTED OUTPUT	614
NEA0930_01.022	FELPO TEST CASE 4.1 INPUT DATA	79
NEA0930_01.023	FELPO TEST CASE 4.1 PRINTED OUTPUT	210

[ top ]

17. CATEGORIES

P. General Mathematical and Computing System Routines

Keywords: finite element method, partial differential equations, two-dimensional.