/[libvob]/libvob/lava/blob/fshape.py
ViewVC logotype

Diff of /libvob/lava/blob/fshape.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch 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)

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

savannah-hackers-public@gnu.org
ViewVC Help
Powered by ViewVC 1.1.26