Пример #1
0
 def __init__(self, f, xmin, xmax, n):
     self.xmin = xmin
     self.xmax = xmax
     self.dx = (xmax - xmin) / float(n)
     self.invdx = 1 / self.dx
     self.unit = get_unit(f(xmin))
     # Tabulation at midpoints
     x = xmin + (.5 + arange(n)) * self.dx
     try:
         self.f = f(x)
     except:
         # If it fails we try passing the values one by one
         self.f = zeros(n) * f(xmin) # for the unit
         for i in xrange(n):
             self.f[i] = f(x[i])
Пример #2
0
 def __init__(self, f, xmin, xmax, n):
     self.xmin = xmin
     self.xmax = xmax
     self.dx = (xmax - xmin) / float(n)
     self.invdx = 1 / self.dx
     self.unit = get_unit(f(xmin))
     # Tabulation at midpoints
     x = xmin + (.5 + arange(n)) * self.dx
     try:
         self.f = f(x)
     except:
         # If it fails we try passing the values one by one
         self.f = zeros(n) * f(xmin)  # for the unit
         for i in xrange(n):
             self.f[i] = f(x[i])
Пример #3
0
 def __init__(self, f, xmin, xmax, n):
     self.xmin = xmin
     self.xmax = xmax
     self.dx = (xmax - xmin) / float(n - 1)
     self.invdx = 1 / self.dx
     # Not at midpoints here
     x = xmin + arange(n) * self.dx
     self.unit = get_unit(f(xmin))
     try:
         self.f = f(x)
     except:
         # If it fails we try passing the values one by one
         self.f = zeros(n) * f(xmin) # for the unit
         for i in xrange(n):
             self.f[i] = f(x[i])
     self.f = array(self.f)
     self.df = (self.f[range(1, n)] - self.f[range(n - 1)]) * float(self.invdx)
Пример #4
0
 def __init__(self, f, xmin, xmax, n):
     self.xmin = xmin
     self.xmax = xmax
     self.dx = (xmax - xmin) / float(n - 1)
     self.invdx = 1 / self.dx
     # Not at midpoints here
     x = xmin + arange(n) * self.dx
     self.unit = get_unit(f(xmin))
     try:
         self.f = f(x)
     except:
         # If it fails we try passing the values one by one
         self.f = zeros(n) * f(xmin)  # for the unit
         for i in xrange(n):
             self.f[i] = f(x[i])
     self.f = array(self.f)
     self.df = (self.f[range(1, n)] - self.f[range(n - 1)]) * float(
         self.invdx)