EigenExpValues.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Dense>
template<class Derived1, class Derived2>
inline auto EigenExpValues(Eigen::MatrixBase<Derived1> const& eigenvectors,
                           Eigen::MatrixBase<Derived2> const& observable)
    -> Eigen::VectorX<typename Derived1::RealScalar> {
    return (eigenvectors.adjoint() * observable * eigenvectors).diagonal().array().real();
}
 
#ifdef GPU
    #include "EigenExpValues.cuh"
#endif