blast_rps_options.cpp
上传用户:yhdzpy8989
上传日期:2007-06-13
资源大小:13604k
文件大小:5k
- /*
- * ===========================================================================
- * PRODUCTION $Log: blast_rps_options.cpp,v $
- * PRODUCTION Revision 1000.2 2004/06/01 18:05:49 gouriano
- * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.6
- * PRODUCTION
- * ===========================================================================
- */
- /* $Id: blast_rps_options.cpp,v 1000.2 2004/06/01 18:05:49 gouriano Exp $
- * ===========================================================================
- *
- * PUBLIC DOMAIN NOTICE
- * National Center for Biotechnology Information
- *
- * This software/database is a "United States Government Work" under the
- * terms of the United States Copyright Act. It was written as part of
- * the author's official duties as a United States Government employee and
- * thus cannot be copyrighted. This software/database is freely available
- * to the public for use. The National Library of Medicine and the U.S.
- * Government have not placed any restriction on its use or reproduction.
- *
- * Although all reasonable efforts have been taken to ensure the accuracy
- * and reliability of the software and data, the NLM and the U.S.
- * Government do not and cannot warrant the performance or results that
- * may be obtained by using this software or data. The NLM and the U.S.
- * Government disclaim all warranties, express or implied, including
- * warranties of performance, merchantability or fitness for any particular
- * purpose.
- *
- * Please cite the author in any work or product based on this material.
- *
- * ===========================================================================
- *
- * Authors: Tom Madden
- *
- */
- /// @file blast_rps_options.cpp
- /// Implements the CBlastRPSOptionsHandle class.
- #include <ncbi_pch.hpp>
- #include <algo/blast/api/blast_rps_options.hpp>
- #include <objects/seqloc/Na_strand.hpp>
- #include "blast_setup.hpp"
- /** @addtogroup AlgoBlast
- *
- * @{
- */
- BEGIN_NCBI_SCOPE
- BEGIN_SCOPE(blast)
- CBlastRPSOptionsHandle::CBlastRPSOptionsHandle(EAPILocality locality)
- : CBlastOptionsHandle(locality)
- {
- if (m_Opts->GetLocality() == CBlastOptions::eRemote) {
- return;
- }
- SetDefaults();
- m_Opts->SetProgram(eRPSBlast);
- }
- void
- CBlastRPSOptionsHandle::SetLookupTableDefaults()
- {
- m_Opts->SetLookupTableType(RPS_LOOKUP_TABLE);
- }
- void
- CBlastRPSOptionsHandle::SetQueryOptionDefaults()
- {
- SetFilterString("F");
- m_Opts->SetStrandOption(objects::eNa_strand_unknown);
- }
- void
- CBlastRPSOptionsHandle::SetInitialWordOptionsDefaults()
- {
- SetXDropoff(BLAST_UNGAPPED_X_DROPOFF_PROT);
- SetWindowSize(BLAST_WINDOW_SIZE_PROT);
- // FIXME: extend_word_method is missing
- m_Opts->SetUngappedExtension();
- }
- void
- CBlastRPSOptionsHandle::SetGappedExtensionDefaults()
- {
- SetGapXDropoff(BLAST_GAP_X_DROPOFF_PROT);
- SetGapXDropoffFinal(BLAST_GAP_X_DROPOFF_FINAL_PROT);
- SetGapTrigger(BLAST_GAP_TRIGGER_PROT);
- m_Opts->SetGapExtnAlgorithm(eDynProgExt);
- }
- void
- CBlastRPSOptionsHandle::SetScoringOptionsDefaults()
- {
- SetGappedMode();
- m_Opts->SetMatrixName("BLOSUM62");
- // set invalid values for options that are not applicable
- m_Opts->SetOutOfFrameMode(false);
- m_Opts->SetFrameShiftPenalty(INT2_MAX);
- m_Opts->SetDecline2AlignPenalty(INT2_MAX);
- }
- void
- CBlastRPSOptionsHandle::SetHitSavingOptionsDefaults()
- {
- SetHitlistSize(500);
- SetPrelimHitlistSize(550);
- SetEvalueThreshold(BLAST_EXPECT_VALUE);
- SetPercentIdentity(0);
- m_Opts->SetSumStatisticsMode(false);
- // set some default here, allow INT4MAX to mean infinity
- SetMaxNumHspPerSequence(0);
- // this is never used... altough it could be calculated
- //SetTotalHspLimit(FIXME);
- SetCutoffScore(0); // will be calculated based on evalue threshold,
- // effective lengths and Karlin-Altschul params in BLAST_Cutoffs_simple
- // and passed to the engine in the params structure
- // not applicable
- m_Opts->SetRequiredStart(0);
- m_Opts->SetRequiredEnd(0);
- }
- void
- CBlastRPSOptionsHandle::SetEffectiveLengthsOptionsDefaults()
- {
- SetDbLength(0);
- SetDbSeqNum(1);
- SetEffectiveSearchSpace(0);
- SetUseRealDbSize();
- }
- void
- CBlastRPSOptionsHandle::SetSubjectSequenceOptionsDefaults()
- {}
- END_SCOPE(blast)
- END_NCBI_SCOPE
- /* @} */
- /*
- * ===========================================================================
- * $Log: blast_rps_options.cpp,v $
- * Revision 1000.2 2004/06/01 18:05:49 gouriano
- * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.6
- *
- * Revision 1.6 2004/05/21 21:41:02 gorelenk
- * Added PCH ncbi_pch.hpp
- *
- * Revision 1.5 2004/05/17 15:32:39 madden
- * Int algorithm_type replaced with enum EBlastPrelimGapExt
- *
- * Revision 1.4 2004/04/23 13:50:16 papadopo
- * derived BlastRPSOptionsHandle from BlastOptions (again)
- *
- * Revision 1.3 2004/04/16 14:27:47 papadopo
- * make this class a derived class of CBlastProteinOptionsHandle, that corresponds to the eRPSBlast program
- *
- * Revision 1.2 2004/03/19 15:13:34 camacho
- * Move to doxygen group AlgoBlast
- *
- * Revision 1.1 2004/03/10 14:52:13 madden
- * Options handle for RPSBlast searches
- *
- *
- * ===========================================================================
- */