1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
import unittest
from pyisda.cdsone import upfront_charge
from pyisda.curve import SpreadCurve
from pyisda.legs import ContingentLeg, FeeLeg
from pyisda.date import roll_date
import datetime
import numpy as np
from quantlib.settings import Settings
from quantlib.time.api import Date
import sys
sys.path.append("..")
from analytics import CreditIndex
from yieldcurve import YC, ql_to_jp, get_curve
class TestUpfront(unittest.TestCase):
index = CreditIndex("ig", 26, "5yr", value_date=datetime.date(2016, 9, 21))
index.notional = 50e6
index.spread = 70
def test_upfront(self):
self.assertAlmostEqual(-self.index.pv, 685292.81, 2)
def test_cdsone(self):
jp_yc = get_curve(self.index.value_date)
fee_dirty = self.index.notional * (
upfront_charge(
self.index.value_date,
self.index._cash_settle_date,
self.index.start_date,
self.index._step_in_date,
self.index.start_date,
self.index.end_date,
self.index.fixed_rate * 1e-4,
jp_yc,
70e-4,
self.index.recovery,
False,
)
)
fee_clean = self.index.notional * (
upfront_charge(
self.index.value_date,
self.index._cash_settle_date,
self.index.start_date,
self.index._step_in_date,
self.index.start_date,
self.index.end_date,
self.index.fixed_rate * 1e-4,
jp_yc,
70e-4,
self.index.recovery,
True,
)
)
self.assertAlmostEqual(-fee_dirty, 685292.81, 2)
self.assertAlmostEqual(fee_clean, self.index.clean_pv)
def test_annuity(self):
self.assertAlmostEqual(
-self.index.clean_pv,
self.index.notional
* self.index.risky_annuity
* (self.index.fixed_rate - self.index.spread)
* 1e-4,
)
class TestSpreadCurve(unittest.TestCase):
def setUp(self):
self.upfront_curve = 1e-2 * np.array(
[
-2.502394,
-4.871879,
-9.329793,
-12.98734,
-15.833254,
-17.622571,
-20.505054,
-24.314297,
]
)
spread_curve = 1e-4 * np.array([500, 500, 500, 500, 500, 500, 500, 500])
recovery_curve = np.full(8, 0.3)
self.trade_date = datetime.date(2018, 6, 18)
self.step_in_date = datetime.date(2018, 6, 19)
self.cash_settle_date = datetime.date(2018, 6, 21)
self.yc = get_curve(self.trade_date, "USD")
self.tenors = np.array((0.5, 1, 2, 3, 4, 5, 7, 10))
self.sc = SpreadCurve(
self.trade_date,
self.yc,
None,
None,
None,
self.tenors,
spread_curve,
self.upfront_curve,
recovery_curve,
ticker="AES",
)
def test_upfront_curves(self):
maturities = [roll_date(self.trade_date, t) for t in self.tenors]
for m, upf in zip(maturities, self.upfront_curve):
pl = ContingentLeg(datetime.date(2017, 9, 20), m, 1)
cl = FeeLeg(datetime.date(2017, 9, 20), m, True, 1, 0.05)
a = pl.pv(
self.trade_date,
self.step_in_date,
self.cash_settle_date,
self.yc,
self.sc,
0.3,
)
b = cl.pv(
self.trade_date,
self.step_in_date,
self.cash_settle_date,
self.yc,
self.sc,
True,
)
self.assertAlmostEqual(a - b, upf)
def test_roundtrip(self):
sc_copy = SpreadCurve.from_bytes(self.sc.as_buffer(True), True)
self.assertEqual(sc_copy.inspect(), self.sc.inspect())
self.assertEqual(sc_copy.full_ticker, self.sc.full_ticker)
sc_copy = SpreadCurve.from_bytes(self.sc.as_buffer(False), False)
self.assertEqual(sc_copy.inspect(), self.sc.inspect())
self.assertEqual(sc_copy.full_ticker, self.sc.full_ticker)
if __name__ == "__main__":
unittest.main()
|
