X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Futil%2FComparison.java;h=08d2c88088f9b7101c3577d677dc6d7fbc4c9975;hb=a8f483d04205bb8273ee311c12968b7e86d205fa;hp=ff4e1a92bbd36b2defaa33407d877a44c38a914d;hpb=a45774ee31d9f35d4eff46d54d7deab719afb092;p=jalview.git diff --git a/src/jalview/util/Comparison.java b/src/jalview/util/Comparison.java index ff4e1a9..08d2c88 100644 --- a/src/jalview/util/Comparison.java +++ b/src/jalview/util/Comparison.java @@ -1,19 +1,20 @@ /* - * Jalview - A Sequence Alignment Editor and Viewer (Version 2.7) - * Copyright (C) 2011 J Procter, AM Waterhouse, G Barton, M Clamp, S Searle + * Jalview - A Sequence Alignment Editor and Viewer (Version 2.8.2) + * Copyright (C) 2014 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.util; @@ -133,19 +134,28 @@ public class Comparison // Another pid with region specification public final static float PID(String seq1, String seq2, int start, int end) { - return PID(seq1, seq2, start, end, true,false); + return PID(seq1, seq2, start, end, true, false); } + /** - * Calculate percent identity for a pair of sequences over a particular range, with different options for ignoring gaps. + * Calculate percent identity for a pair of sequences over a particular range, + * with different options for ignoring gaps. + * * @param seq1 * @param seq2 - * @param start - position in seqs - * @param end - position in seqs - * @param wcGaps - if true - gaps match any character, if false, do not match anything - * @param ungappedOnly - if true - only count PID over ungapped columns + * @param start + * - position in seqs + * @param end + * - position in seqs + * @param wcGaps + * - if true - gaps match any character, if false, do not match + * anything + * @param ungappedOnly + * - if true - only count PID over ungapped columns * @return */ - public final static float PID(String seq1, String seq2, int start, int end, boolean wcGaps, boolean ungappedOnly) + public final static float PID(String seq1, String seq2, int start, + int end, boolean wcGaps, boolean ungappedOnly) { int s1len = seq1.length(); int s2len = seq2.length(); @@ -162,7 +172,7 @@ public class Comparison start = len - 1; // we just use a single residue for the difference } - int elen=len-start,bad = 0; + int elen = len - start, bad = 0; char chr1; char chr2; boolean agap; @@ -184,24 +194,31 @@ public class Comparison // Faster than toUpperCase chr2 -= caseShift; } - + if (chr1 != chr2) { - if (agap) - { - if (ungappedOnly) - { - elen--; - } else if (!wcGaps) { - bad++; - } - } else { - bad++; - } + if (agap) + { + if (ungappedOnly) + { + elen--; + } + else if (!wcGaps) + { + bad++; + } + } + else + { + bad++; + } } - + + } + if (elen < 1) + { + return 0f; } - if (elen<1) { return 0f; } return ((float) 100 * (elen - bad)) / elen; }