def test_qawf(self):
def f2(x, y):
return Numeric.sin(x) / x
sys2 = integrate.gsl_function(f2, None)
def f1(x, y):
return 1 / x
sys1 = integrate.gsl_function(f1, None)
def f3(x, y):
return 1 / -x
sys3 = integrate.gsl_function(f3, None)
pts = Numeric.array((-Numeric.pi, 0, Numeric.pi))
flag, result1, error = integrate.qagp(sys2, pts, 1e-8, 1e-8, 100000,
self.w)
table1 = integrate.qawo_table(1, 100, integrate.SINE, 100)
cyclew = integrate.workspace(1000000)
flag, result2, error = integrate.qawf(sys1, Numeric.pi, 1e-8, 100,
self.w, cyclew, table1)
table2 = integrate.qawo_table(-1, 100, integrate.SINE, 100)
flag, result3, error = integrate.qawf(sys3, Numeric.pi, 1e-8, 100,
self.w, cyclew, table2)
assert (Numeric.absolute(result1 + result2 + result3 - Numeric.pi) <
1e-8)
def run():
def f2(x,y):
return numx.sin(x) / x
sys2 = integrate.gsl_function(f2, None)
def f1(x,y):
return 1 / x
sys1 = integrate.gsl_function(f1, None)
def f3(x,y):
return 1 / -x
sys3 = integrate.gsl_function(f3, None)
w = integrate.workspace(1000000)
cyclew = integrate.workspace(1000000)
table1 = integrate.qawo_table(1, 100, integrate.SINE, 100)
table2 = integrate.qawo_table(-1, 100, integrate.SINE, 100)
# Borders and singualrity for gagp
pts = numx.array((-numx.pi, 0, numx.pi))
flag, result1, error = integrate.qagp(sys2, pts, 1e-8, 1e-8, 100000, w)
flag, result2, error = integrate.qawf(sys1, numx.pi, 1e-8, 100, w,
cyclew, table1)
flag, result3, error = integrate.qawf(sys3, numx.pi, 1e-8, 100, w,
cyclew, table2)
print "Result of integration is :", result1 + result2 + result3
def run():
def f2(x, y):
return numx.sin(x) / x
sys2 = integrate.gsl_function(f2, None)
def f1(x, y):
return 1 / x
sys1 = integrate.gsl_function(f1, None)
def f3(x, y):
return 1 / -x
sys3 = integrate.gsl_function(f3, None)
w = integrate.workspace(1000000)
cyclew = integrate.workspace(1000000)
table1 = integrate.qawo_table(1, 100, integrate.SINE, 100)
table2 = integrate.qawo_table(-1, 100, integrate.SINE, 100)
# Borders and singualrity for gagp
pts = numx.array((-numx.pi, 0, numx.pi))
flag, result1, error = integrate.qagp(sys2, pts, 1e-8, 1e-8, 100000, w)
flag, result2, error = integrate.qawf(sys1, numx.pi, 1e-8, 100, w, cyclew,
table1)
flag, result3, error = integrate.qawf(sys3, numx.pi, 1e-8, 100, w, cyclew,
table2)
print "Result of integration is :", result1 + result2 + result3
def test_qawo(self):
def f(x,y):
return 1
sys = integrate.gsl_function(f, None)
table = integrate.qawo_table(1, 2*Numeric.pi, integrate.COSINE, 100)
flag, result, error = integrate.qawo(sys, -Numeric.pi, 1e-8, 1e-8, 100, self.w, table)
assert(Numeric.absolute(result) < 1e-7)
assert(Numeric.absolute(error) < 1e-7)
def test_qawf(self):
def f2(x,y):
return Numeric.sin(x) / x
sys2 = integrate.gsl_function(f2, None)
def f1(x,y):
return 1 / x
sys1 = integrate.gsl_function(f1, None)
def f3(x,y):
return 1 / -x
sys3 = integrate.gsl_function(f3, None)
pts = Numeric.array((-Numeric.pi, 0, Numeric.pi))
flag, result1, error = integrate.qagp(sys2, pts, 1e-8, 1e-8, 100000, self.w)
table1 = integrate.qawo_table(1, 100, integrate.SINE, 100)
cyclew = integrate.workspace(1000000)
flag, result2, error = integrate.qawf(sys1, Numeric.pi, 1e-8, 100, self.w, cyclew, table1)
table2 = integrate.qawo_table(-1, 100, integrate.SINE, 100)
flag, result3, error = integrate.qawf(sys3, Numeric.pi, 1e-8, 100, self.w, cyclew, table2)
assert(Numeric.absolute(result1+result2+result3 - Numeric.pi) < 1e-8)
