1 | //#---------------------------------------------------------------------------
|
---|
2 | //# SDMathWrapper.h: Wrapper classes to use CountedPtr
|
---|
3 | //#---------------------------------------------------------------------------
|
---|
4 | //# Copyright (C) 2004
|
---|
5 | //# ATNF
|
---|
6 | //#
|
---|
7 | //# This program is free software; you can redistribute it and/or modify it
|
---|
8 | //# under the terms of the GNU General Public License as published by the Free
|
---|
9 | //# Software Foundation; either version 2 of the License, or (at your option)
|
---|
10 | //# any later version.
|
---|
11 | //#
|
---|
12 | //# This program is distributed in the hope that it will be useful, but
|
---|
13 | //# WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
14 | //# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
|
---|
15 | //# Public License for more details.
|
---|
16 | //#
|
---|
17 | //# You should have received a copy of the GNU General Public License along
|
---|
18 | //# with this program; if not, write to the Free Software Foundation, Inc.,
|
---|
19 | //# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
|
---|
20 | //#
|
---|
21 | //# Correspondence concerning this software should be addressed as follows:
|
---|
22 | //# Internet email: Malte.Marquarding@csiro.au
|
---|
23 | //# Postal address: Malte Marquarding,
|
---|
24 | //# Australia Telescope National Facility,
|
---|
25 | //# P.O. Box 76,
|
---|
26 | //# Epping, NSW, 2121,
|
---|
27 | //# AUSTRALIA
|
---|
28 | //#
|
---|
29 | //# $Id:
|
---|
30 | //#---------------------------------------------------------------------------
|
---|
31 | #ifndef SDMATHWRAPPER_H
|
---|
32 | #define SDMATHWRAPPER_H
|
---|
33 |
|
---|
34 | #include <vector>
|
---|
35 | #include <string>
|
---|
36 |
|
---|
37 | #include <boost/python/tuple.hpp>
|
---|
38 |
|
---|
39 | #include "SDMemTableWrapper.h"
|
---|
40 | #include "SDMath.h"
|
---|
41 |
|
---|
42 | namespace asap {
|
---|
43 |
|
---|
44 | namespace SDMathWrapper {
|
---|
45 |
|
---|
46 | // Quotient
|
---|
47 |
|
---|
48 | SDMemTableWrapper quotient(const SDMemTableWrapper& on,
|
---|
49 | const SDMemTableWrapper& off) {
|
---|
50 | return SDMemTableWrapper(SDMath::quotient(on.getCP(),
|
---|
51 | off.getCP()));
|
---|
52 | }
|
---|
53 |
|
---|
54 | // Multiply
|
---|
55 |
|
---|
56 | void scaleInSitu(SDMemTableWrapper& in, casa::Float factor, casa::Bool all)
|
---|
57 | {
|
---|
58 | SDMemTable* sdmt = in.getPtr();
|
---|
59 | SDMath::multiplyInSitu(in.getPtr(),factor, all);
|
---|
60 | }
|
---|
61 | SDMemTableWrapper scale(const SDMemTableWrapper& in,
|
---|
62 | casa::Float factor, casa::Bool all)
|
---|
63 | {
|
---|
64 | return SDMemTableWrapper(SDMath::multiply(in.getCP(), factor, all));
|
---|
65 | }
|
---|
66 |
|
---|
67 | // Add
|
---|
68 |
|
---|
69 | void addInSitu(SDMemTableWrapper& in, casa::Float offset, casa::Bool all)
|
---|
70 | {
|
---|
71 | SDMemTable* sdmt = in.getPtr();
|
---|
72 | SDMath::addInSitu(in.getPtr(), offset, all);
|
---|
73 | }
|
---|
74 | SDMemTableWrapper add(const SDMemTableWrapper& in, casa::Float offset, casa::Bool all)
|
---|
75 | {
|
---|
76 | return SDMemTableWrapper(SDMath::add(in.getCP(), offset, all));
|
---|
77 | }
|
---|
78 |
|
---|
79 | // Hanning
|
---|
80 |
|
---|
81 | SDMemTableWrapper hanning(const SDMemTableWrapper& in) {
|
---|
82 | return SDMemTableWrapper(SDMath::hanning(in.getCP()));
|
---|
83 | }
|
---|
84 |
|
---|
85 | // Bin up
|
---|
86 |
|
---|
87 | SDMemTableWrapper bin(const SDMemTableWrapper& in,
|
---|
88 | int width) {
|
---|
89 | return SDMath::bin(in.getCP(), width);
|
---|
90 | }
|
---|
91 |
|
---|
92 | // Average in time
|
---|
93 |
|
---|
94 | SDMemTableWrapper average (boost::python::tuple tpl,
|
---|
95 | const std::vector<bool>& mask,
|
---|
96 | bool scanAv, const std::string& wt);
|
---|
97 |
|
---|
98 | // Average polarizations
|
---|
99 |
|
---|
100 | void averagePolInSitu (SDMemTableWrapper& in, const std::vector<bool>& mask)
|
---|
101 | {
|
---|
102 | SDMemTable* sdmt = in.getPtr();
|
---|
103 | SDMath::averagePolInSitu(in.getPtr(), mask);
|
---|
104 | }
|
---|
105 | SDMemTableWrapper averagePol (const SDMemTableWrapper& in, const std::vector<bool>& mask)
|
---|
106 | {
|
---|
107 | return SDMemTableWrapper(SDMath::averagePol(in.getCP(), mask));
|
---|
108 | }
|
---|
109 |
|
---|
110 | // Statistics
|
---|
111 |
|
---|
112 | std::vector<float> statistic(const SDMemTableWrapper& in,
|
---|
113 | const std::vector<bool>& mask, const std::string& which) {
|
---|
114 | return SDMath::statistic(in.getCP(), mask, which);
|
---|
115 | }
|
---|
116 |
|
---|
117 | };
|
---|
118 |
|
---|
119 | } // namespace
|
---|
120 | #endif
|
---|