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

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

3G开发

开发平台：

Visual C++

cyclo_partit.cpp：源码内容
  1. //
  2. // primpoly.cpp
  3. //
  5. #include <iostream>
  6. #include <fstream>
  7. #include "cyclo_partit.h"
  8. #include "ipow.h"
  9. //using namespace std;
  10. extern ::ofstream DebugFile;
  13. //=====================================================
  14. CyclotomicPartition::CyclotomicPartition( int prime_field_base,
  15.                                                 int degree )
  16. {
  17.    int first_element;
  18.    int total_count;
  19.    int element_count;
  20.    int i;
  21.    int element;
  22.    int base;
  24.    Specific_Coset = new CyclotomicCoset;
  25.    Modulus = int(ipow(prime_field_base, degree)-1);
  26.    total_count = 0;
  27.    Coset_Start = new int[(Modulus+1)/2];
  28.    Coset_Size = new int[(Modulus+1)/2];
  29.    Element_Vector = new int[Modulus];
  30.    bool* used = new bool[Modulus];
  31.    for(i=0; i<Modulus; i++){
  32.       used[i] = false;
  33.    }
  34.    Coset_Start[0] = 0;
  35.    Coset_Size[0] = 1;
  36.    Element_Vector[0] = 0;
  37.    Num_Cosets = 1;
  38.    total_count++;
  40.    first_element = -1;
  41.    do{
  42.       do{
  43.          first_element += 2;
  44.       }while(used[first_element]);
  45.       element = first_element;
  46.       element_count = 0;
  48.       Coset_Start[Num_Cosets] = Coset_Start[Num_Cosets-1] + Coset_Size[Num_Cosets-1];
  49.       base = Coset_Start[Num_Cosets];
  51.       do{
  52.          Element_Vector[base + element_count] = element;
  53.          used[element] = true;
  54.          element_count++;
  55.          total_count++;
  56.          element = (prime_field_base * element)%Modulus;
  57.       } while(element != first_element);
  58.       Coset_Size[Num_Cosets] = element_count;
  59.       Num_Cosets++;
  60.    }while(total_count <(Modulus-1));
  61. }
  62. //=======================================================
  63. int CyclotomicPartition::GetNumCosets(void)
  64. {
  65.    return(Num_Cosets);
  66. }
  67. //=======================================================
  68. CyclotomicCoset *CyclotomicPartition::GetSpecificCoset(int coset_index)
  69. {
  70.    Specific_Coset->start = &(Element_Vector[Coset_Start[coset_index]]);
  71.    Specific_Coset->size = Coset_Size[coset_index];
  72.    return(Specific_Coset);
  73. }
  74. //===========================================================
  75. int CyclotomicPartition::GetCosetIndex(int element)
  76. {
  77.    int i, j, index;
  78.    i=-1;
  79.    do{
  80.       i++;
  81.    }while( (i<(Modulus-1)) && (element != Element_Vector[i]));
  82.    for(j=0; j<Num_Cosets; j++){
  83.       if( (i>=Coset_Start[j]) && (i<(Coset_Start[j]+Coset_Size[j]))){
  84.          index = j;
  85.          break;
  86.       }
  87.    }
  88.    return(index);
  89. }