Parent Directory | Revision Log | Patch
revision 1.1 by tjl, Fri May 23 12:30:50 2003 UTC | revision 1.2 by jvk, Fri May 23 14:07:51 2003 UTC | |
---|---|---|
# | Line 28 import os | Line 28 import os |
28 | import sys | import sys |
29 | import RandomArray | import RandomArray |
30 | ||
31 | length= 3 | length= 100 |
32 | epsilon = 0.0001 | epsilon = 0.0001 |
33 | ||
34 | # print l | # print l |
# | Line 41 def L(xleft, xright, rleft, rright): | Line 41 def L(xleft, xright, rleft, rright): |
41 | ||
42 | def U(xleft, xright, rleft, rright): | def U(xleft, xright, rleft, rright): |
43 | ||
44 | return A(xleft, xright, rleft, rright) + L(xleft, xright, rleft, rright) | return A(xleft, xright, rleft, rright) +3* L(xleft, xright, rleft, rright) |
45 | ||
46 | ||
47 | def pot(x, r): | def pot(x, r): |
# | Line 70 def pot(x, r): | Line 70 def pot(x, r): |
70 | dudr[:-1] += dudrl | dudr[:-1] += dudrl |
71 | ||
72 | dudx = r * 0 | dudx = r * 0 |
73 | dudx[1:] += dudxr | #dudx[1:] += dudxr |
74 | dudx[:-1] += dudxl | #dudx[:-1] += dudxl |
75 | ||
76 | return (u, dudx, dudr) | return (u, dudx, dudr) |
77 | ||
# | Line 125 class MomGrad: | Line 125 class MomGrad: |
125 | self.mom *= self.brake | self.mom *= self.brake |
126 | ||
127 | sur = [] | sur = [] |
128 | for i in range(1, 4): | for i in range(1, 2): |
129 | l = 3**i | l = 3**i |
130 | tl = l * length | tl = l * length |
131 | print "TL",tl | print "TL",tl |
132 | x = zeros(2*tl, Float) | x = zeros(2*tl, Float) |
133 | x[0] = -l/2.0 | x[0:tl] = arange(tl) * 1.0 / tl * l |
x[1:tl-1] = arange(tl-2) * (l/(tl-3.0)) - l/2.0 | ||
x[tl-1] = l/2.0 | ||
134 | x[tl] = 0.5 | x[tl] = 0.5 |
135 | x[-1] = 0.5 | x[-1] = 0.5 |
136 | mi = MomGrad(minfunc, x) | mi = MomGrad(minfunc, x) |
|
savannah-hackers-public@gnu.org | ViewVC Help |
Powered by ViewVC 1.1.26 |