MonteCarlo.cpp

   1 // Monte Carlo algorithm for computing the denominator Z in a Boltzmann model
   2 // Copyright © 2009 The University of Chicago
   3 #include "MonteCarlo.h"
   4
   5 #include <cstdlib>
   6 #include <Q3TextStream>
   7
   8 MonteCarlo::MonteCarlo(int size)
   9 {
  10         m_Size = size;
  11         m_Values.resize(size);
  12         m_ReadyForAction = false;
  13 };
  14
  15 MonteCarlo::MonteCarlo(int size, QString MyFirstPhone)
  16 {
  17         m_Size = size;
  18         m_Values.resize(size);
  19         m_ReadyForAction = false;
  20         m_MyFirstPhone = MyFirstPhone;
  21 };
  22
  23 void MonteCarlo::Dump( Q3TextStream* out )
  24 {
  25         QString                 FirstPhone,
  26                                         SecondPhone;
  27
  28         int                             Size;
  29
  30         MonteCarlo*             pConditionalMonteCarlo;
  31
  32
  33         *out << endl <<
  34                 "MonteCarlo dump. Size: " << m_Size;
  35         if (m_ReadyForAction )
  36         {
  37                 *out << " Ready for action (normalized). ";
  38         }
  39         else
  40         {
  41                 *out << " Not ready for action (normalized). ";
  42         }
  43         if (m_ModelType == UNIGRAM) { *out << endl <<"Model type is Unigram.\n"; };
  44         if (m_ModelType == BIGRAM)  { *out << endl <<"Model type is Bigram.\n" ;};
  45
  46
  47         *out << endl << "Top level links: ";
  48         for (int i = 0 ; i < m_Size; i++)
  49         {
  50                 *out << i << "\t" << m_Keys[i] << "\t" << m_Values[i];
  51         }
  52
  53
  54         if (m_ModelType == BIGRAM)
  55         {
  56                 for (int i = 0; i < GetSize(); i++)
  57                 {
  58                         Size = GetSize();
  59
  60                         pConditionalMonteCarlo = m_MyBigrams.find( m_MyFirstPhone );
  61                         *out << endl << "<<" << FirstPhone << ">>";
  62
  63                         for (int j = 0; j < pConditionalMonteCarlo->GetSize();j++)
  64                         {
  65                                 *out << "\n\t" << j << "\t" << m_Keys[j] << "\t" << m_Values[j];
  66                         }
  67                 }
  68         }
  69
  70
  71
  72
  73
  74
  75 }
  76
  77
  78 QString MonteCarlo::ReturnCharacter( )
  79 {
  80         using std::rand;
  81
  82         int i= 0;
  83         double random;
  84
  85
  86         if (m_ReadyForAction == false) Normalize();
  87
  88         random = rand()/(double)0x7fff;
  89
  90         for ( i = 0; i < m_Size; i++)
  91         {
  92 //unused variable:
  93 //              double temp ( m_Values[i]) ;
  94                 if ( m_Values[i] > random ) break;
  95         }
  96
  97         return m_Keys[i];
  98
  99 }
 100 QString MonteCarlo::ReturnCharacter( QString Char )
 101 {
 102         using std::rand;
 103
 104         int i= 0;
 105         double random;
 106
 107         MonteCarlo*             pConditionalMonteCarlo;
 108         pConditionalMonteCarlo = m_MyBigrams.find( Char );
 109
 110         random = rand()/(double)0x7fff;
 111
 112         for ( i = 0; i < m_Size; i++)
 113         {
 114 //unused variable:
 115 //              double temp ( pConditionalMonteCarlo->m_Values[i]) ;
 116                 if      ( pConditionalMonteCarlo->m_Values[i] > random ) break;
 117         }
 118
 119         return pConditionalMonteCarlo->m_Keys[i];
 120
 121 }
 122 void    MonteCarlo::Normalize()
 123 {
 124         double total = 0;
 125         int                     i= 0;
 126         for (i = 0; i <m_Size; i++)
 127         {
 128                 total += m_Values[i];
 129         }
 130
 131         if (total != 1.0)
 132         {
 133                 for (i = 0; i <m_Size; i++)
 134                 {
 135                         m_Values[i] /= total;
 136                 }
 137                 for (i = 1; i <m_Size; i++)
 138                 {
 139                         m_Values[i] += m_Values[i-1];
 140                 }
 141         }
 142         //check to see if collection includes #; if not, do not say it's ready for action.
 143
 144         m_ReadyForAction = true;
 145
 146 }
 147
 148 CParse  MonteCarlo::ReturnString(CParse& String)
 149 {
 150         QString Char, Char2;
 151         QString temp;
 152         String.ClearParse();
 153         String.Append(QChar('#'));
 154
 155         Char = "#";
 156         while (TRUE)
 157         {
 158                 Char2 = ReturnCharacter ( Char );
 159                 if (Char2.length() == 0) continue;
 160                 String.Append( Char2);
 161                 temp = String.Display();
 162                 if (Char2 == "#") break;
 163                 Char = Char2;
 164         }
 165
 166         QString temp2 = String.GetPiece(2).Display();
 167
 168         return String;
 169 }