implement custom zero tensor

as long as it is not supported by Eigen

Conflicts:
	src/wrappers.hpp

src/eigen_tensor.hpp

@@ -0,0 +1,68 @@
+// ====================================================================
+// This file is part of FlexibleSUSY.
+//
+// FlexibleSUSY 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 3 of the License,
+// or (at your option) any later version.
+//
+// FlexibleSUSY 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 FlexibleSUSY.  If not, see
+// <http://www.gnu.org/licenses/>.
+// ====================================================================
+
+#ifndef EIGEN_TENSOR_H
+#define EIGEN_TENSOR_H
+
+#include <Eigen/Core>
+
+#if EIGEN_VERSION_AT_LEAST(3,3,0)
+
+#include <unsupported/Eigen/CXX11/Tensor>
+
+namespace flexiblesusy {
+
+template <class Scalar, int N, int M, int K>
+class ZeroTensor3 : public Eigen::TensorFixedSize<Scalar, Eigen::Sizes<N,M,K> > {
+   typedef Eigen::TensorFixedSize<Scalar, Eigen::Sizes<N,M,K> > Base;
+
+   using Base::setZero;
+
+   ZeroTensor3() : Base() {
+      setZero();
+   }
+};
+
+template <class Scalar, int N, int M, int K, int L>
+class ZeroTensor4 : public Eigen::TensorFixedSize<Scalar, Eigen::Sizes<N,M,K,L> > {
+   typedef Eigen::TensorFixedSize<Scalar, Eigen::Sizes<N,M,K,L> > Base;
+
+   using Base::setZero;
+
+   ZeroTensor4() : Base() {
+      setZero();
+   }
+};
+
+} // namespace flexiblesusy
+
+#else
+
+namespace flexiblesusy {
+
+template <class Scalar, int N, int M, int K>
+class ZeroTensor3;
+
+template <class Scalar, int N, int M, int K, int L>
+class ZeroTensor4;
+
+} // namespace flexiblesusy
+
+#endif
+
+#endif

src/module.mk

@@ -63,6 +63,7 @@ LIBFLEXI_HDR := \
 		$(DIR)/derivative.hpp \
 		$(DIR)/dilog.hpp \
 		$(DIR)/eigen_utils.hpp \
+		$(DIR)/eigen_tensor.hpp \
 		$(DIR)/error.hpp \
 		$(DIR)/ew_input.hpp \
 		$(DIR)/ewsb_solver.hpp \

src/wrappers.hpp

@@ -30,6 +30,7 @@
 #include <boost/lexical_cast.hpp>
 
 #include "dilog.hpp"
+#include "eigen_tensor.hpp"
 
 namespace flexiblesusy {
 
@@ -466,11 +467,15 @@ void Symmetrize(Eigen::MatrixBase<Derived>& m)
 
 #define UNITMATRIX(rows)             Eigen::Matrix<double,rows,rows>::Identity()
 #define ZEROMATRIX(rows,cols)        Eigen::Matrix<double,rows,cols>::Zero()
+#define ZEROTENSOR3(d1,d2,d3)        ZeroTensor3<double,d1,d2,d3>
+#define ZEROTENSOR4(d1,d2,d3,d4)     ZeroTensor4<double,d1,d2,d3,d4>
 #define ZEROVECTOR(rows)             Eigen::Matrix<double,rows,1>::Zero()
 #define ZEROARRAY(rows)              Eigen::Array<double,rows,1>::Zero()
 #define UNITMATRIXCOMPLEX(rows)      Eigen::Matrix<std::complex<double>,rows,rows>::Identity()
 #define ZEROMATRIXCOMPLEX(rows,cols) Eigen::Matrix<std::complex<double>,rows,cols>::Zero()
 #define ZEROVECTORCOMPLEX(rows)      Eigen::Matrix<std::complex<double>,rows,1>::Zero()
+#define ZEROTENSOR3COMPLEX(d1,d2,d3) ZeroTensor3<std::complex<double>,d1,d2,d3>
+#define ZEROTENSOR4COMPLEX(d1,d2,d3,d4) ZeroTensor4<std::complex<double>,d1,d2,d3,d4>
 #define ZEROARRAYCOMPLEX(rows)       Eigen::Array<std::complex<double>,rows,1>::Zero()
 
 // MxN matrix projection operator, which projects on the (X,Y)