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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: coord_set.cpp,v $
  4.  * PRODUCTION Revision 1000.2  2004/06/01 18:28:31  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.15
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: coord_set.cpp,v 1000.2 2004/06/01 18:28:31 gouriano Exp $
  10. * ===========================================================================
  11. *
  12. *                            PUBLIC DOMAIN NOTICE
  13. *               National Center for Biotechnology Information
  14. *
  15. *  This software/database is a "United States Government Work" under the
  16. *  terms of the United States Copyright Act.  It was written as part of
  17. *  the author's official duties as a United States Government employee and
  18. *  thus cannot be copyrighted.  This software/database is freely available
  19. *  to the public for use. The National Library of Medicine and the U.S.
  20. *  Government have not placed any restriction on its use or reproduction.
  21. *
  22. *  Although all reasonable efforts have been taken to ensure the accuracy
  23. *  and reliability of the software and data, the NLM and the U.S.
  24. *  Government do not and cannot warrant the performance or results that
  25. *  may be obtained by using this software or data. The NLM and the U.S.
  26. *  Government disclaim all warranties, express or implied, including
  27. *  warranties of performance, merchantability or fitness for any particular
  28. *  purpose.
  29. *
  30. *  Please cite the author in any work or product based on this material.
  31. *
  32. * ===========================================================================
  33. *
  34. * Authors:  Paul Thiessen
  35. *
  36. * File Description:
  37. *      Classes to hold sets of coordinates for atoms and features
  38. *
  39. * ===========================================================================
  40. */
  41. #ifdef _MSC_VER
  42. #pragma warning(disable:4018)   // disable signed/unsigned mismatch warning in MSVC
  43. #endif
  44. #include <ncbi_pch.hpp>
  45. #include <corelib/ncbistd.hpp>
  46. #include <objects/mmdb2/Model_coordinate_set.hpp>
  47. #include <objects/mmdb2/Coordinates.hpp>
  48. #include <objects/mmdb2/Surface_coordinates.hpp>
  49. #include <objects/mmdb2/Model_descr.hpp>
  50. #include <objects/mmdb3/Chem_graph_pntrs.hpp>
  51. #include <objects/mmdb2/Surface_coordinates.hpp>
  52. #include "coord_set.hpp"
  53. #include "atom_set.hpp"
  54. #include "object_3d.hpp"
  55. #include "cn3d_tools.hpp"
  56. USING_NCBI_SCOPE;
  57. USING_SCOPE(objects);
  58. BEGIN_SCOPE(Cn3D)
  59. CoordSet::CoordSet(StructureBase *parent,
  60.                    const CBiostruc_model::TModel_coordinates& modelCoords) :
  61.     StructureBase(parent), atomSet(NULL)
  62. {
  63.     // iterate SEQUENCE OF Model-coordinate-set
  64.     CBiostruc_model::TModel_coordinates::const_iterator j, je=modelCoords.end();
  65.     for (j=modelCoords.begin(); j!=je; ++j) {
  66.         const CModel_coordinate_set::C_Coordinates&
  67.             coordSet = j->GetObject().GetCoordinates();
  68.         if (coordSet.IsLiteral()) {
  69.             const CCoordinates& coords = coordSet.GetLiteral();
  70.             // coordinates of atoms
  71.             if (coords.IsAtomic()) {
  72.                 if (!atomSet)
  73.                     atomSet = new AtomSet(this, coords.GetAtomic());
  74.                 else
  75.                     ERRORMSG("confused by multiple atomic coords");
  76.             }
  77.             // coordinates of 3d-objects
  78.             else if (coords.IsSurface() && j->GetObject().IsSetDescr() &&
  79.                 j->GetObject().GetDescr().front().GetObject().IsOther_comment() &&
  80.                 coords.GetSurface().GetContents().IsResidues()) {
  81.                 const string& descr =
  82.                     j->GetObject().GetDescr().front().GetObject().GetOther_comment();
  83.                 // helix cylinder
  84.                 if (descr == "helix" && coords.GetSurface().GetSurface().IsCylinder()) {
  85.                     Helix3D *helix =
  86.                         new Helix3D(this, coords.GetSurface().GetSurface().GetCylinder(),
  87.                             coords.GetSurface().GetContents().GetResidues());
  88.                     if (helix->moleculeID == Object3D::VALUE_NOT_SET) {
  89.                         this->_RemoveChild(helix);
  90.                         delete helix;
  91.                         WARNINGMSG("bad helix");
  92.                     } else
  93.                         objectMap[helix->moleculeID].push_back(helix);
  94.                 // strand brick
  95.                 } else if (descr == "strand" && coords.GetSurface().GetSurface().IsBrick()) {
  96.                     Strand3D *strand =
  97.                         new Strand3D(this, coords.GetSurface().GetSurface().GetBrick(),
  98.                             coords.GetSurface().GetContents().GetResidues());
  99.                     if (strand->moleculeID == Object3D::VALUE_NOT_SET) {
  100.                         this->_RemoveChild(strand);
  101.                         delete strand;
  102.                         WARNINGMSG("bad strand");
  103.                     } else
  104.                         objectMap[strand->moleculeID].push_back(strand);
  105.                 }
  106.             }
  107.         }
  108.     }
  109. }
  110. END_SCOPE(Cn3D)
  111. /*
  112. * ---------------------------------------------------------------------------
  113. * $Log: coord_set.cpp,v $
  114. * Revision 1000.2  2004/06/01 18:28:31  gouriano
  115. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.15
  116. *
  117. * Revision 1.15  2004/05/21 21:41:39  gorelenk
  118. * Added PCH ncbi_pch.hpp
  119. *
  120. * Revision 1.14  2004/03/15 18:23:01  thiessen
  121. * prefer prefix vs. postfix ++/-- operators
  122. *
  123. * Revision 1.13  2004/02/19 17:04:54  thiessen
  124. * remove cn3d/ from include paths; add pragma to disable annoying msvc warning
  125. *
  126. * Revision 1.12  2003/02/03 19:20:03  thiessen
  127. * format changes: move CVS Log to bottom of file, remove std:: from .cpp files, and use new diagnostic macros
  128. *
  129. * Revision 1.11  2001/05/31 18:54:21  thiessen
  130. * TESTMSG moved
  131. *
  132. * Revision 1.10  2001/02/08 23:01:50  thiessen
  133. * hook up C-toolkit stuff for threading; working PSSM calculation
  134. *
  135. * Revision 1.9  2000/08/27 18:52:21  thiessen
  136. * extract sequence information
  137. *
  138. * Revision 1.8  2000/08/16 14:18:44  thiessen
  139. * map 3-d objects to molecules
  140. *
  141. * Revision 1.7  2000/08/11 12:58:31  thiessen
  142. * added worm; get 3d-object coords from asn1
  143. *
  144. * Revision 1.6  2000/08/07 00:21:17  thiessen
  145. * add display list mechanism
  146. *
  147. * Revision 1.5  2000/08/03 15:12:23  thiessen
  148. * add skeleton of style and show/hide managers
  149. *
  150. * Revision 1.4  2000/07/27 13:30:51  thiessen
  151. * remove 'using namespace ...' from all headers
  152. *
  153. * Revision 1.3  2000/07/16 23:19:10  thiessen
  154. * redo of drawing system
  155. *
  156. * Revision 1.2  2000/07/12 23:27:49  thiessen
  157. * now draws basic CPK model
  158. *
  159. * Revision 1.1  2000/07/11 13:45:29  thiessen
  160. * add modules to parse chemical graph; many improvements
  161. *
  162. */