int main(int argc, char *argv[]) {
-///*
- std::string fname = argv[1];
- std::string fname2 = argv[2];
- PottsHamiltonian H(fname,fname2);
- for (unsigned i=0; i<H.J.size(); i++) {
- for (unsigned j=0; j<H.J[i].size(); j++) {
- std::cout << H.J[i][j] << ' ';
- }
- std::cout << '\n';
- }
- for (unsigned i=0; i<H.aa2state.size(); i++) {
- for (int aa=aa_BEGIN; aa != aa_END; ++aa) {
- printf("%c(%d) ", int2aa.at(aa),H.aa2state[i][aa]);
- //std::cout << H.aa2state[i][aa] << ' ';
- }
- std::cout << '\n';
- }
- //NTVirus v("TTT");
- NTVirus v(H,argv[3]);
- //std::cout << int2aa.at(codon2aa.at(str2codon.at(v.nt_seq))) << '\n';
- std::cout << aa2str.at(v.aa_seq[0]) << '\n';
- std::cout << aaseq2str(v.aa_seq) << '\n';
- for (size_t i=0; i<v.L_aa; i++) {
- std::cout << v.config[i] << ' ';
- }
- std::cout << '\n';
- printf("Nnt=%d Naa=%d E=%f\n",v.L_nt,v.L_aa,v.energy);
- CTLaaSpecies *T = new CTLaaSpecies;
- T->epitopes.push_back(AARecognizer(argv[4],0));
- T->affinity.push_back(1.0);
- T->num_ep = 1;
- if (T->recognized(v) > 0.0)
- std::cout << "got 'im!\n";
- else
- std::cout << "escaped :(\n";
- return 0;
-//*/
-
RunParameters_SS r;
unsigned seed = sim_random_seed();