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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: test_asn_converter.cpp,v $
  4.  * PRODUCTION Revision 1000.2  2004/06/01 19:14:47  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: test_asn_converter.cpp,v 1000.2 2004/06/01 19:14:47 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. * Author:  Aaron Ucko
  35. *
  36. * File Description:
  37. *   Test of <ctools/asn_converter.hpp>
  38. *
  39. * ===========================================================================
  40. */
  41. #include <ncbi_pch.hpp>
  42. #include <corelib/ncbiapp.hpp>
  43. #include <corelib/ncbiargs.hpp>
  44. #include <corelib/ncbienv.hpp>
  45. #include <ctools/asn_converter.hpp>
  46. #include <objects/seqalign/Score.hpp>
  47. #include <objects/seqalign/Seq_align.hpp>
  48. #include <objects/seqalign/Std_seg.hpp>
  49. #include <objects/seqloc/Seq_loc.hpp>
  50. #include <objects/seqset/Seq_entry.hpp>
  51. #include <serial/serial.hpp>
  52. #include <objsset.h>
  53. #include <test/test_assert.h>
  54. USING_NCBI_SCOPE;
  55. USING_SCOPE(objects);
  56. class CTestAsnConverterApp : public CNcbiApplication
  57. {
  58. public:
  59.     void Init();
  60.     int  Run();
  61. };
  62. void CTestAsnConverterApp::Init()
  63. {
  64.     auto_ptr<CArgDescriptions> arg_desc(new CArgDescriptions);
  65.     arg_desc->SetUsageContext
  66.         (GetArguments().GetProgramBasename(),
  67.          "Test of converting objects between C and C++ representations.");
  68.     arg_desc->AddDefaultKey
  69.         ("type", "AsnType", "ASN.1 type of test object",
  70.          CArgDescriptions::eString, "Seq-entry");
  71.     arg_desc->SetConstraint
  72.         ("type", &(*new CArgAllow_Strings, "Seq-entry", "Seq-align", "Score"));
  73.     arg_desc->AddDefaultKey
  74.         ("in", "InputFile",
  75.          "name of file to read from (standard input by default)",
  76.          CArgDescriptions::eString, "-");
  77.     arg_desc->AddDefaultKey("infmt", "InputFormat", "format of input file",
  78.                             CArgDescriptions::eString, "asn");
  79.     arg_desc->SetConstraint
  80.         ("infmt", &(*new CArgAllow_Strings, "asn", "asnb", "xml"));
  81.     SetupArgDescriptions(arg_desc.release());
  82. }
  83. static ESerialDataFormat s_GetFormat(const string& name)
  84. {
  85.     if (name == "asn") {
  86.         return eSerial_AsnText;
  87.     } else if (name == "asnb") {
  88.         return eSerial_AsnBinary;
  89.     } else if (name == "xml") {
  90.         return eSerial_Xml;
  91.     } else {
  92.         THROW1_TRACE(runtime_error, "Bad serial format name " + name);
  93.     }
  94. }
  95. int CTestAsnConverterApp::Run()
  96. {
  97.     CArgs args = GetArgs();
  98.     auto_ptr<CObjectIStream> in;
  99.     if (args["type"].AsString() != "Score") {
  100.         in.reset(CObjectIStream::Open
  101.                  (s_GetFormat(args["infmt"].AsString()),
  102.                   args["in"].AsString(), eSerial_StdWhenDash));
  103.     }
  104.     if (args["type"].AsString() == "Seq-entry") {
  105.         DECLARE_ASN_CONVERTER(CSeq_entry, SeqEntry, convertor);
  106.         CSeq_entry  cpp1, cpp2;
  107.         SeqEntryPtr c1,   c2;
  108.         *in >> cpp1;
  109.         c1 = convertor.ToC(cpp1);
  110.         convertor.FromC(c1, &cpp2);
  111.         c2 = convertor.ToC(cpp2);
  112.         _ASSERT(AsnIoMemComp(c1, c2, (AsnWriteFunc)SeqEntryAsnWrite) == TRUE);
  113.         _ASSERT(cpp1.Equals(cpp2));
  114.     } else if (args["type"].AsString() == "Seq-align") {
  115.         DECLARE_ASN_CONVERTER(CSeq_align, SeqAlign, convertor);
  116.         CSeq_align  cpp1, cpp2;
  117.         SeqAlignPtr c1,   c2;
  118.         *in >> cpp1;
  119.         c1 = convertor.ToC(cpp1);
  120.         convertor.FromC(c1, &cpp2);
  121.         c2 = convertor.ToC(cpp2);
  122.         _ASSERT(AsnIoMemComp(c1, c2, (AsnWriteFunc)SeqAlignAsnWrite) == TRUE);
  123.         _ASSERT(cpp1.Equals(cpp2));
  124.     } else if (args["type"].AsString() == "Score") {
  125.         // special case: treat "in" as a score...
  126.         DECLARE_ASN_CONVERTER(CSeq_align, SeqAlign, convertor);
  127.         CSeq_align  cpp1, cpp2;
  128.         SeqAlignPtr c1,   c2;
  129.         cpp1.SetType(CSeq_align::eType_global);
  130.         {{
  131.             double d = NStr::StringToDouble(args["in"].AsString());
  132.             clog << setprecision(DBL_DIG) << d << endl;
  133.             CRef<CScore> score(new CScore);
  134.             score->SetValue().SetReal(d);
  135.             cpp1.SetScore().push_back(score);
  136.         }}
  137.         {{
  138.             CRef<CStd_seg> dummy_seg(new CStd_seg);
  139.             CRef<CSeq_loc> dummy_loc(new CSeq_loc);
  140.             dummy_loc->SetNull();
  141.             dummy_seg->SetLoc().push_back(dummy_loc);
  142.             cpp1.SetSegs().SetStd().push_back(dummy_seg);
  143.         }}
  144.         c1 = convertor.ToC(cpp1);
  145.         convertor.FromC(c1, &cpp2);
  146.         c2 = convertor.ToC(cpp2);
  147.         clog << cpp1.GetScore().front()->GetValue().GetReal() << ' '
  148.              << c1->score->value.realvalue << ' '
  149.              << cpp2.GetScore().front()->GetValue().GetReal() << ' '
  150.              << c2->score->value.realvalue << endl;
  151.         _ASSERT(AsnIoMemComp(c1, c2, (AsnWriteFunc)SeqAlignAsnWrite) == TRUE);
  152.         _ASSERT(cpp1.Equals(cpp2));
  153.     } else {
  154.         cerr << "Unrecognized ASN.1 type " + args["type"].AsString() << endl;
  155.         return 1;
  156.     }
  157.     cout << "Test succeeded." << endl;
  158.     return 0;
  159. }
  160. int main(int argc, const char* argv[])
  161. {
  162.     // Execute main application function
  163.     return CTestAsnConverterApp().AppMain(argc, argv, 0, eDS_Default, 0);
  164. }
  165. /*
  166. * ===========================================================================
  167. *
  168. * $Log: test_asn_converter.cpp,v $
  169. * Revision 1000.2  2004/06/01 19:14:47  gouriano
  170. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.3
  171. *
  172. * Revision 1.3  2004/05/17 20:59:15  gorelenk
  173. * Added include of PCH ncbi_pch.hpp
  174. *
  175. * Revision 1.2  2003/12/23 19:37:30  ucko
  176. * Add a Score test that treats IN as the real number to try converting.
  177. *
  178. * Revision 1.1  2002/08/08 18:18:46  ucko
  179. * Add test for <ctools/asn_converter.hpp>
  180. *
  181. *
  182. * ===========================================================================
  183. */