Rcpp Version 1.0.9
convolve12_cpp.cpp
Go to the documentation of this file.
1 // -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 8 -*-
2 
3 // This is a rewrite of the 'Writing R Extensions' section 5.10.1 example
4 
5 #include <Rcpp.h>
6 
7 RcppExport SEXP convolve12cpp(SEXP a, SEXP b){
8  Rcpp::NumericVector xa(a), xb(b);
9  int n_xa = xa.size(), n_xb = xb.size();
10  Rcpp::NumericVector xab(n_xa + n_xb - 1);
11 
12  typedef Rcpp::NumericVector::iterator vec_iterator ;
13  vec_iterator ia = xa.begin(), ib = xb.begin();
14  vec_iterator iab = xab.begin();
15  for (int i = 0; i < n_xa; i++)
16  for (int j = 0; j < n_xb; j++)
17  iab[i + j] += ia[i] * ib[j];
18 
19  return xab;
20 }
21 
22 #include "loopmacro.h"
24 
#define RcppExport
Definition: RcppCommon.h:140
R_xlen_t size() const
Definition: Vector.h:276
traits::r_vector_iterator< RTYPE, StoragePolicy >::type iterator
Definition: Vector.h:46
iterator begin()
Definition: Vector.h:334
RcppExport SEXP convolve12cpp(SEXP a, SEXP b)
#define LOOPMACRO_CPP(name)
Definition: loopmacro.h:12