Package org.jquantlib.math.matrixutilities

Examples of org.jquantlib.math.matrixutilities.Array

148
149
150
151
152
153
154
155
156
157
  }

    @Override
  public Operator subtract(final Operator op) {
        final TridiagonalOperator D = (TridiagonalOperator)op;
        final Array low  = lowerDiagonal.sub(D.lowerDiagonal);
        final Array mid  = diagonal.sub(D.diagonal);
        final Array high = upperDiagonal.sub(D.upperDiagonal);
        return new TridiagonalOperator(low, mid, high);
  }
View Full Code Here
156
157
158
159
160
161
162
163
164
165
        return new TridiagonalOperator(low, mid, high);
  }

    @Override
    public Operator multiply(final double a) {
        final Array low  = lowerDiagonal.mul(a);
        final Array mid  = diagonal.mul(a);
        final Array high = upperDiagonal.mul(a);
        return new TridiagonalOperator(low, mid, high);
    }
View Full Code Here
264
265
266
267
268
269
270
271
272
273
   * Identity instance
   */
    @Override
  public TridiagonalOperator identity(final int size) {
    final TridiagonalOperator I = new TridiagonalOperator(
        new Array(size-1),         // lower diagonal
        new Array(size).fill(1.0), // diagonal
        new Array(size-1));        // upper diagonal
    return I;
  }
View Full Code Here
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
  public Array applyTo(final Array v) /*@ReadOnly*/ {
    if (v.size() != size())
      throw new IllegalStateException("vector of the wrong size (" + v.size() + "instead of " + size() + ")");

    // result = diagonal * v
    final Array result = this.diagonal.mul(v);

    // matricial product
    double d = result.get(0) + (this.upperDiagonal.get(0) * v.get(1));
    result.set(0, d);
    for (int j=1; j<=size()-2; j++) {
        d = result.get(j) + (lowerDiagonal.get(j-1) * v.get(j-1)) + (upperDiagonal.get(j) * v.get(j+1));
      result.set(j, d);
    }
    d = result.get(size()-1) + (lowerDiagonal.get(size()-2) * v.get(size()-2));
    result.set(size()-1, d);

    return result;
  }
View Full Code Here
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
  /**
   * Solve linear system for a given right-hand side
   */
    @Override
  public final Array solveFor(final Array rhs) {
    final Array result = new Array(size());
    final Array tmp = new Array(size());

    double bet = diagonal.first();
    if (bet == 0.0) throw new IllegalStateException("division by zero");

    result.set(0, rhs.first() / bet);
    int j;
    for (j = 1; j <= size() - 1; j++) {
      tmp.set(j, upperDiagonal.get(j - 1) / bet);
      bet = diagonal.get(j) - lowerDiagonal.get(j - 1) * tmp.get(j);
      if (bet == 0.0) throw new IllegalStateException("division by zero");
      result.set(j, (rhs.get(j) - lowerDiagonal.get(j - 1) * result.get(j - 1)) / bet);
    }

    // cannot be j>=0 with Size j
    for (j = size() - 2; j > 0; --j) {
            result.set(j, result.get(j) - (tmp.get(j + 1) * result.get(j + 1)));
        }

    result.set(0, result.first() - (tmp.get(1) * result.get(1)));
    return result;
  }
View Full Code Here
78
79
80
81
82
83
84
85
        inArrears_ = (false);
        zeroPayments_ = (false);
    }

    public CmsLeg withNotionals(/* Real */final double notional) {
        notionals_ = new Array(1).fill(notional);
        return this;
    }
View Full Code Here
108
109
110
111
112
113
114
115
        fixingDays_ = fixingDays; // TODO: clone() ?
        return this;
    }

    public CmsLeg withGearings(/* Real */final double gearing) {
        gearings_ = new Array(1).fill(gearing);
        return this;
    }
View Full Code Here
118
119
120
121
122
123
124
125
        gearings_ = gearings;
        return this;
    }

    public CmsLeg withSpreads(/* Spread */final double spread) {
        spreads_ = new Array(1).fill(spread);
        return this;
    }
View Full Code Here
128
129
130
131
132
133
134
135
        spreads_ = spreads;
        return this;
    }

    public CmsLeg withCaps(/* @Rate */final double cap) {
        caps_ = new Array(1).fill(cap);
        return this;
    }
View Full Code Here
138
139
140
141
142
143
144
145
        caps_ = caps;
        return this;
    }

    public CmsLeg withFloors(/* @Rate */final double floor) {
        floors_ = new Array(1).fill(floor);
        return this;
    }
View Full Code Here
TOP

Related Classes of org.jquantlib.math.matrixutilities.Array

  • org.jquantlib.cashflow.AverageBMALeg
  • org.jquantlib.cashflow.CmsLeg
  • org.jquantlib.cashflow.IborLeg
  • org.jquantlib.instruments.bonds.CmsRateBond
  • org.jquantlib.instruments.bonds.ConvertibleFloatingRateBond
  • org.jquantlib.instruments.bonds.FloatingRateBond
  • org.jquantlib.instruments.ConvertibleFloatingRateBond
  • org.jquantlib.instruments.DiscretizedOption
  • org.jquantlib.instruments.FloatingRateBond
  • org.jquantlib.legacy.libormarkets.LfmCovarianceProxy
    • Copyright © 2018 www.massapicom. All rights reserved.
    All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.