docs: private double cost(DMatrixRMaj param )
      {
          function.compute(param, residuals);
  
          double error = NormOps_DDRM.normF(residuals);
  
          return error*error / (double)residuals.numRows;
      }