VectorT_dummy_data.hpp

#pragma once
 
#include <valarray>
#include <complex>
 
template<typename T>
T make_dummy_component() {
    return T();
}
 
template<> std::valarray<double> make_dummy_component();
 
template<typename Vec>
Vec make_dummy_vector() {
    return Vec(make_dummy_component<typename Vec::value_type>());
}
 
template<typename Scalar>
double observe_dummy_component(const Scalar& _s) {
    return _s;
}
 
template<typename T>
double observe_dummy_component(const std::valarray<T>& _s) {
    return _s.sum();
}
 
template<typename T>
double observe_dummy_component(const std::complex<T>& _s) {
    return _s.real() + _s.imag();
}
 
template<typename Vec>
double observe_dummy_vector(const Vec& _vec) {
    double result = 0.0;
    for (int dim = 0; dim < _vec.dim(); ++dim) {
        result += observe_dummy_component(_vec[dim]);
    }
    return result;
}