1+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++++++ ++++++++
++++++++ ++++++++
++++++++ e l l p a c k p r o g r a m ++++++++
++++++++ ++++++++
++++++++ ++++++++
* *********************************************************
* * *
* * example ellpack program 3.a1 *
* * *
* * this is the same equation as the example in *
* * chapter 1. the boundary conditions are changed *
* * and the domain is no longer rectangular. *
* * *
* *********************************************************
option. clockwise = .true.
equation. uxx + uyy + 3*ux - 4*u = exp(x+y)*sin(pi*x)
boundary. ux = 0. on x = 0., y = t for t = -1. to 2.
u = x**2 on x = r, y = 2. for r = 0. to 1.
u = 1. + y/2. on line 1.,2. to 1., 0.
u = x*y/2. on x = 1.-s, y = -s for s = 0. to 1.
grid. 6 x points 0. to 1.
6 y points -1. to 2.
discret. 5 point star
solution. band ge
output. table(u) $ plot(u)
end.
---------------------
preprocessor output
---------------------
2.27 seconds to initialize prepreprocessor
1.30 seconds to parse ellpack program
4.15 seconds to create ellpack control program
7.72 seconds total
ellpack version --- may 1, 1985
1
------------------
domain processor
------------------
boundary points found 19
boundary pieces found 4
grid size 6 by 6
execution sucessful
1
-----------------------
discretization module
-----------------------
5 - p o i n t s t a r
domain non-rectangular
uniform grid 6 x 6
number of equations 19
max no. of unknowns per eq. 5
execution successful
-----------------
solution module
-----------------
b a n d g e
number of equations 19
lower bandwidth 10
upper bandwidth 10
required workspace 608
execution successful
1
----------------
ellpack output
----------------
++++++++++++++++++++++++++++++++++++++++++++++
+ +
+ table of u on 6 x 6 grid +
+ +
++++++++++++++++++++++++++++++++++++++++++++++
x-abscissae are
---------------
0. e+00 2.000000e-01 4.000000e-01 6.000000e-01
8.000000e-01 1.000000e+00
y = 2.000000e+00
----------------
0. e+00 4.000001e-02 1.600000e-01 3.600000e-01
6.400000e-01 2.000000e+00
y = 1.400000e+00
----------------
1.379494e-01 1.847528e-01 3.251630e-01 6.239268e-01
1.098070e+00 1.700000e+00
y = 8.000001e-01
----------------
2.546195e-01 2.956870e-01 4.188895e-01 6.509129e-01
9.915969e-01 1.400000e+00
y = 2.000000e-01
----------------
1.951296e-01 2.248628e-01 3.140624e-01 4.793794e-01
7.290915e-01 1.100000e+00
y =-4.000000e-01
----------------
-5.476528e-02 -6.280090e-02 -8.690773e-02 -1.200000e-01
0. e+00 0. e+00
y =-1.000000e+00
----------------
-5.607122e-02 0. e+00 0. e+00 0. e+00
0. e+00 0. e+00
----------------
ellpack output
----------------
contour plot of u
grid 20 by 20
execution successful
text data bss dec hex
129024 28672 37372 195068 2f9fc
Compiler Messages:
fort.f:
MAIN:
q1pcoe:
r1prhs:
q1bcoe:
r1brhs:
q1bdry:
r0solv:
Wed Aug 14 11:58:43 EST 1985
Contour plot for Problem 2
Solution plot for Problem 2
John Rice's Home Page
ELLPACK Home Page
Purdue University Department of Computer Sciences