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

生物技术

开发平台:

C/C++

  1. %{
  2. /*  $Id: newick.ypp,v 1000.1 2004/06/01 18:09:36 gouriano Exp $
  3.  * ===========================================================================
  4.  *
  5.  *                            PUBLIC DOMAIN NOTICE
  6.  *               National Center for Biotechnology Information
  7.  *
  8.  *  This software/database is a "United States Government Work" under the
  9.  *  terms of the United States Copyright Act.  It was written as part of
  10.  *  the author's official duties as a United States Government employee and
  11.  *  thus cannot be copyrighted.  This software/database is freely available
  12.  *  to the public for use. The National Library of Medicine and the U.S.
  13.  *  Government have not placed any restriction on its use or reproduction.
  14.  *
  15.  *  Although all reasonable efforts have been taken to ensure the accuracy
  16.  *  and reliability of the software and data, the NLM and the U.S.
  17.  *  Government do not and cannot warrant the performance or results that
  18.  *  may be obtained by using this software or data. The NLM and the U.S.
  19.  *  Government disclaim all warranties, express or implied, including
  20.  *  warranties of performance, merchantability or fitness for any particular
  21.  *  purpose.
  22.  *
  23.  *  Please cite the author in any work or product based on this material.
  24.  *
  25.  * ===========================================================================
  26.  *
  27.  * Authors:  Josh Cherry
  28.  *
  29.  * File Description:  bison parser for Newick format phylogenetic trees
  30.  *
  31.  */
  32. /*
  33.  * Meant to be used in conjunction with flex lexer defined
  34.  * by newick.lpp.  Generate code using (essentially)
  35.  *   flex -olex.newick.cpp newick.lpp
  36.  *   bison -d -p newick newick.ypp
  37.  * The line '#include <unistd.h>' must be removed from the flex-generated
  38.  * code for it to compile on some platforms.
  39.  */
  40. #include <ncbi_pch.hpp>
  41. #include <algo/phy_tree/phy_node.hpp>
  42. USING_SCOPE(ncbi);
  43. static int nodecount = 0;
  44. static TPhyTreeNode *g_Tree = 0;
  45. void yyerror(const char *error);
  46. int yylex(void);
  47. CNcbiIstream *g_NewickIstr;
  48. // a list of all allocated nodes so they can be deleted on error
  49. vector<TPhyTreeNode*> g_NodeList;
  50. %}
  51. %union {
  52.     char *strval;
  53.     double dblval;
  54.     TPhyTreeNode *nodeval;
  55. }
  56. %token <strval> LABEL
  57. %token <dblval> NUM
  58. %type  <nodeval> node
  59. %type  <nodeval> nodelist
  60. %%
  61. tree: node ';' { 
  62.     g_Tree = $1;
  63. }
  64. ;
  65. node: '(' nodelist ')' ':' NUM { 
  66.     $$ = $2;
  67.     $$->GetValue().SetDist($5);
  68.     $$->GetValue().SetId(nodecount);
  69.     nodecount++; 
  70. }
  71. |  '(' nodelist ')' { 
  72.     $$ = $2;
  73.     $$->GetValue().SetId(nodecount);
  74.     nodecount++; 
  75. }
  76. | LABEL ':' NUM { 
  77.     $$ = new TPhyTreeNode;
  78.     $$->GetValue().SetLabel($1);
  79.     $$->GetValue().SetDist($3);
  80.     $$->GetValue().SetId(nodecount);
  81.     nodecount++; 
  82. }
  83. | LABEL { 
  84.     $$ = new TPhyTreeNode;
  85.     $$->GetValue().SetLabel($1);
  86.     $$->GetValue().SetId(nodecount);
  87.     nodecount++; 
  88. }
  89. ;
  90. nodelist:  node {
  91.     $$ = new TPhyTreeNode;
  92.     $$->AddNode($1);
  93. }
  94. | nodelist ',' node {
  95.     $$ = $1;
  96.     $$->AddNode($3);
  97. }
  98. ;
  99. %%
  100. void yyerror(const char *s)
  101. {
  102.     if (g_Tree) {
  103.         delete g_Tree;
  104.     } else {
  105.         // delete only those nodes that are nobody's children
  106.         vector<TPhyTreeNode*> del_list;
  107.         ITERATE (vector<TPhyTreeNode*>, iter, g_NodeList) {
  108.             if (!(*iter)->GetParent()) {
  109.                 del_list.push_back(*iter);
  110.             }
  111.         }
  112.         ITERATE (vector<TPhyTreeNode*>, iter, del_list) {
  113.             delete *iter;
  114.         }
  115.     }
  116.     g_NodeList.clear();
  117.     throw runtime_error(string("error parsing Newick format tree file:  ")
  118.                         + s);
  119. }
  120. void newick_flex_reset(void);
  121. BEGIN_NCBI_SCOPE
  122. TPhyTreeNode *ReadNewickTree(CNcbiIstream& is)
  123. {
  124.     g_NewickIstr = &is;
  125.     g_Tree = 0;
  126.     g_NodeList.clear();
  127.     newick_flex_reset();
  128.     yyparse();
  129.     if (!g_Tree) {
  130.         yyerror("unknown");
  131.     }
  132.     g_NodeList.clear();
  133.     return g_Tree;
  134. }
  135. END_NCBI_SCOPE
  136. /*
  137.  * ===========================================================================
  138.  * $Log: newick.ypp,v $
  139.  * Revision 1000.1  2004/06/01 18:09:36  gouriano
  140.  * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.2
  141.  *
  142.  * Revision 1.2  2004/05/25 14:42:45  jcherry
  143.  * #include <ncbi_pch.hpp>
  144.  *
  145.  * Revision 1.1  2004/02/11 17:40:56  jcherry
  146.  * Initial version
  147.  *
  148.  * ===========================================================================
  149.  */