|
5 | 5 | // * that implements a variety of quantum chemistry methods for strongly
|
6 | 6 | // * correlated electrons.
|
7 | 7 | // *
|
8 |
| -// * Copyright (c) 2012-2024 by its authors (see COPYING, COPYING.LESSER, AUTHORS). |
| 8 | +// * Copyright (c) 2012-2025 by its authors (see COPYING, COPYING.LESSER, AUTHORS). |
9 | 9 | // *
|
10 | 10 | // * The copyrights for code used from other parties are included in
|
11 | 11 | // * the corresponding files.
|
|
79 | 79 | // if (not options_->is_none("SCF_TYPE")) {
|
80 | 80 | // auto scf_type = options_->get_str("SCF_TYPE");
|
81 | 81 | // if (scf_type == "PK" or scf_type == "OUT_OF_CORE") {
|
82 |
| -// outfile->Printf("\n\n Please change SCF_TYPE to DIRECT for conventional integrals."); |
83 |
| -// throw psi::PSIEXCEPTION("Please change SCF_TYPE to DIRECT."); |
| 82 | +// outfile->Printf("\n\n Please change SCF_TYPE to DIRECT for conventional |
| 83 | +// integrals."); throw psi::PSIEXCEPTION("Please change SCF_TYPE to DIRECT."); |
84 | 84 | // }
|
85 | 85 | // }
|
86 | 86 |
|
|
249 | 249 | // Ediff = E_mcscf_ - E_mcscf_old;
|
250 | 250 | // if (iter > 1 && std::fabs(Ediff) < econv && g_norm < gconv) {
|
251 | 251 |
|
252 |
| -// outfile->Printf("\n %4d %10.12f %10.12f %10.12f %10.6f s", iter, g_norm, Ediff, |
| 252 | +// outfile->Printf("\n %4d %10.12f %10.12f %10.12f %10.6f s", iter, g_norm, |
| 253 | +// Ediff, |
253 | 254 | // E_mcscf_, mcscf_total_iter.get());
|
254 | 255 |
|
255 |
| -// outfile->Printf("\n\n A miracle has come to pass: CASSCF iterations have converged."); |
256 |
| -// break; |
| 256 | +// outfile->Printf("\n\n A miracle has come to pass: CASSCF iterations have |
| 257 | +// converged."); break; |
257 | 258 | // }
|
258 | 259 |
|
259 | 260 | // Sstep = orbital_optimizer.approx_solve();
|
|
494 | 495 |
|
495 | 496 | // std::shared_ptr<ActiveSpaceIntegrals> CASSCF::get_ci_integrals() {
|
496 | 497 | // std::vector<int> actv_sym = mo_space_info_->symmetry("ACTIVE");
|
497 |
| -// auto fci_ints = std::make_shared<ActiveSpaceIntegrals>(ints_, actv_mos_, actv_sym, rdocc_mos_); |
| 498 | +// auto fci_ints = std::make_shared<ActiveSpaceIntegrals>(ints_, actv_mos_, actv_sym, |
| 499 | +// rdocc_mos_); |
498 | 500 |
|
499 | 501 | // if (!(options_->get_bool("RESTRICTED_DOCC_JK"))) {
|
500 | 502 | // fci_ints->set_active_integrals_and_restricted_docc();
|
|
509 | 511 | // size_t nactv3 = nactv2 * nactv_;
|
510 | 512 |
|
511 | 513 | // active_ab.iterate([&](const std::vector<size_t>& i, double& value) {
|
512 |
| -// value = tei_paaa_data[actv_mos_[i[0]] * nactv3 + i[1] * nactv2 + i[2] * nactv_ + i[3]]; |
| 514 | +// value = tei_paaa_data[actv_mos_[i[0]] * nactv3 + i[1] * nactv2 + i[2] * nactv_ + |
| 515 | +// i[3]]; |
513 | 516 | // });
|
514 | 517 |
|
515 | 518 | // active_aa.copy(active_ab);
|
516 | 519 | // active_aa("u,v,x,y") -= active_ab("u, v, y, x");
|
517 | 520 |
|
518 | 521 | // fci_ints->set_active_integrals(active_aa, active_ab, active_aa);
|
519 | 522 | // if (mcscf_debug_print_) {
|
520 |
| -// outfile->Printf("\n\n tei_active_aa: %8.8f, tei_active_ab: %8.8f", active_aa.norm(2), |
| 523 | +// outfile->Printf("\n\n tei_active_aa: %8.8f, tei_active_ab: %8.8f", |
| 524 | +// active_aa.norm(2), |
521 | 525 | // active_ab.norm(2));
|
522 | 526 | // }
|
523 | 527 |
|
|
683 | 687 | // // D_{uv}^{restricted} = \sum_{i}^{restricted} C_{ui} * C_{vi}
|
684 | 688 |
|
685 | 689 | // // grab part of Ca for inactive docc
|
686 |
| -// auto Cdocc = std::make_shared<psi::Matrix>("C_INACTIVE", nirrep_, nsopi_, inactive_docc_dim_); |
687 |
| -// for (size_t h = 0; h < nirrep_; h++) { |
| 690 | +// auto Cdocc = std::make_shared<psi::Matrix>("C_INACTIVE", nirrep_, nsopi_, |
| 691 | +// inactive_docc_dim_); for (size_t h = 0; h < nirrep_; h++) { |
688 | 692 | // for (int i = 0; i < inactive_docc_dim_[h]; i++) {
|
689 | 693 | // Cdocc->set_column(h, i, Ca->get_column(h, i));
|
690 | 694 | // }
|
|
0 commit comments