seqdbgeneral.hpp
上传用户:yhdzpy8989
上传日期:2007-06-13
资源大小:13604k
文件大小:5k
源码类别:

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: seqdbgeneral.hpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 19:46:42  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. #ifndef CORELIB__SEQDB__SEQDBGENERAL_HPP
  10. #define CORELIB__SEQDB__SEQDBGENERAL_HPP
  11. /*  $Id: seqdbgeneral.hpp,v 1000.1 2004/06/01 19:46:42 gouriano Exp $
  12.  * ===========================================================================
  13.  *
  14.  *                            PUBLIC DOMAIN NOTICE
  15.  *               National Center for Biotechnology Information
  16.  *
  17.  *  This software/database is a "United States Government Work" under the
  18.  *  terms of the United States Copyright Act.  It was written as part of
  19.  *  the author's official duties as a United States Government employee and
  20.  *  thus cannot be copyrighted.  This software/database is freely available
  21.  *  to the public for use. The National Library of Medicine and the U.S.
  22.  *  Government have not placed any restriction on its use or reproduction.
  23.  *
  24.  *  Although all reasonable efforts have been taken to ensure the accuracy
  25.  *  and reliability of the software and data, the NLM and the U.S.
  26.  *  Government do not and cannot warrant the performance or results that
  27.  *  may be obtained by using this software or data. The NLM and the U.S.
  28.  *  Government disclaim all warranties, express or implied, including
  29.  *  warranties of performance, merchantability or fitness for any particular
  30.  *  purpose.
  31.  *
  32.  *  Please cite the author in any work or product based on this material.
  33.  *
  34.  * ===========================================================================
  35.  *
  36.  * Author:  Kevin Bealer
  37.  *
  38.  */
  39. #include <objtools/readers/seqdb/seqdbcommon.hpp>
  40. BEGIN_NCBI_SCOPE
  41. // Temporary development tools/tricks
  42. extern int seqdb_debug_class;
  43. enum seqdb_debug_bits {
  44.     debug_rh    = 1,
  45.     debug_rhsum = 2,
  46.     debug_mvol  = 4,
  47.     debug_alias = 8,
  48.     debug_oid   = 16
  49. };
  50. #define ifdebug_rh    if (seqdb_debug_class & debug_rh)    cerr
  51. #define ifdebug_rhsum if (seqdb_debug_class & debug_rhsum) cerr
  52. #define ifdebug_mvol  if (seqdb_debug_class & debug_mvol)  cerr
  53. #define ifdebug_alias if (seqdb_debug_class & debug_alias) cerr
  54. #define ifdebug_oid   if (seqdb_debug_class & debug_oid)   cerr
  55. // Byte-order-nonspecific (long) versions
  56. template<typename T>
  57. inline T SeqDB_GetStdOrdUnaligned(const T * stdord_obj)
  58. {
  59.     unsigned char * stdord =
  60. (unsigned char*)(stdord_obj);
  61.     
  62.     unsigned char * pend = stdord + sizeof(T) - 1;
  63.     unsigned char * pcur = stdord;
  64.     
  65.     T retval = *pcur;
  66.     
  67.     while(pcur < pend) {
  68. retval <<= 8;
  69. retval += *++pcur;
  70.     }
  71.     
  72.     return retval;
  73. }
  74. template<typename T>
  75. inline T SeqDB_GetBrokenUnaligned(const T * stdord_obj)
  76. {
  77.     unsigned char * stdord =
  78. (unsigned char*)(stdord_obj);
  79.     
  80.     unsigned char * pend = stdord;
  81.     unsigned char * pcur = stdord + sizeof(T) - 1;
  82.     
  83.     T retval = *pcur;
  84.     
  85.     while(pcur > pend) {
  86. retval <<= 8;
  87. retval += *--pcur;
  88.     }
  89.     
  90.     return retval;
  91. }
  92. // Macro Predicates for binary qualities
  93. #define IS_POWER_OF_TWO(x)    ((x) & ((x)-1))
  94. #define ALIGNED_TO_POW2(x,y)  (! ((x) & (-y)))
  95. #define PTR_ALIGNED_TO_SELF_SIZE(x) 
  96.     (IS_POWER_OF_TWO(sizeof(*x)) && ALIGNED_TO_POW2(size_t(x), sizeof(*x)))
  97. // Portable byte swapping from marshalled version
  98. #ifdef WORDS_BIGENDIAN
  99. template<typename T>
  100. inline T SeqDB_GetStdOrd(const T * stdord_obj)
  101. {
  102.     if (PTR_ALIGNED_TO_SELF_SIZE(stdord_obj)) {
  103.         return *stdord_obj;
  104.     } else {
  105.         return SeqDB_GetStdOrdUnaligned(stdord_obj);
  106.     }
  107. }
  108. template<typename T>
  109. inline T SeqDB_GetBroken(const T * stdord_obj)
  110. {
  111.     return SeqDB_GetBrokenUnaligned(stdord_obj);
  112. }
  113. #else
  114. template<typename T>
  115. inline T SeqDB_GetStdOrd(const T * stdord_obj)
  116. {
  117.     return SeqDB_GetStdOrdUnaligned(stdord_obj);
  118. }
  119. template<typename T>
  120. inline T SeqDB_GetBroken(const T * stdord_obj)
  121.     if (PTR_ALIGNED_TO_SELF_SIZE(stdord_obj)) {
  122.         return *stdord_obj;
  123.     } else {
  124.         return SeqDB_GetBrokenUnaligned(stdord_obj);
  125.     }
  126. }
  127. #endif
  128. // Combine two paths, trying a little to avoid duplicated delimiters.
  129. // If either string is empty, the other is returned.  Conceptually,
  130. // the first path is "cwd" and the second path is the filename.  So,
  131. // if the second path starts with "/", the first path is ignored.
  132. string SeqDB_CombinePath(const string & path, const string & file);
  133. // Returns a path minus filename (actually, last component).
  134. string SeqDB_GetDirName (string s);
  135. // Returns a filename minus greedy path.
  136. string SeqDB_GetFileName(string s);
  137. // Returns a filename minus non-greedy extension.
  138. string SeqDB_GetBasePath(string s);
  139. // Composition of the above two functions; returns a filename minus
  140. // greedy path and non-greedy extension.
  141. string SeqDB_GetBaseName(string s);
  142. // Find the full name, minus extension, of a ".?al" or ".?in" file,
  143. // and return it.  If not found, return null.
  144. string SeqDB_FindBlastDBPath(const string & file_name, char dbtype);
  145. END_NCBI_SCOPE
  146. #endif // CORELIB__SEQDB__SEQDBGENERAL_HPP