@@ -32,7 +32,7 @@ SIDD_Base::SIDD_Base()
3232 plasmid_seq = 0 ;
3333 en_cruciforms = 0 ;
3434 profile = 0 ;
35- for (int m = 0 ; m <= MaxInitialWindowSize; m++)
35+ for (int m = 0 ; m <= MaxInitialWindowSize; m++)
3636 promatrix[m] = 0 ;
3737 R = 8.314 /(4.2 *1000.0 ); // gas constant
3838 C = 3.6 ; // stiffness constant
@@ -240,8 +240,8 @@ double SIDD_Base::prepare_sequence(std::string sequence) {
240240void SIDD_Base::prepare_cruciforms () {
241241 string str1 = cr_string;
242242 size_t found1 = 0 ;
243- int IR [2 ] = {0 , 0 };
244- double energy = 0.0 ;
243+ int IR [2 ] = {0 , 0 };
244+ double energy = 0.0 ;
245245 int pos1 = 0 ;
246246 int pos2 = 0 ;
247247 string str2;
@@ -267,16 +267,16 @@ void SIDD_Base::prepare_cruciforms() {
267267 if (found2!=string::npos) {
268268 str3 = str2.substr (pos2,found2-pos2);
269269 pos2 = int (found2)+1 ;
270- if (i == 2 )
271- energy = atof (str3.c_str ());
272- else if (i < 2 )
270+ if (i == 2 )
271+ energy = atof (str3.c_str ());
272+ else if (i < 2 )
273273 IR [i] = atoi (str3.c_str ());
274274 i++;
275275 }
276276 }
277- if (i == 3 && IR [0 ] > 0 && IR [0 ] <= length_seq &&
278- IR [1 ] > 0 && IR [1 ] <= MaxWindowSize)
279- en_cruciforms[IR [0 ]-1 ][IR [1 ]] = energy;
277+ if (i == 3 && IR [0 ] > 0 && IR [0 ] <= length_seq &&
278+ IR [1 ] > 0 && IR [1 ] <= MaxWindowSize)
279+ en_cruciforms[IR [0 ]-1 ][IR [1 ]] = energy;
280280 pos1 = int (found1)+1 ;
281281 }
282282 }
@@ -837,17 +837,17 @@ void SIDD_Base::calc_profile()
837837
838838void SIDD_Base::sum_open ()
839839{
840- double sumn=0.0 ;
841- for (int i = 0 ; i < length_seq; i++)
842- sumn+=profile[i].get_px ();
843- if (results) {
844- if (EnergyType == Z_DNA )
845- cout << " Number of Z-DNA bases = " << sumn << endl;
846- else if (EnergyType == Cruciform)
847- cout << " Number of cruciform bases = " << sumn << endl;
848- else
849- cout << " Number of melted bases = " << sumn << endl;
850- }
840+ double sumn=0.0 ;
841+ for (int i = 0 ; i < length_seq; i++)
842+ sumn+=profile[i].get_px ();
843+ if (results) {
844+ if (EnergyType == Z_DNA )
845+ cout << " Number of Z-DNA bases = " << sumn << endl;
846+ else if (EnergyType == Cruciform)
847+ cout << " Number of cruciform bases = " << sumn << endl;
848+ else
849+ cout << " Number of melted bases = " << sumn << endl;
850+ }
851851}
852852
853853void SIDD_Base::get_column_header ()
@@ -956,7 +956,7 @@ void SIDD_Base::show_parameter()
956956 cout << " Energetics: Melting - nearest neighbor energetics\n " ;
957957 else if (EnergyType == Z_DNA )
958958 cout << " Transition: Z-DNA\n " ;
959- else if (EnergyType == Cruciform)
959+ else if (EnergyType == Cruciform)
960960 cout << " Transition: Cruciform\n " ;
961961 if (MoleculeType == Circular)
962962 cout << " Molecule type: circular DNA\n " ;
0 commit comments