[Rcpp-commits] r1968 - in pkg/Rcpp: . inst/include/Rcpp/stats

[noreply at r-forge.r-project.org]

Author: romain
Date: 2010-08-10 10:50:10 +0200 (Tue, 10 Aug 2010)
New Revision: 1968

Added:
   pkg/Rcpp/inst/include/Rcpp/stats/cauchy.h
Modified:
   pkg/Rcpp/TODO
   pkg/Rcpp/inst/include/Rcpp/stats/stats.h
Log:
dpq cauchy

Modified: pkg/Rcpp/TODO
===================================================================
--- pkg/Rcpp/TODO	2010-08-10 08:46:37 UTC (rev 1967)
+++ pkg/Rcpp/TODO	2010-08-10 08:50:10 UTC (rev 1968)
@@ -68,7 +68,7 @@
     o	other statistical distribution functions : 
     	nchisq, nbinom, nbinom_mu, nbeta, nf, nt,
     	multinom, 
-    	cauchy, geom, hyper, tukey, wilcox, signrank
+    	geom, hyper, tukey, wilcox, signrank
     
     o	sum : deal with missing values
 

Added: pkg/Rcpp/inst/include/Rcpp/stats/cauchy.h
===================================================================
--- pkg/Rcpp/inst/include/Rcpp/stats/cauchy.h	                        (rev 0)
+++ pkg/Rcpp/inst/include/Rcpp/stats/cauchy.h	2010-08-10 08:50:10 UTC (rev 1968)
@@ -0,0 +1,118 @@
+
+// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 4 -*-
+//
+// auto generated file (from script/stats.R) 
+//
+// cauchy.h: Rcpp R/C++ interface class library -- 
+//
+// Copyright (C) 2010 Douglas Bates, 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/>.
+
+#ifndef Rcpp__stats__cauchy_h
+#define Rcpp__stats__cauchy_h
+
+namespace Rcpp {
+namespace stats {
+namespace impl {
+
+	template <bool NA, typename T>
+	class DCauchy : public Rcpp::VectorBase< REALSXP, NA, DCauchy<NA,T> >{
+	public:
+		typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+	
+		DCauchy( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 , bool log_ = false ) : 
+			vec(vec_), location(location_), scale(scale_) , log(log_) {}
+		
+		inline double operator[]( int i) const {
+			return ::dcauchy( vec[i], location, scale , log );
+		}
+		
+		inline int size() const { return vec.size(); }
+		
+	private:
+		const VEC_TYPE& vec;
+		double location; double scale; 
+		int log;
+	
+	};
+
+	template <bool NA, typename T>
+	class PCauchy : public Rcpp::VectorBase< REALSXP, NA, PCauchy<NA,T> >{
+	public:
+		typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+	
+		PCauchy( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 ,
+			   bool lower_tail = true, bool log_ = false ) : 
+			vec(vec_), location(location_), scale(scale_) , lower(lower_tail), log(log_) {}
+		
+		inline double operator[]( int i) const {
+			return ::pcauchy( vec[i], location, scale, lower, log );
+		}
+		
+		inline int size() const { return vec.size(); }
+		
+	private:
+		const VEC_TYPE& vec;
+		double location; double scale; 
+		int lower, log;
+	
+	};
+
+	template <bool NA, typename T>
+	class QCauchy : public Rcpp::VectorBase< REALSXP, NA, QCauchy<NA,T> >{
+	public:
+		typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+	
+		QCauchy( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 ,
+			   bool lower_tail = true, bool log_ = false ) : 
+			vec(vec_), location(location_), scale(scale_), lower(lower_tail), log(log_) {}
+		
+		inline double operator[]( int i) const {
+			return ::qcauchy( vec[i], location, scale, lower, log );
+		}
+		
+		inline int size() const { return vec.size(); }
+		
+	private:
+		const VEC_TYPE& vec;
+		double location; double scale; 
+		int lower, log;
+	
+	};
+	
+} // impl
+
+template <bool NA, typename T>
+inline impl::DCauchy<NA,T> dcauchy( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool log = false ) {
+	return impl::DCauchy<NA,T>( x, location_, scale_, log ); 
+}
+
+template <bool NA, typename T>
+inline impl::PCauchy<NA,T> pcauchy( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool lower = true, bool log = false ) {
+	return impl::PCauchy<NA,T>( x, location_, scale_, lower, log ); 
+}
+
+template <bool NA, typename T>
+inline impl::QCauchy<NA,T> qcauchy( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool lower = true, bool log = false ) {
+	return impl::QCauchy<NA,T>( x, location_, scale_, lower, log ); 
+}
+	
+}
+}
+
+#endif
+

Modified: pkg/Rcpp/inst/include/Rcpp/stats/stats.h
===================================================================
--- pkg/Rcpp/inst/include/Rcpp/stats/stats.h	2010-08-10 08:46:37 UTC (rev 1967)
+++ pkg/Rcpp/inst/include/Rcpp/stats/stats.h	2010-08-10 08:50:10 UTC (rev 1968)
@@ -34,6 +34,8 @@
 #include <Rcpp/stats/weibull.h>
 #include <Rcpp/stats/logis.h>
 #include <Rcpp/stats/f.h>
+#include <Rcpp/stats/exp.h>
+#include <Rcpp/stats/cauchy.h>
 
 #include <Rcpp/stats/binom.h>
 #include <Rcpp/stats/pois.h>