/// <summary>
/// Value a caplet or floorlet under the 1 factor Hull-White model.
/// </summary>
public override void Value(Vector pv, Vector cash, double baseDate, double valueDate, ISACCRResult saccrResult,
IIntraValuationDiagnosticsWriter intraValuationDiagnosticsWriter)
{
int count = fCashflows.Count();
bool forecastIsDiscount = ReferenceEquals(fForecastRate, fDiscountRate);
// time of dfStart and dfEnd
double tDfStart = double.NegativeInfinity;
double tDfEnd = double.NegativeInfinity;
using (var cache = Vector.CacheLike(pv))
{
// Shared between loops
Vector dfStart = cache.Get();
Vector dfEnd = cache.Get();
VectorEngine.For(0, count, LoopDirection.Backwards, i =>
{
using (var innerCache = Vector.CacheLike(pv))
{
CFFloatingInterest cashflow = fCashflows[i];
if (cashflow.Payment_Date < valueDate || cashflow.Payment_Date <= fCutoffDate)
{
return(LoopAction.Break);
}
Vector rate = innerCache.Get();
Vector dfPay = innerCache.Get();
Vector stdDev = innerCache.GetClear();
Vector amount = innerCache.GetClear();
GeneralCashflowProperties properties = fCashflows.GetCashflowProperties(i);
double tPay = CalcUtils.DaysToYears(cashflow.Payment_Date - baseDate);
bool haveDfPay = false;
if (forecastIsDiscount && tPay == tDfStart)
{
dfPay.Assign(dfStart);
haveDfPay = true;
}
using (IntraValuationDiagnosticsHelper.StartCashflow(intraValuationDiagnosticsWriter))
using (var volatilitiesAtDateStore = IntraValuationDiagnosticsHelper.CreateVolatilitiesAtDateStore(intraValuationDiagnosticsWriter, pv.Count))
{
cashflow.AddPropertiesToIntraValuationDiagnostics(intraValuationDiagnosticsWriter);
// Standard Libor implies single reset.
var reset = cashflow.Resets.Single();
if (reset.IsKnown(baseDate))
{
rate.Assign(reset.Known_Rate);
}
else
{
double tValue = CalcUtils.DaysToYears(valueDate - baseDate);
double tReset = CalcUtils.DaysToYears(reset.Reset_Date - baseDate);
double tStart = CalcUtils.DaysToYears(reset.Rate_Start_Date - baseDate);
double tEnd = CalcUtils.DaysToYears(reset.Rate_End_Date - baseDate);
// Reset is a historical or forward Libor rate.
InterestRateUtils.LiborRate(rate, fForecastRate, tValue, tReset, tStart, tEnd, reset.Rate_Year_Fraction,
dfStart, ref tDfStart, dfEnd, ref tDfEnd);
if (tReset > tValue)
{
GetStandardDeviation(stdDev, tValue, tReset, tStart, tEnd);
volatilitiesAtDateStore.Add(valueDate, reset.Reset_Date, stdDev);
}
}
if (!haveDfPay && forecastIsDiscount && tPay == tDfEnd)
{
dfPay.Assign(dfEnd);
haveDfPay = true;
}
// Add swaplet value
amount.AddProduct(properties.Swap_Multiplier, rate);
double tau = reset.Rate_Year_Fraction;
rate.Assign(1.0 + rate * tau);
// Add cap and floor option values.
AddOptionValue(amount, OptionType.Call, rate, properties.Cap_Strike, stdDev, tau, properties.Cap_Multiplier);
AddOptionValue(amount, OptionType.Put, rate, properties.Floor_Strike, stdDev, tau, properties.Floor_Multiplier);
amount.Assign(fBuySellSign * (cashflow.Fixed_Amount + cashflow.Notional * (amount + cashflow.Margin) * cashflow.Accrual_Year_Fraction));
IntraValuationDiagnosticsHelper.AddImpliedVolatilities(intraValuationDiagnosticsWriter, volatilitiesAtDateStore);
CFFixedList.RoundCashflow(amount, Cashflow_Rounding);
CFFixedList.UpdatePvAndCash(cashflow, baseDate, valueDate, haveDfPay ? null : fDiscountRate, null, amount,
dfPay, pv, cash, intraValuationDiagnosticsWriter);
}
}
return(LoopAction.Continue);
});
}
}
/// <summary>
/// Calculate valuation profiles.
/// </summary>
public override void Value(ValuationResults valuationResults, PriceFactorList factors, BaseTimeGrid baseTimes)
{
PreValue(factors);
TimeGridIterator tgi = new TimeGridIterator(fT);
PVProfiles result = valuationResults.Profile;
CashAccumulators cashAccumulators = valuationResults.Cash;
BondOptionDeal deal = (BondOptionDeal)Deal;
double baseDate = factors.BaseDate;
double notional = deal.Notional;
double interval = deal.Coupon_Interval;
double buySign = (deal.Buy_Sell == BuySell.Buy) ? +1 : -1;
double paySign = (deal.Option_Type == OptionType.Call) ? +1 : -1;
double coupon = Percentage.PercentagePoint * deal.Coupon_Rate;
double tExpiry = CalcUtils.DaysToYears(deal.Expiry_Date - baseDate);
double tMaturity = CalcUtils.DaysToYears(deal.Bond_Maturity_Date - baseDate);
IInterestYieldVol interestYieldVol = InterestVolBase.GetYieldVol(factors, deal.Yield_Volatility, fCurrency);
if ((deal.Amortisation) != null && (deal.Amortisation.Count > 0))
{
notional = deal.Amortisation.GetPrincipal(notional, deal.Expiry_Date);
}
bool respectDefault = Respect_Default == YesNo.Yes && fCreditRating != null;
using (IntraValuationDiagnosticsHelper.StartDeal(fIntraValuationDiagnosticsWriter, Deal))
{
using (var pricerCache = Vector.Cache(factors.NumScenarios))
{
Vector defaultTime = null;
Vector bondIsAlive = null;
Vector historicalRecovery = null;
if (respectDefault)
{
defaultTime = pricerCache.Get();
bondIsAlive = pricerCache.Get(1.0);
historicalRecovery = pricerCache.GetClear();
fCreditRating.DefaultTime(defaultTime);
}
var defaultedBeforeBaseDate = respectDefault && CreditRating.DefaultedBeforeBaseDate(fCreditRating, baseDate);
VectorEngine.For(tgi, () =>
{
using (IntraValuationDiagnosticsHelper.StartValuation(fIntraValuationDiagnosticsWriter, tgi.Date))
{
using (var cache = Vector.Cache(factors.NumScenarios))
{
Vector optionValue = cache.GetClear();
Vector stdDev = cache.Get(); // Std.Dev of Price
Vector stdDevYield = cache.Get(); //Std.Dev of Yield
Vector price = cache.Get();
Vector yield = cache.Get();
Vector macaulayDuration = cache.Get();
Vector bondValue = cache.Get();
Vector df = cache.Get();
Vector dfr = fRepoIsDiscount ? null : cache.Get();
if (defaultedBeforeBaseDate)
{
result.AppendVector(tgi.Date, optionValue);
return(LoopAction.Break);
}
// This BondPrice function returns the value of the bond cashflows after ExpiryDate, including accrual, discounted back to T.date
double accrual, cash;
PricingFunctions.BondPrice(bondValue, out accrual, out cash, baseDate, tgi.Date, deal.Expiry_Date, deal.Issue_Date, deal.Bond_Maturity_Date, notional, coupon, fPayDates, fAccruals, fDiscountRate, deal.Amortisation, fPrincipals, fFinalPrincipal, fSurvivalProb, +1.0);
// Now check scenario by scenario for defaults, overwriting bondValue as necessary
if (respectDefault)
{
AdjustBondValueForDefault(notional, tExpiry, bondValue, bondIsAlive, historicalRecovery, defaultTime, tgi.T, fDiscountRate, fRecoveryRate);
}
// convert price and duration to forward (tExpiry) basis
if (tgi.Date == deal.Expiry_Date)
{
optionValue.Assign(buySign * VectorMath.Max(0.0, paySign * (bondValue - notional * fStrike)));
cashAccumulators.Accumulate(fFxRate, tgi.Date, optionValue);
}
else
{
fDiscountRate.GetValue(df, tgi.T, tExpiry);
if (fRepoIsDiscount)
{
dfr = df;
}
else
{
fRepoRate.GetValue(dfr, tgi.T, tExpiry);
}
// Need yield and duration to convert yield vol to price vol.
PricingFunctions.BondForwardPriceAndAdjustedMacaulayDuration(price, macaulayDuration, tgi.T, tExpiry, tMaturity, coupon, interval, df, fDiscountRate, fSurvivalProb);
PricingFunctions.BondYieldFromPrice(yield, tExpiry, tMaturity, coupon, interval, price);
// Calculate Modified Duration from Macaulay Duration.
Vector modifiedDuration = cache.GetClear();
PricingFunctions.GetModifiedDuration(modifiedDuration, macaulayDuration, yield, interval);
// Calculate Std.Dev of Yield and Price
interestYieldVol.GetStdDev(stdDevYield, tgi.T, yield, fStrikeYield, tExpiry, tMaturity - tExpiry);
stdDev.Assign(modifiedDuration * stdDevYield);
if (interestYieldVol.GetDistributionType() == ProbabilityDistribution.Lognormal)
{
stdDev.MultiplyBy(yield);
}
price.AssignQuotient(bondValue, df);
PricingFunctions.BlackFunction(optionValue, deal.Option_Type, price, notional * fStrike, stdDev);
optionValue.MultiplyBy(buySign * dfr);
if (fIntraValuationDiagnosticsWriter.Level > IntraValuationDiagnosticsLevel.None)
{
// Add Intra-valuation Diagnostics
using (var volatilitiesAtDateStore = IntraValuationDiagnosticsHelper.CreateVolatilitiesAtDateStore(fIntraValuationDiagnosticsWriter, factors.NumScenarios))
using (var volatilitiesYieldAtDateStore = IntraValuationDiagnosticsHelper.CreateVolatilitiesAtDateStore(fIntraValuationDiagnosticsWriter, factors.NumScenarios))
{
volatilitiesAtDateStore.Add(tgi.Date, tgi.TimeGrid.fEndDate, stdDev);
volatilitiesYieldAtDateStore.Add(tgi.Date, tgi.TimeGrid.fEndDate, stdDevYield);
IntraValuationDiagnosticsHelper.AddBondOptionProperties(fIntraValuationDiagnosticsWriter, price, dfr, bondValue, accrual,
volatilitiesAtDateStore, volatilitiesYieldAtDateStore);
IntraValuationDiagnosticsHelper.AddCashflowsPV(fIntraValuationDiagnosticsWriter, optionValue);
}
}
}
result.AppendVector(tgi.Date, fFxRate.Get(tgi.T) * optionValue);
return(LoopAction.Continue);
}
}
});
}
result.Complete(fT);
}
}