/// <inheritdoc/> /// <exception cref="NumberIsTooLargeException"> if <c>m</c> is an /// <c>OpenMapRealMatrix</c>, and the total number of entries of the product /// is larger than <c>Int32.MaxValue</c>.</exception> public new RealMatrix multiply(RealMatrix m) { try { return(multiply((OpenMapRealMatrix)m)); } catch (InvalidCastException) { MatrixUtils.checkMultiplicationCompatible(this, m); int outCols = m.getColumnDimension(); BlockRealMatrix outp = new BlockRealMatrix(rows, outCols); for (OpenIntToDoubleHashMap.Iterator iterator = entries.iterator(); iterator.MoveNext();) { double value = iterator.Current.Value; int key = iterator.Current.Key; int i = key / columns; int k = key % columns; for (int j = 0; j < outCols; ++j) { outp.addToEntry(i, j, value * m.getEntry(k, j)); } } return(outp); } }
/// <summary> /// Subtract <c>m</c> from this matrix. /// </summary> /// <param name="m">Matrix to be subtracted.</param> /// <returns><c>this</c> - <c>m</c>.</returns> /// <exception cref="MatrixDimensionMismatchException"> if <c>m</c> is not the same /// size as <c>this</c>.</exception> public OpenMapRealMatrix subtract(OpenMapRealMatrix m) { MatrixUtils.checkAdditionCompatible(this, m); OpenMapRealMatrix outp = new OpenMapRealMatrix(this); for (OpenIntToDoubleHashMap.Iterator iterator = m.entries.iterator(); iterator.MoveNext();) { int row = iterator.Current.Key / columns; int col = iterator.Current.Key - row * columns; outp.setEntry(row, col, getEntry(row, col) - iterator.Current.Value); } return(outp); }
/// <summary> /// Postmultiply this matrix by <c>m</c>. /// </summary> /// <param name="m">Matrix to postmultiply by.</param> /// <returns><c>this</c> * <c>m</c>.</returns> /// <exception cref="DimensionMismatchException"> if the number of rows of <c>m</c> /// differ from the number of columns of <c>this</c> matrix.</exception> /// <exception cref="NumberIsTooLargeException"> if the total number of entries of the /// product is larger than <c>Int32.MaxValue</c>.</exception> public OpenMapRealMatrix multiply(OpenMapRealMatrix m) { // Safety check. MatrixUtils.checkMultiplicationCompatible(this, m); int outCols = m.getColumnDimension(); OpenMapRealMatrix outp = new OpenMapRealMatrix(rows, outCols); for (OpenIntToDoubleHashMap.Iterator iterator = entries.iterator(); iterator.MoveNext();) { double value = iterator.Current.Value; int key = iterator.Current.Key; int i = key / columns; int k = key % columns; for (int j = 0; j < outCols; ++j) { int rightKey = m.computeKey(k, j); if (m.entries.containsKey(rightKey)) { int outKey = outp.computeKey(i, j); double outValue = outp.entries.get(outKey) + value * m.entries.get(rightKey); if (outValue == 0.0) { outp.entries.remove(outKey); } else { outp.entries.put(outKey, outValue); } } } } return(outp); }