void expectedColIn_range() { cout << "- Compute expectedColIn_range() ... "; if(_genColVec.in_range(_rowIndRange, _colIndRange)) { save<double>("Col.in_range", Mat<double>({1})); } else { save<double>("Col.in_range", Mat<double>({0})); } cout << "done." << endl; }
void expectedRowsElemDivide() { for(int n = 0; n < _elemInds.n_elem; n++) { if(!_genColVec.in_range(_elemInds.at(n))) { return; } } if(_genRowVec.n_elem != _elemInds.n_elem) { return; } cout << "- Compute expectedElemElemDivide() ... "; _genColVec.rows(_elemInds) /= _genRowVec; save<double>("Col.rowsElemDivide", _genColVec); cout << "done." << endl; }
void expectedColVecSubvecElemDivide() { if(_elemIndRange.whole) { return; } if (!_genColVec.in_range(_elemIndRange)) { return; } if (_genRowVec.n_cols != _genColVec.n_cols || _genRowVec.n_rows != _elemIndRange.b - _elemIndRange.a + 1) { return; } cout << "- Compute expectedColVecSubvecElemDivide() ... "; _genColVec.subvec(_elemIndRange.a, _elemIndRange.b) /= _genRowVec; save<double>("Col.subvecElemDivide", _genColVec); cout << "done." << endl; }
void expectedColVecRowsElemTimes() { if(_elemIndRange.whole) { return; } if (!_genColVec.in_range(_elemIndRange)) { return; } if(_genRowVec.n_cols != _genColVec.n_cols || _genRowVec.n_rows != _elemIndRange.b - _elemIndRange.a + 1) { return; } cout << "- Compute expectedColVecRowsElemTimes() ... "; _genColVec.rows(_elemIndRange.a, _elemIndRange.b) %= _genRowVec; save<double>("Col.rowsElemTimes", _genColVec); cout << "done." << endl; }