// This file is part of Eigen, a lightweight C++ template library // for linear algebra. Eigen itself is part of the KDE project. // // Copyright (C) 2006-2007 Benoit Jacob // // Eigen 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 2 or (at your option) any later version. // // Eigen 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 Eigen; if not, write to the Free Software Foundation, Inc., 51 // Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. // // As a special exception, if other files instantiate templates or use macros // or functions from this file, or you compile this file and link it // with other works to produce a work based on this file, this file does not // by itself cause the resulting work to be covered by the GNU General Public // License. This exception does not invalidate any other reasons why a work // based on this file might be covered by the GNU General Public License. #ifndef EIGEN_DOT_H #define EIGEN_DOT_H template struct DotUnroller { static void run(const Derived1 &v1, const Derived2& v2, typename Derived1::Scalar &dot) { DotUnroller::run(v1, v2, dot); dot += v1.coeff(Index) * conj(v2.coeff(Index)); } }; template struct DotUnroller<0, Size, Derived1, Derived2> { static void run(const Derived1 &v1, const Derived2& v2, typename Derived1::Scalar &dot) { dot = v1.coeff(0) * conj(v2.coeff(0)); } }; template struct DotUnroller { static void run(const Derived1&, const Derived2&, typename Derived1::Scalar&) {} }; // prevent buggy user code from causing an infinite recursion template struct DotUnroller { static void run(const Derived1&, const Derived2&, typename Derived1::Scalar&) {} }; template template Scalar MatrixBase::dot(const OtherDerived& other) const { assert(IsVector && OtherDerived::IsVector && size() == other.size()); Scalar res; if(EIGEN_UNROLLED_LOOPS && SizeAtCompileTime != Dynamic && SizeAtCompileTime <= 16) DotUnroller ::run(*static_cast(this), other, res); else { res = (*this).coeff(0) * conj(other.coeff(0)); for(int i = 1; i < size(); i++) res += (*this).coeff(i)* conj(other.coeff(i)); } return res; } template typename NumTraits::Real MatrixBase::norm2() const { return real(dot(*this)); } template typename NumTraits::Real MatrixBase::norm() const { return sqrt(norm2()); } template ScalarMultiple MatrixBase::normalized() const { return (*this) / norm(); } #endif // EIGEN_DOT_H