var_strategy.hpp

Go to the documentation of this file.

/*
 * Copyright (C) 1998-2018 ALPS Collaboration. See COPYRIGHT.TXT
 * All rights reserved. Use is subject to license terms. See LICENSE.TXT
 * For use in publications, see ACKNOWLEDGE.TXT
 */
#pragma once

#include <Eigen/Core>

#include <alps/alea/complex_op.hpp>

// Forward declarations
namespace alps { namespace alea {
    struct circular_var;
    struct elliptic_var;
    template <typename Strategy, typename T> struct bind;
}}

// Actual declarations

namespace alps { namespace alea {

struct circular_var { };

struct elliptic_var { };


template <typename T>
struct bind<circular_var, T>
{
    typedef T value_type;
    typedef T var_type;
    typedef T cov_type;

    typedef Eigen::internal::scalar_abs2_op<T> abs2_op;

    static cov_type outer(T x, T y) { return x * y; }
};

template <typename T>
struct bind<circular_var, std::complex<T> >
{
    typedef std::complex<T> value_type;
    typedef T var_type;
    typedef std::complex<T> cov_type;

    typedef Eigen::internal::scalar_abs2_op<std::complex<T> > abs2_op;

    static cov_type outer(std::complex<T> x, std::complex<T> y)
    { return x * std::conj(y); }
};

template <typename T>
struct bind<elliptic_var, T>
{
    typedef T value_type;
    typedef T var_type;
    typedef T cov_type;

    typedef Eigen::internal::scalar_abs2_op<T> abs2_op;

    static cov_type outer(T x, T y) { return x * y; }
};

template <typename T>
struct bind<elliptic_var, std::complex<T> >
{
    typedef std::complex<T> value_type;
    typedef complex_op<T> var_type;
    typedef complex_op<T> cov_type;

    struct abs2_op
    {
        typedef complex_op<T> result_type;
        const complex_op<T> operator() (const std::complex<T> &x) const
        { return complex_op<T>::outer(x, x); }
    };

    static complex_op<T> outer(std::complex<T> x, std::complex<T> y)
    { return complex_op<T>::outer(x, y); }
};


}} /* namespace alps::alea */