/* @Integer */final int lengths[] = { 1, 2 };
/* @Volatility */final double vols[] = { 0.11, 0.50, 1.20 };
final DayCounter dc = new Actual360();
final SimpleQuote spot = new SimpleQuote(0.0);
final SimpleQuote qRate = new SimpleQuote(0.0);
final YieldTermStructure qTS = Utilities.flatRate(qRate, dc);
final SimpleQuote rRate = new SimpleQuote(0.0);
final YieldTermStructure rTS = Utilities.flatRate(rRate, dc);
final SimpleQuote vol = new SimpleQuote(0.0);
final BlackVolTermStructure volTS = Utilities.flatVol(vol, dc);
final BlackScholesMertonProcess stochProcess = new BlackScholesMertonProcess(
new Handle<Quote>(spot),
new Handle<YieldTermStructure>(qTS),
new Handle<YieldTermStructure>(rTS),
new Handle<BlackVolTermStructure>(volTS));
final PricingEngine engine = new AnalyticContinuousGeometricAveragePriceAsianEngine(stochProcess);
final Date today = new Settings().evaluationDate();
for (final Type type : types) {
for (final double strike : strikes) {
for (final int length : lengths) {
final Date exerciseDate = new Date(today.dayOfMonth(), today.month(), today.year() + length);
final EuropeanExercise maturity = new EuropeanExercise(exerciseDate);
final PlainVanillaPayoff payoff = new PlainVanillaPayoff(type, strike);
final ContinuousAveragingAsianOption option = new ContinuousAveragingAsianOption(
AverageType.Geometric, payoff, maturity);
option.setPricingEngine(engine);
/* @Size */final int pastFixings = Integer.MAX_VALUE;
/* @Real */final double runningAverage = Double.NaN;
for (final double u : underlyings) {
for (final double q : qRates) {
for (final double r : rRates) {
for (final double v : vols) {
spot.setValue(u);
qRate.setValue(q);
rRate.setValue(r);
vol.setValue(v);
/* @Real */final double value = option.NPV();
final Map<String, Double> calculated = new HashMap<String, Double>();
calculated.put("delta", option.delta());
calculated.put("gamma", option.gamma());
calculated.put("theta", option.theta());
calculated.put("rho", option.rho());
calculated.put("divRho", option.dividendRho());
calculated.put("vega", option.vega());
final Map<String, Double> expected = new HashMap<String, Double>();
if (value > spot.value() * 1.0e-5) {
// perturb spot and get delta and gamma
/* @Real */final double du = u * 1.0e-4;
spot.setValue(u + du);
/* @Real */double value_p = option.NPV();
/* @Real */final double delta_p = option.delta();
spot.setValue(u - du);
/* @Real */double value_m = option.NPV();
/* @Real */final double delta_m = option.delta();
spot.setValue(u);
expected.put("delta", (value_p - value_m) / (2 * du));
expected.put("gamma", (delta_p - delta_m) / (2 * du));
// perturb rates and get rho and dividend rho
/* @Spread */final double dr = r * 1.0e-4;
rRate.setValue(r + dr);
value_p = option.NPV();
rRate.setValue(r - dr);
value_m = option.NPV();
rRate.setValue(r);
expected.put("rho", (value_p - value_m) / (2 * dr));
/* @Spread */final double dq = q * 1.0e-4;
qRate.setValue(q + dq);
value_p = option.NPV();
qRate.setValue(q - dq);
value_m = option.NPV();
qRate.setValue(q);
expected.put("divRho", (value_p - value_m) / (2 * dq));
// perturb volatility and get vega
/* @Volatility */final double dv = v * 1.0e-4;
vol.setValue(v + dv);
value_p = option.NPV();
vol.setValue(v - dv);
value_m = option.NPV();
vol.setValue(v);
expected.put("vega", (value_p - value_m) / (2 * dv));
// perturb date and get theta
final Date yesterday = today.sub(1);
final Date tomorrow = today.add(1);