// RcppMatrixExample.cpp: Rcpp R/C++ interface class library RcppMatrix example

//

// Copyright (C) 2005 - 2006 Dominick Samperi

// Copyright (C) 2008 Dirk Eddelbuettel

// Copyright (C) 2009 - 2023 Dirk Eddelbuettel and Romain Francois

//

// This file is part of Rcpp.

//

// Rcpp is free software: you can redistribute it and/or modify it

// under the terms of the GNU General Public License as published by

// the Free Software Foundation, either version 2 of the License, or

// (at your option) any later version.

//

// Rcpp is distributed in the hope that it will be useful, but

// WITHOUT ANY WARRANTY; without even the implied warranty of

// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the

// GNU General Public License for more details.

//

// You should have received a copy of the GNU General Public License

// along with Rcpp. If not, see <http://www.gnu.org/licenses/>.

#include <RcppClassic.h>

#include <cmath>

/* suncc needs help to disambiguate between sqrt( float ) and sqrt(double) */

inline static double sqrt_double( double x ){ return ::sqrt( x ) ; }

RcppExport SEXP classicRcppMatrixExample(SEXP matrix) {

SEXP rl = R_NilValue; // Use this when there is nothing to be returned.

char *exceptionMesg = NULL;

try {

// Get parameters in params.

RcppMatrix<int> orig(matrix);

int n = orig.rows(), k = orig.cols();

RcppMatrix<double> mat(n, k); // reserve n by k matrix

for (int i=0; i<n; i++) {

for (int j=0; j<k; j++) {

mat(i,j) = sqrt_double(orig(i,j));

}

}

// Build result set to be returned as a list to R.

RcppResultSet rs;

rs.add("result", mat);

rs.add("original", orig);

// Get the list to be returned to R.

rl = rs.getReturnList();

} catch(std::exception& ex) {

exceptionMesg = copyMessageToR(ex.what());

} catch(...) {

exceptionMesg = copyMessageToR("unknown reason");

}

if (exceptionMesg != NULL)

Rcpp::stop("%s", exceptionMesg);

return rl;

}