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

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

3G开发

开发平台：

Visual C++

colvec_pf.cpp：源码内容
  1. //
  2. //  File = colvec_pf.cpp
  3. //
  5. #include "colvec_pf.h"
  6. #include "rowvec_pf.h"
  7. #include "matrix_pf.h"
  8. #include "pfelem.h"
  10. #ifdef _DEBUG
  11.   #include <fstream>
  12.   //#define _VEC_DEBUG 1
  13.   extern ofstream DebugFile;
  14. #endif
  16. //------------------------------------------------
  17. // transpose operator
  18. rowvec_pf& colvec_pf::operator!( void )
  19. {
  20.   //rowvec<T> *rv = new rowvec<T>(pV->orig_indx,pV->length);
  21.   rowvec_pf *rv = new rowvec_pf();
  22.   #ifdef _VEC_DEBUG
  23.     DebugFile << "ncv::op!(): new rowvec at " 
  24.               << (void*)rv << endl;
  25.     DebugFile << "ncv::op!(): hook vrep "
  26.               << (void*)pV << " to rowvec " 
  27.               << (void*)rv << endl;
  28.   #endif
  29.   rv->pV = pV;
  30.   *(rv->pV)=*pV;
  31.   ((rv->pV)->refcnt)++;
  32.   rv->Is_Temp = 1;
  33.   if(Is_Temp)
  34.     {
  35.     #ifdef _VEC_DEBUG
  36.       DebugFile << "ncv::op!(): deleting colvec at " 
  37.                 << (void*)this << endl;
  38.     #endif
  39.     delete this;
  40.     }
  41.   return(*rv);
  42. }
  44. //------------------------------------------
  45. //  constructor for column vector
  46. colvec_pf::colvec_pf( void )
  47.           :vector_pf()
  48.   {
  49.   #ifdef _VEC_DEBUG
  50.     DebugFile << "nshallow ctor for colvec at " << (void*)this 
  51.               << endl;
  52.   #endif
  53.   //pV->is_row_vec = 0;
  54.   }
  56. //------------------------------------------
  57. //  constructor for column vector
  58. colvec_pf::colvec_pf(int origin, int size)
  59.           :vector_pf(origin, size)
  60.   {
  61.   #ifdef _VEC_DEBUG
  62.     DebugFile << "nctor for colvec at " << (void*)this 
  63.               << "  (vrep = " << (void*)pV << ")" << endl;
  64.   #endif
  65.   }
  67. //---------------------------------------------------
  68. //  method to multiply column vector times row vector
  69. matrix_pf& colvec_pf::operator*( rowvec_pf &v2)
  70.   {
  71.   // get origin and length of column vector
  72.   int v1_orig = pV->orig_indx;
  73.   int v1_len = pV->length;
  75.   // get origin and length of row vector
  76.   int v2_orig = v2.pV->orig_indx;
  77.   int v2_len = v2.pV->length;
  79.   // allocate matrix for result
  80.   matrix_pf *m_res = 
  81.                  new matrix_pf(v1_orig,v1_len,v2_orig,v2_len);
  82.   m_res->Is_Temp = 1;
  83.   #ifdef _VEC_DEBUG
  84.     DebugFile << "ncv::op*(rv): new matrix at " 
  85.               << (void*)m_res << endl;
  86.   #endif
  87.   for(int row=0; row<v1_len; row++)
  88.     {
  89.     for(int col=0; col<v2_len; col++)
  90.       {
  91.       ((((m_res->pM)->f[row])->pV)->f[col]) = 
  92.                                 (pV->f[row]) * (v2.pV->f[col]);
  93.       }
  94.     }
  95.   if(v2.Is_Temp)
  96.     {
  97.     #ifdef _VEC_DEBUG
  98.       DebugFile << "ncv::op*(rv): deleting rowvec at " 
  99.                 << (void*)(&v2) << endl;
  100.     #endif
  101.     delete (&v2);
  102.     }
  103.   if(Is_Temp)
  104.     {
  105.     #ifdef _VEC_DEBUG
  106.       DebugFile << "ncv::op*(rv): deleting colvec at " 
  107.                 << (void*)this << endl;
  108.     #endif
  109.     delete this;
  110.     }
  111.   return(*m_res);
  112.   };
  113. //-----------------------------------------
  114. // destructor
  115. colvec_pf::~colvec_pf()
  116.   {
  117.   #ifdef _VEC_DEBUG
  118.     DebugFile << "ndtor for colvec at " << (void*)this << endl;
  119.   #endif
  120.   //~vector();
  121.   };
  123. //----------------------------------------------
  124. // assignment x = y
  125. colvec_pf& colvec_pf::operator=(vector_pf &vec)
  126. {
  127.   vec.pV->refcnt++;
  128.   if(--pV->refcnt == 0)
  129.     {
  130.     #ifdef _VEC_DEBUG
  131.       DebugFile << "ncv::op=(v): deleting vrep at " 
  132.                 << (void*)pV << endl;
  133.     #endif
  134.     delete[] pV->f;
  135.     delete pV;
  136.     }
  137.   #ifdef _VEC_DEBUG
  138.     DebugFile << "ncv::op=(v): hook vrep "
  139.           << (void*)(vec.pV) << " to vector " 
  140.           << (void*)this << endl;
  141.   #endif
  142.   pV = vec.pV;
  143.   if(vec.Is_Temp)
  144.     {
  145.     #ifdef _VEC_DEBUG
  146.       DebugFile << "ncv::op=(v): deleting vector at " 
  147.                 << (void*)(&vec) << endl;
  148.     #endif
  149.     delete &vec;
  150.     }
  151.   return *this;
  152. }
  153. //---------------------------------