l_ew.h

Go to the documentation of this file.

/* $Id: l_ew.h 17249 2014-08-07 16:03:32Z predoehl $ */
/* =========================================================================== *
   |
   |  Copyright (c) 1994-2010 by Kobus Barnard (author)
   |
   |  Personal and educational use of this code is granted, provided that this
   |  header is kept intact, and that the authorship is not misrepresented, that
   |  its use is acknowledged in publications, and relevant papers are cited.
   |
   |  For other use contact the author (kobus AT cs DOT arizona DOT edu).
   |
   |  Please note that the code in this file has not necessarily been adequately
   |  tested. Naturally, there is no guarantee of performance, support, or fitness
   |  for any particular task. Nonetheless, I am interested in hearing about
   |  problems that you encounter.
   |
   |  Author:  Ernesto Brau
 * =========================================================================== */

#ifndef KJB_L_EW_DEFINED
#define KJB_L_EW_DEFINED

namespace kjb
{

template<class Indexable>
Indexable ew_multiply(const Indexable& I, const Indexable& J)
{
    Indexable K = I;
    for(int i = 0; i < I.get_length(); i++)
    {
        K(i) = K(i) * J(i);
    }

    return K;
}

template<class Indexable>
void ew_abs_ow(Indexable& I)
{
    for(int i = 0; i < I.get_length(); i++)
    {
        I(i) = abs(I(i));
    }
}

template<class Indexable>
void ew_square_ow(Indexable& I)
{
    for(int i = 0; i < I.get_length(); i++)
    {
        I(i) = I(i) * I(i);
    }
}

template<class Indexable>
Indexable ew_square(const Indexable& I)
{
    Indexable M = I;
    ew_square_ow(M);
    return M;
}

template<class Indexable>
void ew_sqrt_ow(Indexable& I)
{
    for(int i = 0; i < I.get_length(); i++)
    {
        using ::sqrt;
        I(i) = sqrt(I(i));
    }
}

template<class Indexable>
Indexable ew_sqrt(const Indexable& I)
{
    Indexable M = I;
    ew_sqrt_ow(M);
    return M;
}


} //namespace kjb

#endif /* KJB_L_EW_DEFINED */