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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: object_3d.cpp,v $
  4.  * PRODUCTION Revision 1000.2  2004/06/01 18:28:50  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.9
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: object_3d.cpp,v 1000.2 2004/06/01 18:28:50 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. *      holds 3d-objects - helix cylinders and strand bricks
  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/mmdb3/Residue_interval_pntr.hpp>
  47. #include <objects/mmdb3/Model_space_point.hpp>
  48. #include <objects/mmdb1/Molecule_id.hpp>
  49. #include <objects/mmdb1/Residue_id.hpp>
  50. #include "object_3d.hpp"
  51. #include "opengl_renderer.hpp"
  52. #include "style_manager.hpp"
  53. #include "structure_set.hpp"
  54. #include "cn3d_tools.hpp"
  55. USING_NCBI_SCOPE;
  56. USING_SCOPE(objects);
  57. BEGIN_SCOPE(Cn3D)
  58. const int Object3D::VALUE_NOT_SET = -1;
  59. Object3D::Object3D(StructureBase *parent, const CResidue_pntrs& residues) :
  60.     StructureBase(parent),
  61.     moleculeID(VALUE_NOT_SET), fromResidueID(VALUE_NOT_SET), toResidueID(VALUE_NOT_SET)
  62. {
  63.     if (!residues.IsInterval() || residues.GetInterval().size() != 1) {
  64.         ERRORMSG("Object3D::Object3D() - can't handle this type of Residue-pntrs (yet)!");
  65.         return;
  66.     }
  67.     const CResidue_interval_pntr& interval = residues.GetInterval().front().GetObject();
  68.     moleculeID = interval.GetMolecule_id().Get();
  69.     fromResidueID = interval.GetFrom().Get();
  70.     toResidueID = interval.GetTo().Get();
  71. }
  72. static inline void ModelPoint2Vector(const CModel_space_point& msp, Vector *v)
  73. {
  74.     v->x = msp.GetX();
  75.     v->y = msp.GetY();
  76.     v->z = msp.GetZ();
  77.     *v /= msp.GetScale_factor();
  78. }
  79. Helix3D::Helix3D(StructureBase *parent, const CCylinder& cylinder, const CResidue_pntrs& residues) :
  80.     Object3D(parent, residues)
  81. {
  82.     if (moleculeID == Object3D::VALUE_NOT_SET) return;
  83.     ModelPoint2Vector(cylinder.GetAxis_bottom(), &Nterm);
  84.     ModelPoint2Vector(cylinder.GetAxis_top(), &Cterm);
  85. }
  86. bool Helix3D::Draw(const AtomSet *data) const
  87. {
  88.     if (!parentSet->renderer) {
  89.         ERRORMSG("Helix3D::Draw() - no renderer");
  90.         return false;
  91.     }
  92.     // get object parent
  93.     const StructureObject *object;
  94.     if (!GetParentOfType(&object)) return false;
  95.     // get Style
  96.     HelixStyle helixStyle;
  97.     if (!parentSet->styleManager->GetHelixStyle(object, *this, &helixStyle))
  98.         return false;
  99.     // draw the Helix
  100.     if (helixStyle.style != StyleManager::eNotDisplayed)
  101.         parentSet->renderer->DrawHelix(Nterm, Cterm, helixStyle);
  102.     return true;
  103. }
  104. Strand3D::Strand3D(StructureBase *parent, const CBrick& brick, const CResidue_pntrs& residues) :
  105.     Object3D(parent, residues)
  106. {
  107.     if (moleculeID == Object3D::VALUE_NOT_SET) return;
  108.     Vector c1, c2;
  109.     ModelPoint2Vector(brick.GetCorner_000(), &c1);
  110.     ModelPoint2Vector(brick.GetCorner_011(), &c2);
  111.     Nterm = (c1 + c2) / 2;
  112.     ModelPoint2Vector(brick.GetCorner_100(), &c1);
  113.     ModelPoint2Vector(brick.GetCorner_111(), &c2);
  114.     Cterm = (c1 + c2) / 2;
  115.     ModelPoint2Vector(brick.GetCorner_010(), &c1);
  116.     ModelPoint2Vector(brick.GetCorner_000(), &c2);
  117.     unitNormal = c1 - c2;
  118.     unitNormal.normalize();
  119. }
  120. bool Strand3D::Draw(const AtomSet *data) const
  121. {
  122.     if (!parentSet->renderer) {
  123.         ERRORMSG("Strand3D::Draw() - no renderer");
  124.         return false;
  125.     }
  126.     // get object parent
  127.     const StructureObject *object;
  128.     if (!GetParentOfType(&object)) return false;
  129.     // get Style
  130.     StrandStyle strandStyle;
  131.     if (!parentSet->styleManager->GetStrandStyle(object, *this, &strandStyle))
  132.         return false;
  133.     // draw the Strand
  134.     if (strandStyle.style != StyleManager::eNotDisplayed)
  135.         parentSet->renderer->DrawStrand(Nterm, Cterm, unitNormal, strandStyle);
  136.     return true;
  137. }
  138. END_SCOPE(Cn3D)
  139. /*
  140. * ---------------------------------------------------------------------------
  141. * $Log: object_3d.cpp,v $
  142. * Revision 1000.2  2004/06/01 18:28:50  gouriano
  143. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.9
  144. *
  145. * Revision 1.9  2004/05/21 21:41:39  gorelenk
  146. * Added PCH ncbi_pch.hpp
  147. *
  148. * Revision 1.8  2004/02/19 17:04:59  thiessen
  149. * remove cn3d/ from include paths; add pragma to disable annoying msvc warning
  150. *
  151. * Revision 1.7  2003/02/03 19:20:04  thiessen
  152. * format changes: move CVS Log to bottom of file, remove std:: from .cpp files, and use new diagnostic macros
  153. *
  154. * Revision 1.6  2001/12/12 14:04:14  thiessen
  155. * add missing object headers after object loader change
  156. *
  157. * Revision 1.5  2001/02/08 23:01:50  thiessen
  158. * hook up C-toolkit stuff for threading; working PSSM calculation
  159. *
  160. * Revision 1.4  2000/08/17 14:24:06  thiessen
  161. * added working StyleManager
  162. *
  163. * Revision 1.3  2000/08/16 14:18:45  thiessen
  164. * map 3-d objects to molecules
  165. *
  166. * Revision 1.2  2000/08/13 02:43:01  thiessen
  167. * added helix and strand objects
  168. *
  169. * Revision 1.1  2000/08/11 18:24:57  thiessen
  170. * add 3-d objects code
  171. *
  172. */