| 676 | //} |
| 677 | |
| 678 | double norm(VectorXd& b) { |
| 679 | double norm = 0.0; |
| 680 | #pragma omp parallel for reduction(+:norm) |
| 681 | for (int i = 0; i < b.size(); i++) { |
| 682 | norm += b(i) * (i); |
| 683 | } |
| 684 | return norm; |
| 685 | } |
| 686 | |
| 687 | void get_dR_dC(vector<double>& dR, vector<double>& dC, SparseQMatrix<double>& m) { |
| 688 | int size = m.size; |