Changeset 6141 for trunk/STAT/ST_ml.cxx
- Timestamp:
- 14/08/09 16:29:27 (3 years ago)
- Files:
-
- 1 modified
-
trunk/STAT/ST_ml.cxx (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/STAT/ST_ml.cxx
r6010 r6141 307 307 if (distr > maxm) 308 308 distr = maxm; 309 st_ml->rate_matri zes[distl].mult(lb, dest);310 st_ml->rate_matri zes[distr].mult(rb, &hbv);309 st_ml->rate_matrices[distl].mult(lb, dest); 310 st_ml->rate_matrices[distr].mult(rb, &hbv); 311 311 dest->mult(&hbv); 312 312 dest->check_overflow(); … … 337 337 if (distl > maxm) 338 338 distl = maxm; 339 st_ml->rate_matri zes[distl].mult(lefts, out);339 st_ml->rate_matrices[distl].mult(lefts, out); 340 340 341 341 // correct frequencies … … 372 372 delete[]base_frequencies; 373 373 delete[]inv_base_frequencies; 374 delete[]rate_matri zes;374 delete[]rate_matrices; 375 375 if (!awt_csp) { 376 376 delete rates; … … 476 476 } 477 477 478 void ST_ML::create_matri zes(double max_disti, int nmatrizes) {478 void ST_ML::create_matrices(double max_disti, int nmatrices) { 479 479 max_dist = max_disti; 480 if (rate_matri zes)481 delete rate_matri zes;482 rate_matri zes = new ST_rate_matrix[nmatrizes];483 max_matr = nmatri zes;480 if (rate_matrices) 481 delete rate_matrices; 482 rate_matrices = new ST_rate_matrix[nmatrices]; 483 max_matr = nmatrices; 484 484 step_size = max_dist / max_matr; 485 485 int i; 486 486 for (i = 0; i < max_matr; i++) { 487 rate_matri zes[i].set((i + 1) * step_size, 0); // ttratio[i]487 rate_matrices[i].set((i + 1) * step_size, 0); // ttratio[i] 488 488 } 489 489 } … … 611 611 tree_root->tree->load_sequences_rek(alignment_name, GB_TRUE, GB_TRUE); 612 612 613 /* create matri zes */614 create_matri zes(2.0, 1000);613 /* create matrices */ 614 create_matrices(2.0, 1000); 615 615 616 616 ST_sequence_ml::tmp_out = new ST_base_vector[alignment_len];
