/usr/include/mia-2.4/mia/template/multicost.cxx is in libmia-2.4-dev 2.4.6-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 | /* -*- mia-c++ -*-
*
* This file is part of MIA - a toolbox for medical image analysis
* Copyright (c) Leipzig, Madrid 1999-2017 Gert Wollny
*
* MIA 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 3 of the License, or
* (at your option) any later version.
*
* This program 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 MIA; if not, see <http://www.gnu.org/licenses/>.
*
*/
#include <sstream>
NS_MIA_BEGIN
template <typename T>
TFullCostList<T>::TFullCostList():
TFullCost<T>(1.0)
{
}
template <typename T>
void TFullCostList<T>::push(typename TFullCost<T>::Pointer cost)
{
m_costs.push_back(cost);
}
template <typename T>
bool TFullCostList<T>::do_has(const char *property) const
{
TRACE_FUNCTION;
bool result = !m_costs.empty();
if (!result)
cvwarn() << "No cost functions given\n";
auto ic = m_costs.begin();
while (result && ic != m_costs.end()) {
result &= (*ic)->has(property);
if (!result) {
cvwarn() << "Cost '"<< (*ic)->get_init_string()
<< "' doesn't provide property '"
<< property << "'\n";
}
++ic;
}
return result;
}
template <typename T>
double TFullCostList<T>::do_evaluate(const T& t, CDoubleVector& gradient) const
{
double result = 0;
CDoubleVector tmp(gradient.size());
std::stringstream msg;
msg << "Cost: ";
for (auto i = m_costs.begin(); i != m_costs.end(); ++i) {
std::fill(tmp.begin(), tmp.end(), 0.0);
double h = (*i)->evaluate(t, tmp);
msg << h << "("<< (*i)->get_init_string() << ") ";
result += h;
std::transform(gradient.begin(), gradient.end(), tmp.begin(), gradient.begin(),
[](double x, double y){return x+y;});
}
cvinfo() << msg.str() << " = " << result << "\n";
return result;
}
template <typename T>
double TFullCostList<T>::do_value(const T& t) const
{
double result = 0;
std::stringstream msg;
msg << "Cost: ";
for (auto i = m_costs.begin(); i != m_costs.end(); ++i) {
double h = (*i)->cost_value(t);
msg << h << "("<< (*i)->get_init_string() << ") ";
result += h;
}
cvinfo() << msg.str() << " = " << result << "\n";
return result;
}
template <typename T>
double TFullCostList<T>::do_value() const
{
double result = 0;
for (auto i = m_costs.begin(); i != m_costs.end(); ++i) {
result += (*i)->cost_value();
}
return result;
}
template <typename T>
void TFullCostList<T>::do_set_size()
{
for (auto i = m_costs.begin(); i != m_costs.end(); ++i)
(*i)->set_size(this->get_current_size());
}
template <typename T>
bool TFullCostList<T>::do_get_full_size(Size& size) const
{
bool result = true;
for (auto i = m_costs.begin(); i != m_costs.end() && result; ++i)
result = (*i)->get_full_size(size);
return result;
}
template <typename T>
void TFullCostList<T>::do_reinit()
{
for (auto i = m_costs.begin(); i != m_costs.end(); ++i)
(*i)->reinit();
}
NS_MIA_END
|