@@ -29,47 +29,47 @@ namespace sugar{
2929template <int RTYPE, bool NA, typename T>
3030class Var : public Lazy < double , Var<RTYPE,NA,T> > {
3131public:
32- typedef typename Rcpp::VectorBase<RTYPE,NA,T> VEC_TYPE ;
32+ typedef typename Rcpp::VectorBase<RTYPE,NA,T> VEC_TYPE ;
3333
34- Var ( const VEC_TYPE& object_ ) : object(object_){}
34+ Var ( const VEC_TYPE& object_ ) : object(object_){}
3535
36- double get () const {
36+ double get () const {
3737 double
3838 m = mean (object).get (),
3939 ssq = sum ( pow (object, 2.0 ) ).get ();
40- return (ssq - m * m * object.size ()) / (object.size () - 1 );
41- }
40+ return (ssq - m * m * object.size ()) / (object.size () - 1 );
41+ }
4242
4343private:
44- const VEC_TYPE& object ;
44+ const VEC_TYPE& object ;
4545} ;
4646
4747template <bool NA, typename T>
4848class Var <CPLXSXP,NA,T> : public Lazy< double , Var<CPLXSXP,NA,T> > {
4949public:
5050 typedef typename Rcpp::VectorBase<CPLXSXP,NA,T> VEC_TYPE ;
5151
52- Var ( const VEC_TYPE& object_ ) : object(object_){}
52+ Var ( const VEC_TYPE& object_ ) : object(object_){}
5353
54- double get () const {
54+ double get () const {
5555 double sq = 0 , ssq = 0 ;
5656 for (int i = 0 ;i < object.size ();i++) {
5757 Rcomplex z = object[i];
5858 sq += z.r ;
5959 ssq += z.r * z.r ;
6060 }
61- return (ssq - sq * sq / object.size ()) / (object.size () - 1 );
62- }
61+ return (ssq - sq * sq / object.size ()) / (object.size () - 1 );
62+ }
6363
6464private:
65- const VEC_TYPE& object ;
65+ const VEC_TYPE& object ;
6666} ;
6767
6868} // sugar
6969
7070template <bool NA, typename T>
7171inline sugar::Var<REALSXP,NA,T> var ( const VectorBase<REALSXP,NA,T>& t){
72- return sugar::Var<REALSXP,NA,T>( t ) ;
72+ return sugar::Var<REALSXP,NA,T>( t ) ;
7373}
7474
7575template <bool NA, typename T>
0 commit comments