def test_analytical_rho(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S, K, t, r, sigma = row['S'], row['K'], row['t'], row['R'], row[
'v']
self.assertTrue(
almost_equal(analytical.rho('c', S, K, t, r, sigma),
row['CR'] * .01,
epsilon=.000000001))
self.assertTrue(
almost_equal(analytical.rho('p', S, K, t, r, sigma),
row['PR'] * .01,
epsilon=.000000001))
def test_prices(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S, K, t, r, sigma = row['S'], row['K'], row['t'], row['R'], row[
'v']
self.assertTrue(
almost_equal(black_scholes('c', S, K, t, r, sigma),
row['bs_call'],
epsilon=.000001))
self.assertTrue(
almost_equal(black_scholes('p', S, K, t, r, sigma),
row['bs_put'],
epsilon=.000001))
def test_analytical_rho(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S,K,t,r,sigma = row['S'],row['K'],row['t'],row['R'],row['v']
self.assertTrue(
almost_equal(
analytical.rho('c', S, K, t, r, sigma), row['CR']*.01, epsilon=.000000001
)
)
self.assertTrue(
almost_equal(
analytical.rho('p', S, K, t, r, sigma), row['PR']*.01, epsilon=.000000001
)
)
def test_prices(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S,K,t,r,sigma = row['S'],row['K'],row['t'],row['R'],row['v']
self.assertTrue(
almost_equal(
black_scholes('c', S, K, t, r, sigma), row['bs_call'], epsilon=.000001
)
)
self.assertTrue(
almost_equal(
black_scholes('p', S, K, t, r, sigma), row['bs_put'], epsilon=.000001
)
)
def test_implied_volatility(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S,K,t,r,sigma = row['S'],row['K'],row['t'],row['R'],row['v']
C,P = black_scholes('c', S, K, t, r, sigma), black_scholes('p', S, K, t, r, sigma)
try:
iv = implied_volatility(C, S, K, t, r, 'c')
self.assertTrue(almost_equal(sigma, iv, epsilon = .0001))
except:
print 'could not calculate iv for ', C, S, K, t, r, 'c'
iv = implied_volatility(P, S, K, t, r, 'p')
self.assertTrue(almost_equal(sigma, iv, epsilon = .001) or (iv ==0.0))
def test_implied_volatility(self):
while self.tdi.has_next():
row = self.tdi.next_row()
S, K, t, r, sigma = row['S'], row['K'], row['t'], row['R'], row[
'v']
C, P = black_scholes('c', S, K, t, r,
sigma), black_scholes('p', S, K, t, r, sigma)
try:
iv = implied_volatility(C, S, K, t, r, 'c')
self.assertTrue(almost_equal(sigma, iv, epsilon=.0001))
except:
print('could not calculate iv for ', C, S, K, t, r, 'c')
iv = implied_volatility(P, S, K, t, r, 'p')
self.assertTrue(
almost_equal(sigma, iv, epsilon=.001) or (iv == 0.0))
