| 395 | } |
| 396 | |
| 397 | void AddNZeroIndex(vector<int>& TotalIndex, int jadd) { |
| 398 | auto index_add = No0Index[jadd]; |
| 399 | TotalIndex.insert(TotalIndex.end(), index_add.begin(), index_add.end()); |
| 400 | sort(TotalIndex.begin(), TotalIndex.end()); |
| 401 | TotalIndex.erase(unique(TotalIndex.begin(), TotalIndex.end()), TotalIndex.end()); |
| 402 | } |
| 403 | |
| 404 | void FastGetSubMatrixRight(MatrixXd &subA, const vector<int> &TotalIndex, const vector<int> &Sp, SparseQMatrix<double>& Ae) { |
| 405 |