void expectedMatIn_range() {
        cout << "- Compute expectedMatIn_range() ... ";

        if(_genMat.in_range(_elemInd)) {
          save<double>("Mat.in_range", Mat<double>({1}));
        } else {
          save<double>("Mat.in_range", Mat<double>({0}));
        }

        cout << "done." << endl;
      }
      void expectedMatSubmatDivide() {
        if(!_genMat.in_range(_rowInd, _colInd, size(_matSize.first, _matSize.second))) {
          return;
        }

        _genMat.submat(_rowInd, _colInd, size(_matSize.first, _matSize.second)) /= _genDouble;

        cout << "- Compute expectedMatSubmatDivide() ... ";
        save<double>("Mat.submatDivide", _genMat);
        cout << "done." << endl;
      }
      void expectedRowsDivide() {
        for(int n = 0; n < _elemInds.n_elem; n++) {
          if(!_genColVec.in_range(_elemInds.at(n))) {
            return;
          }
        }

        cout << "- Compute expectedRowsDivide() ... ";

        _genColVec.rows(_elemInds) /= _genDouble;

        save<double>("Col.rowsDivide", _genColVec);
        cout << "done." << endl;
      }
      void expectedElemTimes() {
        for(int n = 0; n < _elemInds.n_elem; n++) {
          if(!_genColVec.in_range(_elemInds.at(n))) {
            return;
          }
        }

        cout << "- Compute expectedElemTimes() ... ";

        _genColVec.elem(_elemInds) *= _genDouble;

        save<double>("Col.elemTimes", _genColVec);
        cout << "done." << endl;
      }
      void expectedElemDivide() {
        for(int n = 0; n < _elemInds.n_elem; n++) {
          if(!_genMat.in_range(_elemInds.at(n))) {
            return;
          }
        }

        cout << "- Compute expectedElemDivide() ... ";

        _genMat.elem(_elemInds) /= _genDouble;
        save<double>("Mat.elemDivide", _genMat);

        cout << "done." << endl;
      }
      void expectedMatSwapCols() {
        if(_colIndRange.whole) {
          return;
        }

        if(!_genMat.in_range(span::all, _colIndRange)) {
          return;
        }

        cout << "- Compute expectedMatSwapCols() ... ";

        _genMat.swap_cols(_colIndRange.a, _colIndRange.b);
        save<double>("Mat.swap_cols", _genMat);

        cout << "done." << endl;
      }
      void expectedMatSwapRows() {
        if(_rowIndRange.whole) {
          return;
        }

        if(!_genMat.in_range(_rowIndRange, span::all)) {
          return;
        }

        cout << "- Compute expectedMatSwapRows() ... ";

        _genMat.swap_rows(_rowIndRange.a, _rowIndRange.b);
        save<double>("Mat.swap_rows", _genMat);

        cout << "done." << endl;
      }