X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fmath%2FMatrixI.java;h=2c786532782465107b22ccf181ef48837f8c5e75;hb=187269da0130408d24901c88acb10c217f6284f9;hp=05fe22ffa887fee646501d628d81d1764fbc7614;hpb=d6a30c6516e8e2f923c8c2f6dce7592cb3b7d974;p=jalview.git diff --git a/src/jalview/math/MatrixI.java b/src/jalview/math/MatrixI.java index 05fe22f..2c78653 100644 --- a/src/jalview/math/MatrixI.java +++ b/src/jalview/math/MatrixI.java @@ -1,7 +1,31 @@ +/* + * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) + * Copyright (C) $$Year-Rel$$ The Jalview Authors + * + * This file is part of Jalview. + * + * Jalview is free software: you can redistribute it and/or + * modify it under the terms of the GNU General Public License + * as published by the Free Software Foundation, either version 3 + * of the License, or (at your option) any later version. + * + * Jalview is distributed in the hope that it will be useful, but + * WITHOUT ANY WARRANTY; without even the implied warranty + * of MERCHANTABILITY or FITNESS FOR A PARTICULAR + * PURPOSE. See the GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with Jalview. If not, see . + * The Jalview Authors are detailed in the 'AUTHORS' file. + */ package jalview.math; import java.io.PrintStream; +/** + * An interface that describes a rectangular matrix of double values and + * operations on it + */ public interface MatrixI { /** @@ -42,19 +66,60 @@ public interface MatrixI * @return */ double[] getRow(int i); - + + /** + * Answers a new matrix with a copy of the values in this one + * + * @return + */ MatrixI copy(); + /** + * Returns a new matrix which is the transpose of this one + * + * @return + */ MatrixI transpose(); + /** + * Returns a new matrix which is the result of premultiplying this matrix by + * the supplied argument. If this of size AxB (A rows and B columns), and the + * argument is CxA (C rows and A columns), the result is of size CxB. + * + * @param in + * + * @return + * @throws IllegalArgumentException + * if the number of columns in the pre-multiplier is not equal to + * the number of rows in the multiplicand (this) + */ MatrixI preMultiply(MatrixI m); + /** + * Returns a new matrix which is the result of postmultiplying this matrix by + * the supplied argument. If this of size AxB (A rows and B columns), and the + * argument is BxC (B rows and C columns), the result is of size AxC. + *

+ * This method simply returns the result of in.preMultiply(this) + * + * @param in + * + * @return + * @throws IllegalArgumentException + * if the number of rows in the post-multiplier is not equal to the + * number of columns in the multiplicand (this) + * @see #preMultiply(Matrix) + */ MatrixI postMultiply(MatrixI m); double[] getD(); double[] getE(); + void setD(double[] v); + + void setE(double[] v); + void print(PrintStream ps, String format); void printD(PrintStream ps, String format); @@ -65,12 +130,43 @@ public interface MatrixI void tred(); - double getMaxValue(); + /** + * Reverses the range of the matrix values, so that the smallest values become + * the largest, and the largest become the smallest. This operation supports + * using a distance measure as a similarity measure, or vice versa. + *

+ * If parameter maxToZero is true, then the maximum value becomes + * zero, i.e. all values are subtracted from the maximum. This is consistent + * with converting an identity similarity score to a distance score - the most + * similar (identity) corresponds to zero distance. However note that the + * operation is not reversible (unless the original minimum value is zero). + * For example a range of 10-40 would become 30-0, which would reverse a + * second time to 0-30. Also note that a general similarity measure (such as + * BLOSUM) may give different 'identity' scores for different sequences, so + * they cannot all convert to zero distance. + *

+ * If parameter maxToZero is false, then the values are reflected + * about the average of {min, max} (effectively swapping min and max). This + * operation is reversible. + * + * @param maxToZero + */ + void reverseRange(boolean maxToZero); /** - * Update each value in the matrix by subtracting it from the given value + * Multiply all entries by the given value * - * @param val + * @param d + */ + void multiply(double d); + + /** + * Answers true if the two matrices have the same dimensions, and corresponding values all differ by no + * more than delta (which should be a positive value), else false + * + * @param m2 + * @param delta + * @return */ - void subtractAllFrom(double val); + boolean equals(MatrixI m2, double delta); }