开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 3G开发 > 查看源码
syst_graph.cpp
资源名称：simulation_wireless_c++.rar　[点击查看]
上传用户：jtjnyq9001
上传日期：2014-11-21
资源大小：3974k
文件大小：41k
源码类别：

3G开发

开发平台：

Visual C++

syst_graph.cpp：源码内容
  1. //
  2. //  File = syst_graph.cpp
  3. //
  5. #include <stdlib.h>
  6. #include <fstream>
  7. #include "syst_graph.h"
  8. #include "globals.h"
  9. #include "sigplot.h"
  10. //#include "sig_mgr.h"
  11. //#include "cntl_mgr.h"
  12. //#include "actsysgp.h"
  13. //#include "cmsg.h"
  14. //#include "reinit.h"
  16. #ifdef _DEBUG
  17.   extern ofstream *DebugFile;
  18. #endif
  20. extern SignalPlotter SigPlot;
  21. //extern Reinitializer Reinit;
  22. extern int PassNumber;
  23. extern int EnclaveNumber;
  24. extern int EnclaveOffset[10];
  25. extern PracSimModel *ActiveModel;
  27. //============================================
  28. // constructor
  30. SystemGraph::SystemGraph()
  31. {
  32.   Sdg_Vert_Descr = new std::vector<sdg_sig_desc_type*>;
  33.   Sdg_Edge_Descr = new std::vector<sdg_edge_desc_type*>;
  34.   Samp_Intvl = new std::vector<double>;
  35.   Block_Size = new std::vector<int>;
  36.   Syst_Lev_Models = new std::vector<PracSimModel*>;
  37.   Sig_Dep_Graph = new DirectedGraph;
  38.   Num_Regular_Sigs = 0;
  39.   return;
  40. }
  41. //============================================
  42. // destructor
  44. SystemGraph::~SystemGraph()
  45. {
  46. }
  47. //============================================
  49. void SystemGraph::MergeCurrModelGraph( ModelGraph* curr_mod_graph )
  50. {
  51.    std::vector<int> merged_verts;
  52.    int num_cmg_verts;
  53.    int cmg_vert_num;
  54.    int sdg_sig_num;
  55.    obj_id_type sdg_sig_id, cmg_sig_id;
  56.    SignalKinds_type cmg_vert_kind;
  57.    bool node_not_found;
  58.    double cmg_samp_intvl;
  59.    int cmg_block_size;
  61.    //------------------------------------------------------------
  62.    // First section of code merges the vertices from the current
  63.    // model graph (CMG) into the system dependency graph (SDG)
  64.    //
  65.    //  Check each signal from CSG to see if it's already in SDG
  67.    num_cmg_verts = curr_mod_graph->GetNumVerts();
  69.    //-----------------------------------------------------
  70.    //  This loop checks each signal from CMG to see if it's 
  71.    //  already in SDG
  73.    for(  cmg_vert_num = 0; 
  74.          cmg_vert_num < num_cmg_verts;
  75.          cmg_vert_num++) {
  76.       cmg_vert_kind = curr_mod_graph->GetVertexKind(cmg_vert_num);
  77.       cmg_sig_id = curr_mod_graph->GetVertexId(cmg_vert_num);
  79.       node_not_found = true;
  80.       int num_sdg_verts = Sig_Dep_Graph->GetNumVerts();
  82.       for( sdg_sig_num = 0; 
  83.            (sdg_sig_num < num_sdg_verts) & node_not_found; 
  84.            sdg_sig_num++) {
  86.          sdg_sig_id = Sig_Dep_Graph->GetVertexId(sdg_sig_num);
  88.          if( sdg_sig_id == cmg_sig_id ) {
  89.             //---------------------------------------
  90.             // this CMSG signal is already in SDG, so add its
  91.             // index (in SDG) to list of merged signals
  93.             merged_verts.push_back(sdg_sig_num);
  95.             //----------------------------------------
  96.             // If sampling rate is undefined for signal in SDG,
  97.             // copy value from CMG signal (which may or may not be 
  98.             // defined).  If sampling rate is defined for signal in SDG,
  99.             // check to see if there is a different value defined in
  100.             // CMG.  If there is, a fatal error condition exists.
  102.             cmg_samp_intvl = curr_mod_graph->GetSampIntvl(cmg_vert_num);
  104.             if(((Sdg_Vert_Descr->at(sdg_sig_num))->samp_intvl) == 0.0) {
  105.                ((Sdg_Vert_Descr->at(sdg_sig_num))->samp_intvl) = cmg_samp_intvl;
  106.             }
  107.             else {
  108.                if( (cmg_samp_intvl != 0.0) &&
  109.                    (cmg_samp_intvl != 
  110.                    ((Sdg_Vert_Descr->at(sdg_sig_num))->samp_intvl)) ) {
  112.                   // -- ERROR --
  113.                   *DebugFile << "Fatal error -- "
  114.                              << "incomsistent sampling rate for signal "
  115.                              << ((GenericSignal*)sdg_sig_id)->GetName() << endl;
  116.                   exit(1);
  117.                }
  118.             }
  120.             //---------------------------------------------------
  121.             //  If block size is undefined for signal in SDG,
  122.             //  copy value from CMG signal (which may or may not be
  123.             //  defined.  If block size is defined for signal in SDG,
  124.             //  check to see if there is a different value in CMG.
  125.             //  If there is, a fatal error condition exists.
  127.             cmg_block_size = curr_mod_graph->GetBlockSize(cmg_vert_num);
  128.             if( ((Sdg_Vert_Descr->at(sdg_sig_num))->block_size) == 0 ) {
  129.                ((Sdg_Vert_Descr->at(sdg_sig_num))->block_size) = cmg_block_size;
  130.             }
  131.             else {
  132.                if( (cmg_block_size != 0 ) &&
  133.                    (cmg_block_size != 
  134.                    ((Sdg_Vert_Descr->at(sdg_sig_num))->block_size) )) {
  135.                // -- ERROR --
  136.                *DebugFile << "Fatal error -- "
  137.                           << "inconsistent block size for signal "
  138.                           << ((GenericSignal*)sdg_sig_id)->GetName() << endl;
  139.                exit(1);
  140.                }
  141.             }
  142.             // no need to check remaining SDG signals for this CMSG signal
  143.             node_not_found = false;
  145.          } // end of if( sdg_sig_id == cmg_sig_id )
  146.       } // end of loop over sdg_sig_num
  148.       //-----------------------------------------------------
  149.       // If loop above completed without a match, we need
  150.       // to add CMG vertex to SDG
  152.       if( node_not_found ) {
  154.          // maybe...we should use the sig_desc_type structure
  155.          // in the CMG, then this would become a simple
  156.          // structure to structure copy
  157.          sdg_sig_desc_type* new_sig_desc = new sdg_sig_desc_type;
  158.          new_sig_desc->signal_id = 
  159.                      curr_mod_graph->GetSignalId(cmg_vert_num);
  160.          new_sig_desc->block_size = 
  161.                      curr_mod_graph->GetBlockSize(cmg_vert_num);
  162.          new_sig_desc->samp_intvl = 
  163.                      curr_mod_graph->GetSampIntvl(cmg_vert_num);
  164.          new_sig_desc->kind_of_signal = 
  165.                      curr_mod_graph->GetVertexKind(cmg_vert_num);
  167.          if(new_sig_desc->kind_of_signal == SK_REGULAR_SIGNAL) Num_Regular_Sigs++;
  168.          Sdg_Vert_Descr->push_back(new_sig_desc);
  169.          sdg_sig_num = Sig_Dep_Graph->AddVertex(cmg_sig_id);
  170.          merged_verts.push_back(sdg_sig_num);
  171.       }
  172.    } // end of loop over cmg_vert_num
  175.    //----------------------------------------------------------
  176.    //  at this point we are done adding vertices, and now begin
  177.    //  to add edges between these vertices as specified in the CMG
  178.    //
  179.    //  new edges in the SDG can only connect to vertices that were
  180.    //  just added or which matched vertices in the CMG
  182.    int sdg_fm_idx, sdg_to_idx;
  183.    int cmg_edge_num;
  184.    int existing_edge;
  185.    bool conn_is_feedback;
  186.    obj_id_type cmg_edge_id;
  188.    int num_merged_verts;
  189.    num_merged_verts = int(merged_verts.size());
  190.    for(  int cmg_fm_idx = 0; 
  191.          cmg_fm_idx < num_merged_verts; 
  192.          cmg_fm_idx++)  {
  194.       sdg_fm_idx = merged_verts.at(cmg_fm_idx);
  195.       conn_is_feedback = curr_mod_graph->ConnIsFeedback(cmg_fm_idx);
  197.       for( int cmg_to_idx = 0; 
  198.            cmg_to_idx < num_merged_verts; 
  199.            cmg_to_idx++) {
  201.          sdg_to_idx = merged_verts.at(cmg_to_idx);
  202.          cmg_edge_num = curr_mod_graph->GetEdgeNum(cmg_fm_idx, cmg_to_idx);
  204.          if(cmg_edge_num >= 0) {
  205.             //---------------------------------
  206.             // CMG contains an edge from vertex 'cmg_fm_idx'
  207.             // to vertex 'cmg_to_idx'
  208.             // See if ASG-SDG already contains an edge between
  209.             // the corresponding vertices.
  211.             existing_edge = 
  212.                   Sig_Dep_Graph->GetEdgeNum(sdg_fm_idx, sdg_to_idx);
  214.             if(existing_edge >= 0 ) {
  215.                double cmg_resamp = curr_mod_graph->GetResampRate(cmg_edge_num);
  216.               // double sdg_resamp = Sig_Dep_Graph->GetResampRate(existing_edge);
  217.                double sdg_resamp = (Sdg_Edge_Descr->at(existing_edge))->resamp_rate;;
  219.                if(cmg_resamp != sdg_resamp)
  220.                {
  221.                   // parallel edge exists in ASG-SDG -- not sure how to
  222.                   // handle -- see old listing for handling of edges
  223.                   // due to subordinate models
  224.                   *DebugFile << "parallel edge error" << endl;
  225.                   *DebugFile << "existing_edge = " << existing_edge << endl;
  226.                   *DebugFile << "cmg_fm_idx = " << cmg_fm_idx << endl;
  227.                   *DebugFile << "cmg_to_idx = " << cmg_to_idx << endl;
  228.                   exit(1);
  229.                }
  230.             }
  231.             //-----------------------------------------
  232.             // OK to add edge
  233.       
  234.             cmg_edge_id = curr_mod_graph->GetEdgeId(cmg_edge_num);
  235.             Sig_Dep_Graph->AddEdge( cmg_edge_id, 
  236.                                     sdg_fm_idx, 
  237.                                     sdg_to_idx);
  239.             sdg_edge_desc_type* new_edge_desc = new sdg_edge_desc_type;
  240.             new_edge_desc->model_id = (PracSimModel*)cmg_edge_id;
  241.             new_edge_desc->resamp_rate = 
  242.                      curr_mod_graph->GetResampRate(cmg_edge_num);
  243.             new_edge_desc->is_const_intvl = 
  244.                      curr_mod_graph->GetIsConstIntvl(cmg_edge_num);
  245.             new_edge_desc->delta_delay =
  246.                      curr_mod_graph->GetDelay(cmg_edge_num);
  247.             Sdg_Edge_Descr->push_back(new_edge_desc);
  249.          } // end of if(cmg_edge_num)
  251.       } // end of loop over cmg_to_idx
  253.    } // end of loop over cmg_fm_idx
  254.    return;
  255. }
  256. //========================================================
  257. void SystemGraph::DumpSDGraph(void)
  258. {
  259.    int sig_num, sig_num2, edge_num;
  260.    GenericSignal *sig_id;
  261.    PracSimModel *model_id;
  262.    int num_nodes = Sig_Dep_Graph->GetNumVerts();
  263.    int num_edges = Sig_Dep_Graph->GetNumEdges();
  265.    *DebugFile << "nIn SystemGraph::DumpSDGraph" << endl;
  266.    *DebugFile << "Num_SDG_Nodes = " << num_nodes << endl;
  267.    *DebugFile << "Vertices:" << endl;
  269.    *DebugFile << "index, Sdg_Vtx_Id, Signal Name, Samp Intvl, "
  270.               << "Block Size" << endl;
  272.    for( sig_num = 0; sig_num < num_nodes; sig_num++) {
  273.       sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  274.       *DebugFile << sig_num << ")  "
  275.                  << sig_id << "  "
  276.                  << sig_id->GetName() << "  "
  277.                  << (Sdg_Vert_Descr->at(sig_num))->samp_intvl << "  "
  278.                  << (Sdg_Vert_Descr->at(sig_num))->block_size << endl;
  279.    }
  281.    *DebugFile << "nEdges:" << endl;
  282.    *DebugFile << "index, SDG edge model, Model Name, Resamp Rate, "
  283.               << "Delay" << endl;
  284.    for( edge_num = 0; edge_num < num_edges; edge_num ++) {
  285.       model_id = (Sdg_Edge_Descr->at(edge_num))->model_id;
  286.       *DebugFile << edge_num << ")  "
  287.                  << model_id << "  "
  288.                  << model_id->GetModelName() << ":" 
  289.                  << model_id->GetInstanceName() << "  "
  290.                  << (Sdg_Edge_Descr->at(edge_num))->resamp_rate << "  "
  291.                  << (Sdg_Edge_Descr->at(edge_num))->delta_delay << endl;
  292.    }
  294.    *DebugFile << "nAdjacency Matrix:" << endl;
  295.    for( sig_num = 0; sig_num < num_nodes; sig_num++) {
  296.       for( sig_num2 = 0; sig_num2 < num_nodes; sig_num2++) {
  297.          edge_num = Sig_Dep_Graph->GetEdgeNum(sig_num, sig_num2);
  299.          if( edge_num >= 0 ) {
  300.         *DebugFile << "Edge " << edge_num
  301.                    << " connects from vertex " << sig_num
  302.                    << " to vertex " << sig_num2 << endl;
  303.          }
  304.       }
  305.    }
  306. }
  307. //===================================================================
  309. void SystemGraph::ResolveSignalParms(void)
  310. {
  311.   int sig_num, unsorted_sig_num, base_sig_num, edge_num;
  312.   int num_base_sigs_used = 0;
  313.   bool base_sig_not_found;
  314.   double samp_rate;
  315.   int block_size;
  316.   GenericSignal *sig_id;
  317.   bool is_forward;
  319.   int num_nodes = Sig_Dep_Graph->GetNumVerts();
  320.   bool *used_as_base = new bool[num_nodes];
  322.   for( sig_num = 0; sig_num < num_nodes; sig_num++)
  323.     {
  324.     used_as_base[sig_num] = false;
  325.     }
  326.   this->TopoSortSDG();
  328.   while(num_base_sigs_used < Num_Regular_Sigs)
  329.     {
  330.     //---------------------------------------------------------------
  331.     // find base signal
  333.     base_sig_not_found = true;
  334.     for(unsorted_sig_num = 0; unsorted_sig_num < num_nodes; unsorted_sig_num++)
  335.       {
  336.       sig_num = Sorted_Sig_Nodes[unsorted_sig_num];
  338.       sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  339.       samp_rate = (Sdg_Vert_Descr->at(sig_num))->samp_intvl;
  340.       if(samp_rate <= 0.0)
  341.         {
  342.         #ifdef _DEBUG
  343.         *DebugFile << sig_id->GetName()
  344.                   << " rejected because samp rate <= 0 ( "
  345.                   << samp_rate << " )" << endl;
  346.         #endif
  347.         continue;
  348.         }
  350.       block_size = (Sdg_Vert_Descr->at(sig_num))->block_size;
  351.       if(block_size <= 0.0)
  352.         {
  353.         #ifdef _DEBUG
  354.         *DebugFile << sig_id->GetName()
  355.                   << " rejected because block size <= 0 ( "
  356.                   << block_size << " )" << endl;
  357.         #endif
  358.         continue;
  359.         }
  361.       if(used_as_base[sig_num])
  362.         {
  363.         #ifdef _DEBUG
  364.         *DebugFile << sig_id->GetName()
  365.                   << " rejected because already used as base"
  366.                   << endl;
  367.         #endif
  368.         continue;
  369.         }
  370.       //
  371.       // if execution gets to here, we have found a useable base signal
  372.       //
  373.       base_sig_num = sig_num;
  374.       used_as_base[base_sig_num] = true;
  375.       num_base_sigs_used++;
  376.       base_sig_not_found = false;
  377.       break;
  378.       } // end of loop over unsorted_sig_num
  380.     if(base_sig_not_found)
  381.       {
  382.       (*DebugFile) << "Unable to find suitable base signal" << endl;
  383.       #ifdef _DEBUG
  384.         *DebugFile << "The following signals have not"
  385.                   << " been used as base signals:" << endl;
  386.         for(unsorted_sig_num = 0; unsorted_sig_num < num_nodes; unsorted_sig_num++)
  387.           {
  388.           sig_num = Sorted_Sig_Nodes[unsorted_sig_num];
  389.           if(used_as_base[sig_num]) continue;
  390.             sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  391.             *DebugFile << "   " << sig_id->GetName() << endl;
  392.           }
  393.       #endif
  394.       exit(0);
  395.       }
  397.     #ifdef _DEBUG
  398.       sig_id = (Sdg_Vert_Descr->at(base_sig_num))->signal_id;
  399.       *DebugFile << "nbase signal is " << sig_id->GetName() << endl;
  400.     #endif
  402.     //---------------------------------------------------------
  403.     //  Work down the column corresponding to the base signal to
  404.     //  find inbound edges.  For each edge found, attempt to 
  405.     //  propagate parameters back to the edge's source signal
  407.     for( sig_num = 0; sig_num < num_nodes; sig_num++)
  408.       {
  409.       edge_num = Sig_Dep_Graph->GetEdgeNum( sig_num, base_sig_num );
  410.       if(edge_num < 0 ) continue;
  411.         //
  412.         //else
  413.         is_forward = false;
  414.         #ifdef _DEBUG
  415.           sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  416.           *DebugFile << "looking backward to " << sig_id->GetName() << endl;
  417.         #endif
  418.         Propagate( base_sig_num, edge_num, sig_num, is_forward);
  419.       }
  421.     //-------------------------------------------------------------
  422.     //  Work across the row corresponding to the base signal to
  423.     //  find outbound edges.  For each edge found, attempt to
  424.     //  propagate parameters forward to the edge's destination signal.
  426.     for( sig_num = 0; sig_num < num_nodes; sig_num++)
  427.       {
  428.       edge_num = Sig_Dep_Graph->GetEdgeNum( base_sig_num, sig_num );
  429.       if(edge_num < 0 ) continue;
  430.         //
  431.         //else
  432.         is_forward = true;
  433.         #ifdef _DEBUG
  434.           sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  435.           *DebugFile << "looking forward to " << sig_id->GetName() << endl;
  436.         #endif
  437.         Propagate( base_sig_num, edge_num, sig_num, is_forward);
  438.       }
  439.     } // end of while(num_base_sigs_used < Num_Regular_Sigs)
  440.   return;
  441. }
  442. //============================================================
  443. void SystemGraph::TopoSortSDG(void)
  444. {
  445.   int num_nodes = Sig_Dep_Graph->GetNumVerts();
  446.   Sorted_Sig_Nodes = new int[num_nodes];
  447.   for(int i=0; i<num_nodes; i++)
  448.     {
  449.     Sorted_Sig_Nodes[i] = i;
  450.     }
  451.   return;
  452. }
  453. //=======================================================================
  454. void SystemGraph::Propagate( int base_sig_num,
  455.                                    int edge_num,
  456.                                    int sig_num,
  457.                                    bool is_forward )
  458. {
  459.   double t_base, t_sig;
  460.   double r_edge;
  461.   int n_base, n_sig;
  462.   int r_def, t_def, n_def;
  463.   int config_key;
  464.   bool is_const_intvl;
  466.   t_base = (Sdg_Vert_Descr->at(base_sig_num))->samp_intvl;
  467.   n_base = (Sdg_Vert_Descr->at(base_sig_num))->block_size;
  468.   r_edge = (Sdg_Edge_Descr->at(edge_num))->resamp_rate;
  469.   is_const_intvl = (Sdg_Edge_Descr->at(edge_num))->is_const_intvl;
  470.   t_sig = (Sdg_Vert_Descr->at(sig_num))->samp_intvl;
  471.   n_sig = (Sdg_Vert_Descr->at(sig_num))->block_size;
  473.   r_def = 0;
  474.   t_def = 0;
  475.   n_def = 0;
  476.   if( r_edge > 0.0 ) r_def = 1;
  477.   if( t_sig > 0.0) t_def = 2;
  478.   if( n_sig >0 ) n_def = 4;
  480.   config_key = r_def + t_def + n_def;
  481.   if(is_forward) config_key +=8;
  482.   if(is_const_intvl) config_key +=16;
  484.   switch (config_key)
  485.     {
  486.     case 0: //backward, nothing defined
  487.       // no action possible
  488.       #ifdef _DEBUG
  489.         *DebugFile << "  case 0 - no action possible" << endl;
  490.       #endif
  491.       return;
  492.     case 1:  //backward, resampling rate is defined
  493.       // compute sampling rate : Tn = R * Tb
  494.       // compute block size : Nn = Nb / R
  496.       t_sig = r_edge * t_base;
  497.       n_sig = int(n_base / r_edge);
  498.       #ifdef _DEBUG
  499.         *DebugFile << " case 1:n"
  500.                   << "   new samp intvl = " << t_sig
  501.                   << "   new block size = " << n_sig
  502.                   << endl;
  503.       #endif
  505.       break;
  506.     case 2:   //backward, sampling interval defined
  507.       // compute resampling rate : R = Tn / Tb
  508.       // compute block size : Nn = Nb / R
  510.       r_edge = t_sig/t_base;
  511.       n_sig = int(n_base/r_edge);
  512.       #ifdef _DEBUG
  513.         *DebugFile << " case 2:n"
  514.                   << "   new block size = " << n_sig
  515.                   << "   new resamp rate = " << r_edge
  516.                   << endl;
  517.       #endif
  518.       break;
  519.     case 3:   //backward, resampling rate and 
  520.               //sampling interval are defined
  521.       // compute block size : Nn = Nb / R
  522.       // check sampling interval : (Tn == R * Tb)?
  524.       n_sig = int(n_base / r_edge);
  525.       if( t_sig != (r_edge * t_base))
  526.         { // error
  527.         #ifdef _DEBUG
  528.           *DebugFile << " error in case 3" << endl;
  529.         #endif
  530.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  531.         ErrorStream << "resampling rate in instance '"
  532.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  533.                       << "' of model '" 
  534.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  535.                       << "'n is inconsistent input/output sampling intervals" << endl;
  536.         exit(0);
  537.         }
  538.       #ifdef _DEBUG
  539.         *DebugFile << " case 3:n"
  540.                   << "   new block size = " << n_sig
  541.                   << endl;
  542.       #endif
  543.       break;
  544.     case 4:   //backward, block size is defined
  545.       // compute resampling rate : R = Nb / Nn
  546.       // compute sampling interval : Tn = R * Tb
  548.       r_edge = float(n_base) / n_sig;
  549.       t_sig = r_edge * t_base;
  550.       #ifdef _DEBUG
  551.         *DebugFile << " case 4:n"
  552.                   << "   new samp interval = " << t_sig
  553.                   << "   new resamp rate = " << r_edge
  554.                   << endl;
  555.       #endif
  556.       break;
  557.     case 5:   //backward, resampling rate and block size defined
  558.       // compute sampling interval : Tn = R * Tb
  559.       // check block size : (Nb == R * Nn)?
  561.       t_sig = r_edge * t_base;
  562.       if( n_base != (r_edge * n_sig))
  563.         { // error
  564.         #ifdef _DEBUG
  565.           *DebugFile << " error in case 5" << endl;
  566.         #endif
  567.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  568.         ErrorStream << "resampling rate in instance '"
  569.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  570.                       << "' of model '" 
  571.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  572.                       << "'n is inconsistent with input/output block sizes" << endl;
  573.         exit(0);
  574.         }
  575.       #ifdef _DEBUG
  576.         *DebugFile << " case 5:n"
  577.                   << "   new samp interval = " << t_sig
  578.                   << endl;
  579.       #endif
  580.       break;
  581.     case 6:   //backward, sampling interval and block size defined
  582.       // compute resampling rate : R = Tn / Tb
  583.       // check block size : (Nb == R * Nn)?
  585.       r_edge = t_sig / t_base;
  586.       if( n_base != (r_edge * n_sig))
  587.         { //error
  588.         #ifdef _DEBUG
  589.           *DebugFile << " error in case 6" << endl;
  590.         #endif
  591.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  592.         ErrorStream << "resampling rate in instance '"
  593.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  594.                       << "' of model '" 
  595.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  596.                       << "'n is inconsistent with input/output block sizes" << endl;
  597.         exit(0);
  598.         }
  599.       #ifdef _DEBUG
  600.         *DebugFile << " case 6:n"
  601.                   << "   new resamp rate = " << r_edge
  602.                   << endl;
  603.       #endif
  604.       break;
  605.     case 7:   //backward, all parameters are defined
  606.       // check : (Tn == R * Tb)?
  607.       // check : (Nb == R * Nn)?
  609.       if( t_sig != (r_edge * t_base))
  610.         { // error
  611.         #ifdef _DEBUG
  612.           *DebugFile << " error in case 7a" << endl;
  613.         #endif
  614.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  615.         ErrorStream << "resampling rate in instance '"
  616.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  617.                       << "' of model '" 
  618.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  619.                       << "'n is inconsistent with input/output sampling intervals" << endl;
  620.         exit(0);
  621.         }
  622.       if( n_base != ceil(r_edge * n_sig))
  623.         { // error
  624.         #ifdef _DEBUG
  625.           *DebugFile << " error in case 7b" << endl;
  626.         #endif
  627.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  628.         ErrorStream << "resampling rate in instance '"
  629.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  630.                       << "' of model '" 
  631.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  632.                       << "'n is inconsistent with input/output block sizes" << endl;
  633.         ErrorStream << "n_base = " << n_base << endl;
  634.         ErrorStream << "r_edge = " << r_edge << endl;
  635.         ErrorStream << "n_sig = " << n_sig << endl;
  636.         exit(0);
  637.         }
  638.       #ifdef _DEBUG
  639.         *DebugFile << " case 7 - all parameters defined" << endl;
  640.       #endif
  641.       break;
  642.     case 8:   //forward, nothing defined
  643.       // no action possible
  644.       #ifdef _DEBUG
  645.         *DebugFile << " case 8 - no action possible" << endl;
  646.       #endif
  647.       return;
  648.     case 9:   //forward, resampling rate is defined
  649.       // compute sampling interval : Tn = Tb / R
  650.       // compute block size : Nn = R * Nb
  652.       t_sig = t_base / r_edge;
  653.       n_sig = int(ceil(r_edge * n_base));
  654.       #ifdef _DEBUG
  655.         *DebugFile << " case 9:n"
  656.                   << "   new samp interval = " << t_sig
  657.                   << "   new block size = " << n_sig
  658.                   << endl;
  659.       #endif
  660.       break;
  661.     case 10:    //forward, sampling interval defined
  662.       // compute resampling rate : R = Tb / Tn
  663.       // compute block size : Nn = R * Nb
  665.       r_edge = t_base / t_sig;
  666.       n_sig = int(r_edge * n_base);
  667.       #ifdef _DEBUG
  668.         *DebugFile << " case 10:n"
  669.                   << "   new block size = " << n_sig
  670.                   << "   new resamp rate = " << r_edge
  671.                   << endl;
  672.       #endif
  673.       break;
  674.     case 11:    //forward, resampling rate and
  675.                 // sampling interval defined
  676.       // compute block size : Nn = R * Nb
  677.       // check sampling interval : (Tb == R * Tn)?
  679.       n_sig = int(r_edge * n_base);
  680.       if(t_base != r_edge * t_sig)
  681.         { // error
  682.         #ifdef _DEBUG
  683.           *DebugFile << " error in case 11" << endl;
  684.         #endif
  685.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  686.         ErrorStream << "resampling rate in instance '"
  687.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  688.                       << "' of model '" 
  689.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  690.                       << "'n is inconsistent with input/output sampling intervals" << endl;
  691.         exit(0);
  692.         }
  693.       #ifdef _DEBUG
  694.         *DebugFile << " case 11:n"
  695.                   << "   new samp interval = " << t_sig
  696.                   << "   new block size = " << n_sig
  697.                   << endl;
  698.       #endif
  699.       break;
  700.     case 12:    //forward, block size is defined
  701.       // compute resampling rate : R = Nn / Nb
  702.       // compute sampling interval : Tn = Tb / R
  704.       r_edge = n_sig / n_base;
  705.       t_sig = t_base / r_edge;
  706.       #ifdef _DEBUG
  707.         *DebugFile << " case 12:n"
  708.                   << "   new samp interval = " << t_sig
  709.                   << "   new resamp rate = " << r_edge
  710.                   << endl;
  711.       #endif
  712.       break;
  713.     case 13:    //forward, resampling rate and
  714.                 // block size are defined
  715.       // compute sampling interval : Tn = Tb / R
  716.       // check block size : (Nn == R * Nb)?
  718.       t_sig = t_base / r_edge;
  719.       if(n_sig != (r_edge * n_base))
  720.         { // error
  721.         #ifdef _DEBUG
  722.           *DebugFile << " error in case 13" << endl;
  723.         #endif
  724.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  725.         ErrorStream << "resampling rate in instance '"
  726.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  727.                       << "' of model '" 
  728.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  729.                       << "'n is inconsistent with input/output block sizes" << endl;
  730.         exit(0);
  731.         }
  732.       #ifdef _DEBUG
  733.         *DebugFile << " case 13:n"
  734.                   << "   new samp interval = " << t_sig
  735.                   << endl;
  736.       #endif
  737.       break;
  738.     case 14:    //forward, sampling interval and
  739.                 // block size are defined
  740.       // compute resampling rate : R = Tb / Tn
  741.       // check block size : (Nn == R * Nb)?
  743.       r_edge = t_base / t_sig;
  744.       if( n_sig != (r_edge * n_base))
  745.         { // error
  746.         #ifdef _DEBUG
  747.           *DebugFile << " error in case 14" << endl;
  748.           *DebugFile << "   ( n_sig != (r_edge * n_base))" << endl;
  749.           *DebugFile << "    n_sig = " << n_sig << endl;
  750.           *DebugFile << "    r_edge = " << r_edge << endl;
  751.           *DebugFile << "    n_base = " << n_base << endl;
  752.         #endif
  753.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  754.         ErrorStream << "resampling rate in instance '"
  755.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  756.                       << "' of model '" 
  757.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  758.                       << "'n is inconsistent with input/output block sizes" << endl;
  759.         exit(0);
  760.         }
  761.       #ifdef _DEBUG
  762.         *DebugFile << " case 14:n"
  763.                   << "   new resamp rate = " << r_edge
  764.                   << endl;
  765.       #endif
  766.       break;
  767.     case 15:    //forward, all parameters are defined
  768.       // check : (Tb == R * Tn)?
  769.       // check : (Nn == R * Nb)?
  771.       if( t_base != (r_edge * t_sig))
  772.         { // error
  773.         #ifdef _DEBUG
  774.           *DebugFile << " error in case 15a" << endl;
  775.         #endif
  776.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  777.         ErrorStream << "resampling rate in instance '"
  778.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  779.                       << "' of model '" 
  780.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  781.                       << "'n is inconsistent with input/output sampling intervals" << endl;
  782.         exit(0);
  783.         }
  784.       if( n_sig != (r_edge * n_base))
  785.         { // error
  786.         #ifdef _DEBUG
  787.           *DebugFile << " error in case 15b" << endl;
  788.         #endif
  789.         ErrorStream << "nn********* FATAL ERROR ************" << endl;
  790.         ErrorStream << "resampling rate in instance '"
  791.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetInstanceName() 
  792.                       << "' of model '" 
  793.                       << ((Sdg_Edge_Descr->at(edge_num))->model_id)->GetModelName() 
  794.                       << "'n is inconsistent with input/output block sizes" << endl;
  795.         ErrorStream << "n_base = " << n_base << endl;
  796.         ErrorStream << "r_edge = " << r_edge << endl;
  797.         ErrorStream << "n_sig = " << n_sig << endl;
  798.         exit(0);
  799.         }
  800.       #ifdef _DEBUG
  801.         *DebugFile << " case 15 - all parameters defined" << endl;
  802.       #endif
  803.       break;
  804.     case 16: //backward, constant interval,  nothing defined
  805.       // no action possible
  806.       #ifdef _DEBUG
  807.         *DebugFile << "  case 16 - no action possible" << endl;
  808.       #endif
  809.       return;
  810.     case 17:  //backward, constant interval, resampling rate is defined
  811.       // constant sampling rate : Tn = Tb
  812.       // compute block size : Nn = ceil(Nb / R)
  814.       t_sig = t_base;
  815.       n_sig = int(ceil(n_base / r_edge));
  816.       #ifdef _DEBUG
  817.         *DebugFile << " case 17:n"
  818.                   << "   new block size = " << n_sig
  819.                   << endl;
  820.       #endif
  822.       break;
  823.     case 18:   //backward, constant interval, sampling interval defined
  824.       // compute resampling rate : R = Tn / Tb
  825.       //      this is not really constant interval if Tn != Tb
  826.       // compute block size : Nn = ceil(Nb / R)
  828.       r_edge = t_sig/t_base;
  829.       n_sig = int(ceil(n_base/r_edge));
  830.       #ifdef _DEBUG
  831.         *DebugFile << " case 18:n"
  832.                   << "   new block size = " << n_sig
  833.                   << "   new resamp rate = " << r_edge
  834.                   << endl;
  835.       #endif
  836.       break;
  837.     case 19:   //backward, constant interval, resampling rate and 
  838.               //sampling interval are defined
  839.       // compute block size : Nn = ceil(Nb / R)
  841.       n_sig = int(ceil(n_base / r_edge));
  842.       #ifdef _DEBUG
  843.         *DebugFile << " case 19:n"
  844.                   << "   new block size = " << n_sig
  845.                   << endl;
  846.       #endif
  847.       break;
  848.     case 20:   //backward, constant interval, block size is defined
  849.       // compute resampling rate : R = Nb / Nn
  850.       // constant sampling interval : Tn = Tb
  852.       r_edge = float(n_base) / n_sig;
  853.       t_sig = t_base;
  854.       #ifdef _DEBUG
  855.         *DebugFile << " case 20:n"
  856.                   << "   new samp interval = " << t_sig
  857.                   << "   new resamp rate = " << r_edge
  858.                   << endl;
  859.       #endif
  860.       break;
  861.     case 21:   //backward, constant interval, resampling rate and block size defined
  862.       // constant sampling interval : Tn = Tb
  864.       t_sig = t_base;
  865.       #ifdef _DEBUG
  866.         *DebugFile << " case 21:n"
  867.                   << "   new samp interval = " << t_sig
  868.                   << endl;
  869.       #endif
  870.       break;
  871.     case 22:   //backward, constant interval, sampling interval and block size defined
  872.       // compute resampling rate : R = Nb / Nn
  874.       r_edge = n_base / n_sig;
  875.       #ifdef _DEBUG
  876.         *DebugFile << " case 22:n"
  877.                   << "   new resamp rate = " << r_edge
  878.                   << endl;
  879.       #endif
  880.       break;
  881.     case 23:   //backward, constant interval, all parameters are defined
  882.       #ifdef _DEBUG
  883.         *DebugFile << " case 23 - all parameters defined" << endl;
  884.       #endif
  885.       break;
  886.     case 24:   //forward, constant interval, nothing defined
  887.       // no action possible
  888.       #ifdef _DEBUG
  889.         *DebugFile << " case 24 - no action possible" << endl;
  890.       #endif
  891.       return;
  892.     case 25:   //forward, constant interval, resampling rate is defined
  893.       // constant sampling interval : Tn = Tb
  894.       // compute block size : Nn = ceil(R * Nb)
  896.       t_sig = t_base;
  897.       n_sig = int(ceil(r_edge * n_base));
  898.       #ifdef _DEBUG
  899.         *DebugFile << " case 25:n"
  900.                   << "   new samp interval = " << t_sig
  901.                   << "   new block size = " << n_sig
  902.                   << endl;
  903.       #endif
  904.       break;
  905.     case 26:    //forward, constant interval, sampling interval defined
  906.       // compute resampling rate : R = Tb / Tn
  907.       // compute block size : Nn = ceil(R * Nb)
  909.       r_edge = t_base / t_sig;
  910.       n_sig = int(ceil(r_edge * n_base));
  911.       #ifdef _DEBUG
  912.         *DebugFile << " case 26:n"
  913.                   << "   new block size = " << n_sig
  914.                   << "   new resamp rate = " << r_edge
  915.                   << endl;
  916.       #endif
  917.       break;
  918.     case 27:    //forward, constant interval, resampling rate and
  919.                 // sampling interval defined
  920.       // compute block size : Nn = ceil(R * Nb)
  922.       n_sig = int(ceil(r_edge * n_base));
  923.       #ifdef _DEBUG
  924.         *DebugFile << " case 27:n"
  925.                   << "   new block size = " << n_sig
  926.                   << endl;
  927.       #endif
  928.       break;
  929.     case 28:    //forward, constant interval, block size is defined
  930.       // compute resampling rate : R = Nn / Nb
  931.       // constant sampling interval : Tn = Tb
  933.       r_edge = n_sig / n_base;
  934.       t_sig = t_base;
  935.       #ifdef _DEBUG
  936.         *DebugFile << " case 28:n"
  937.                   << "   new samp interval = " << t_sig
  938.                   << "   new resamp rate = " << r_edge
  939.                   << endl;
  940.       #endif
  941.       break;
  942.     case 29:    //forward, constant interval, resampling rate and
  943.                 // block size are defined
  944.       // constant sampling interval : Tn = Tb
  945.       // check block size : (Nn == R * Nb)?
  947.       t_sig = t_base;
  948.       #ifdef _DEBUG
  949.         *DebugFile << " case 29:n"
  950.                   << "   new samp interval = " << t_sig
  951.                   << endl;
  952.       #endif
  953.       break;
  954.     case 30:    //forward, constant interval, sampling interval and
  955.                 // block size are defined
  956.       // compute resampling rate : R = Nn / Nb
  958.       r_edge = n_sig / n_base;
  959.       #ifdef _DEBUG
  960.         *DebugFile << " case 30:n"
  961.                   << "   new resamp rate = " << r_edge
  962.                   << endl;
  963.       #endif
  964.       break;
  965.     case 31:    //forward, constant interval, all parameters are defined
  967.       #ifdef _DEBUG
  968.         *DebugFile << " case 31 - all parameters defined" << endl;
  969.       #endif
  970.       break;
  971.     default:
  972.       *DebugFile << "Fatal Error in Propagate()" << endl;
  973.       exit(0);
  974.       
  975.     } // end of switch on config_key
  976.       // check for integer block size
  978. //  if( (int(r_edge * n_sig)) != n_base )
  979. //    { // errror
  980. //    *DebugFile << "base block size = " << n_base << endl;
  981. //    *DebugFile << "resampling rate = " << r_edge << endl;
  982. //    *DebugFile << "calculated block size = " << (float(n_base)/r_edge) 
  983. //              << endl;
  984. //    }
  985.   (Sdg_Vert_Descr->at(sig_num))->samp_intvl = t_sig;
  986.   (Sdg_Vert_Descr->at(sig_num))->block_size = n_sig;
  987.   (Sdg_Edge_Descr->at(edge_num))->resamp_rate = r_edge;
  988.   return;
  989. }
  990. //====================================================================
  991. //  Set the propagated values of signal parameters within each
  992. //  individual signal object
  993. void SystemGraph::DistributeSignalParms(void)
  994. {
  995.   int sig_num, block_size;
  996.   GenericSignal *sig_id;
  997.   double samp_intvl;
  998.   int num_nodes = Sig_Dep_Graph->GetNumVerts();
  1000.   for( sig_num = 0; sig_num < num_nodes; sig_num++)
  1001.     {
  1002.     if( ((Sdg_Vert_Descr->at(sig_num))->kind_of_signal) != SK_REGULAR_SIGNAL) continue;
  1003.       //
  1004.       // else
  1005.       sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1007.       block_size = (Sdg_Vert_Descr->at(sig_num))->block_size;
  1008.       sig_id->SetBlockSize( block_size );
  1010.       //buffer_size = (Sdg_Vert_Descr->at(sig_num))->buffer_size;
  1011.       //sig_id->SetBufferSize( buffer_size );
  1013.       samp_intvl = (Sdg_Vert_Descr->at(sig_num))->samp_intvl;
  1014.       sig_id->SetSampIntvl( samp_intvl );
  1015.     }
  1016.   return;
  1017. }
  1018. //====================================================================
  1019. void SystemGraph::AllocateStorageBuffers(void)
  1020. {
  1021.   int sig_num;
  1022.   GenericSignal *sig_id;
  1023.   int num_nodes = Sig_Dep_Graph->GetNumVerts();
  1025.   for(sig_num = 0; sig_num < num_nodes; sig_num++)
  1026.     {
  1027.     //------------------------------------------------
  1028.     // don't consider dummy signals
  1030.     if( ((Sdg_Vert_Descr->at(sig_num))->kind_of_signal) != SK_REGULAR_SIGNAL) continue;
  1031.       //
  1032.       // else
  1033.       sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1034.       sig_id->AllocateSignalBuffer();
  1035.     }
  1036. }
  1037. //============================================================
  1038. void SystemGraph::InitializeReadPtrs(void)
  1039. {
  1040.   GenericSignal* sig_id;
  1041.   int num_sigs = Sig_Dep_Graph->GetNumVerts();
  1042.   for(int sig_num = 0; sig_num < num_sigs; sig_num++)
  1043.     {
  1044.     if( ((Sdg_Vert_Descr->at(sig_num))->kind_of_signal) != SK_REGULAR_SIGNAL ) continue;
  1045.         //
  1046.         // else
  1047.         sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1048.         sig_id->InitializeReadPtrs();
  1049.     }
  1050. }
  1051. //===================================================================
  1052. void SystemGraph::InitializeModels(void)
  1053. {
  1054.   PracSimModel *model_id;
  1055.   *DebugFile << "In SystemGraph::InitializeModels()" << endl;
  1056.   EnclaveNumber = 0;
  1057.   int num_edges = Sig_Dep_Graph->GetNumEdges();
  1058.   for( int edge_num = 0; edge_num < num_edges; edge_num ++)
  1059.     {
  1060.     model_id = (Sdg_Edge_Descr->at(edge_num))->model_id;
  1061.     // note: this might pick up models more than once!
  1063.     model_id->Initialize();
  1064.     *DebugFile << "Initialized model " << model_id->GetModelName() << ":"
  1065.                << model_id->GetInstanceName() << endl;
  1066.     }
  1067. }
  1068. //==============================================================================
  1069. void SystemGraph::DeleteModels(void)
  1070. {
  1071.   PracSimModel *model_id;
  1072.   *DebugFile << "In SystemGraph::InitializeModels()" << endl;
  1073.   int num_edges = Sig_Dep_Graph->GetNumEdges();
  1074.   for(int model_num = 0; model_num < Num_Sys_Lev_Models; model_num++)
  1075.     {
  1076.     model_id = Syst_Lev_Models->at(model_num);
  1078.     #ifdef _DEBUG
  1079.     *DebugFile << "Deleting model " << model_id->GetModelName() << ":"
  1080.                << model_id->GetInstanceName() << endl;
  1081.     #endif
  1082.     delete model_id;
  1083.     }
  1084. }
  1085. //==============================================================================
  1086. void SystemGraph::RunSimulation(void)
  1087. {
  1088.   int model_num;
  1089.   int model_exec_status;
  1091.   EnclaveNumber = 0;
  1092.   EnclaveOffset[0] = 0;
  1093.   for(model_num = 0; model_num < Num_Sys_Lev_Models; model_num++)
  1094.     {
  1095.     ActiveModel = Syst_Lev_Models->at(model_num);
  1096.     #ifdef _DEBUG
  1097.       *DebugFile << "ASG launching Execute for " << ActiveModel->GetModelName() 
  1098.                   << ":" << ActiveModel->GetInstanceName() << endl;
  1099.     #endif
  1100.     model_exec_status = ActiveModel->Execute();
  1101.     if( model_exec_status == _MES_AOK ) continue;
  1102.       // else take action appropriate for returned status code
  1103.       switch (model_exec_status)
  1104.         {
  1105.         case _MES_RESTART:
  1106.           *DebugFile << "Restarting at top of Enclave 0" << endl;
  1107.           goto signal_update;
  1108.         default: //error
  1109.           cout << "unrecognized status returned by model execute method" << endl;
  1110.           exit(1);
  1111.         } // end of switch on model_exec_status
  1112.     }
  1113.   SigPlot.CollectData();
  1114.   //----------------------------------------------------------------------
  1115.   // Update signals
  1116. signal_update:
  1117.   int num_nodes = Sig_Dep_Graph->GetNumVerts();
  1118.   GenericSignal *sig_id;
  1120.   for( int sig_num = 0; sig_num < num_nodes; sig_num++)
  1121.     {
  1122.     if( ((Sdg_Vert_Descr->at(sig_num))->kind_of_signal) != SK_REGULAR_SIGNAL) continue;
  1123.       //
  1124.       // else
  1125.       sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1126.       sig_id->PassUpdate();
  1127.     }
  1128. }
  1129. //=========================================================================
  1130. void SystemGraph::RegisterModel(PracSimModel* model)
  1131. {
  1132.   if(model->GetNestDepth() == 1)
  1133.     {
  1134.     // add model to system level list
  1135.     Syst_Lev_Models->push_back(model);
  1136.     Num_Sys_Lev_Models++;
  1137.     }
  1138. }
  1139. //================================================================
  1140. void SystemGraph::AllocatePlotPointers(void)
  1141. {
  1142.   GenericSignal* sig_id;
  1143.   int num_sigs = Sig_Dep_Graph->GetNumVerts();
  1144.   for(int sig_num = 0; sig_num < num_sigs; sig_num++)
  1145.     {
  1146.     sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1147.     sig_id->SetupPlotSignal();
  1148.     }
  1149. }
  1150. //==========================================================
  1151. GenericSignal* SystemGraph::GetSignalId( char* sig_name)
  1152. {
  1153.   GenericSignal *sig_id;
  1154.   int num_sigs = Sig_Dep_Graph->GetNumVerts();
  1155.   for(int sig_num = 0; sig_num < num_sigs; sig_num++)
  1156.     {
  1157.     sig_id = (Sdg_Vert_Descr->at(sig_num))->signal_id;
  1158.     if( strcmp( sig_name, sig_id->GetName())) continue;
  1159.       return(sig_id);      
  1160.     }
  1161.   return(NULL);
  1162. }