Cheetah - SKA - PSS - Prototype Time Domain Search Pipeline
Spectrum.cpp
1 /*
2  * The MIT License (MIT)
3  *
4  * Copyright (c) 2016 The SKA organisation
5  *
6  * Permission is hereby granted, free of charge, to any person obtaining a copy
7  * of this software and associated documentation files (the "Software"), to deal
8  * in the Software without restriction, including without limitation the rights
9  * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
10  * copies of the Software, and to permit persons to whom the Software is
11  * furnished to do so, subject to the following conditions:
12  *
13  * The above copyright notice and this permission notice shall be included in all
14  * copies or substantial portions of the Software.
15  *
16  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
17  * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
18  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
19  * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
20  * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
21  * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
22  * SOFTWARE.
23  */
24 #include "cheetah/rfim/ampp/Spectrum.h"
25 #include <algorithm>
26 #include <numeric>
27 
28 
29 namespace ska {
30 namespace cheetah {
31 namespace rfim {
32 namespace ampp {
33 
34 
35 template<typename NumericalRep, typename Alloc>
36 Spectrum<NumericalRep, Alloc>::Spectrum(pss::astrotypes::DimensionSize<data::Frequency> size)
37  : BaseT(size)
38 {
39 }
40 
41 template<typename NumericalRep, typename Alloc>
42 Spectrum<NumericalRep, Alloc>::~Spectrum()
43 {
44 }
45 
46 template<typename NumericalRep, typename Alloc>
48 {
49  //std::transform(_dataBandPass.begin(), _dataBandPass.end(), _dataBandPass.begin(), bind2nd(std::plus<float>(), offset));
50  float mean=0.0;
51  auto delta = v - _mean;
52  for(auto& val : *this) {
53  val += delta;
54  mean += val;
55  }
56  _mean = mean/this->template dimension<data::Frequency>();
57 }
58 
59 template<typename NumericalRep, typename Alloc>
61 {
62  return _mean;
63 }
64 
65 template<typename NumericalRep, typename Alloc>
67 {
68  _mean = std::accumulate(this->begin(), this->end(), 0.0)/this->template dimension<data::Frequency>();
69  return _mean;
70 }
71 
72 } // namespace ampp
73 } // namespace rfim
74 } // namespace cheetah
75 } // namespace ska
Some limits and constants for FLDO.
Definition: Brdz.h:35
A single dimension representation of the bandPass.
Definition: Spectrum.h:42