From d1def995a46f69c346b9ffbbfb466ac944010d99 Mon Sep 17 00:00:00 2001
From: young <young@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 13 Feb 2014 09:51:05 +0000
Subject: [PATCH] A patch by Denis Davydov: Multiply a real tensor by a complex
 number.

git-svn-id: https://svn.dealii.org/trunk@32475 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/include/deal.II/base/tensor.h      | 42 ++++++++++++++++++++++
 deal.II/include/deal.II/base/tensor_base.h | 41 +++++++++++++++++++++
 2 files changed, 83 insertions(+)

diff --git a/deal.II/include/deal.II/base/tensor.h b/deal.II/include/deal.II/base/tensor.h
index 5ddc7dc752..c3c6dab702 100644
--- a/deal.II/include/deal.II/base/tensor.h
+++ b/deal.II/include/deal.II/base/tensor.h
@@ -2093,6 +2093,48 @@ operator / (const Tensor<rank,dim> &t,
   return tt;
 }
 
+
+
+/**
+ * Multiplication of a tensor of general rank by a scalar
+ * complex<double> from the left.
+ *
+ * @relates Tensor
+ */
+template <int rank, int dim>
+  inline
+  Tensor<rank,dim,std::complex<double> >
+  operator * (const std::complex<double>  factor,
+	      const Tensor<rank,dim>     &t)
+{
+  Tensor<rank,dim,std::complex<double> > tt;
+  for (unsigned int d=0; d<dim; ++d)
+    tt[d] = factor * t[d];
+  return tt;
+}
+
+
+
+/**
+ * Multiplication of a tensor of general rank by a scalar
+ * complex<double> from the right.
+ *
+ * @relates Tensor
+ */
+template <int rank, int dim>
+  inline
+  Tensor<rank,dim,std::complex<double> >
+  operator * (const Tensor<rank,dim>     &t,
+	      const std::complex<double>  factor)
+{
+  Tensor<rank,dim,std::complex<double> > tt;
+  for (unsigned int d=0; d<dim; ++d)
+    tt[d] = t[d] * factor;
+  return tt;
+}
+
+
+
 DEAL_II_NAMESPACE_CLOSE
 
 #endif
diff --git a/deal.II/include/deal.II/base/tensor_base.h b/deal.II/include/deal.II/base/tensor_base.h
index 94c8e1b3e9..1b3c58fe60 100644
--- a/deal.II/include/deal.II/base/tensor_base.h
+++ b/deal.II/include/deal.II/base/tensor_base.h
@@ -1463,6 +1463,47 @@ operator / (const Tensor<1,dim> &t,
 }
 
 
+
+/**
+ * Multiplication of a tensor of rank 1 by a scalar complex<double>
+ * from the right.
+ *
+ * @relates Tensor<1,dim,Number>
+ */
+template <int dim>
+inline
+Tensor<1,dim,std::complex<double> >
+  operator * (const Tensor<1,dim>        &t,
+	      const std::complex<double>  factor)
+{                                  
+  Tensor<1,dim,std::complex<double> > tt (false);
+  for (unsigned int d=0; d<dim; ++d)
+    tt[d] = t[d] * factor;
+  return tt;
+}
+
+
+
+/**
+ * Multiplication of a tensor of rank 1 by a scalar complex<double>
+ * from the left.
+ *
+ * @relates Tensor<1,dim,Number>
+ */
+template <int dim>
+inline
+Tensor<1,dim,std::complex<double> >
+  operator * (const std::complex<double>  factor,
+	      const Tensor<1,dim>        &t)
+{                                  
+  Tensor<1,dim,std::complex<double> > tt (false);
+  for (unsigned int d=0; d<dim; ++d)
+    tt[d] = t[d] * factor;
+  return tt;
+}
+
+
+
 DEAL_II_NAMESPACE_CLOSE
 
 #endif
-- 
2.39.5