From: Maien Hamed <tomaien@hotmail.com>
Date: Wed, 19 Aug 2015 14:29:23 +0000 (+0200)
Subject: Generalized conversions between DerivativeForm<rank,dim,spacedim> and
X-Git-Tag: v8.4.0-rc2~586^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F1389%2Fhead;p=dealii.git

Generalized conversions between DerivativeForm<rank,dim,spacedim> and
Tensor<rank,spacedim> to general ranks.
---

diff --git a/doc/news/changes.h b/doc/news/changes.h
index 73000a826a..2f59813705 100644
--- a/doc/news/changes.h
+++ b/doc/news/changes.h
@@ -134,6 +134,13 @@ inconvenience this causes.
 
 
 <ol>
+  <li> Improved: Generalized conversion between Tensor<order+1,dim> and
+  DerivativeForm<order,dim,dim> to general order using converting constructor
+  and assignment operator.
+  <br>
+  (Maien Hamed, 2015/08/01-2015/08/09)
+  </li>
+
   <li> Changed: The function Vector::add() that adds a scalar number to all
   elements of a vector has been deprecated. The same is true for the
   Vector::ratio() function, and for the corresponding functions in other
diff --git a/include/deal.II/base/derivative_form.h b/include/deal.II/base/derivative_form.h
index 961e900ef1..a93c843630 100644
--- a/include/deal.II/base/derivative_form.h
+++ b/include/deal.II/base/derivative_form.h
@@ -63,9 +63,9 @@ public:
   DerivativeForm ();
 
   /**
-   * Constructor from a second order tensor.
+   * Constructor from a tensor.
    */
-  DerivativeForm (const Tensor<2,dim,Number> &);
+  DerivativeForm (const Tensor<order+1,dim,Number> &);
 
   /**
    * Read-Write access operator.
@@ -86,7 +86,7 @@ public:
   /**
    * Assignment operator.
    */
-  DerivativeForm   &operator = (const Tensor<2,dim, Number> &);
+  DerivativeForm   &operator = (const Tensor<order+1,dim, Number> &);
 
   /**
    * Assignment operator.
@@ -94,10 +94,11 @@ public:
   DerivativeForm   &operator = (const Tensor<1,dim, Number> &);
 
   /**
-   * Converts a DerivativeForm <1,dim, dim> to Tensor<2,dim,Number>. If the
-   * derivative is the Jacobian of F, then Tensor[i] = grad(F^i).
+   * Converts a DerivativeForm <order,dim,dim> to Tensor<order+1,dim,Number>.
+   * In particular, if order==1 and the derivative is the Jacobian of F, then
+   * Tensor[i] = grad(F^i).
    */
-  operator Tensor<2,dim,Number>() const;
+  operator Tensor<order+1,dim,Number>() const;
 
   /**
    * Converts a DerivativeForm <1, dim, 1> to Tensor<1,dim,Number>.
@@ -180,11 +181,11 @@ DerivativeForm<order, dim, spacedim, Number>::DerivativeForm  ()
 
 template <int order, int dim, int spacedim, typename Number>
 inline
-DerivativeForm<order, dim, spacedim, Number>::DerivativeForm(const Tensor<2,dim,Number> &T)
+DerivativeForm<order, dim, spacedim, Number>::DerivativeForm(const Tensor<order+1,dim,Number> &T)
 {
-  Assert( (dim == spacedim) && (order==1),
-          ExcMessage("Only allowed for square tensors."));
-  if ((dim == spacedim) && (order==1))
+  Assert( (dim == spacedim),
+          ExcMessage("Only allowed for forms with dim==spacedim."));
+  if (dim == spacedim)
     for (unsigned int j=0; j<dim; ++j)
       (*this)[j] = T[j];
 }
@@ -207,12 +208,12 @@ operator = (const DerivativeForm<order, dim, spacedim, Number> &ta)
 template <int order, int dim, int spacedim, typename Number>
 inline
 DerivativeForm<order, dim, spacedim, Number> &DerivativeForm<order, dim, spacedim, Number>::
-operator = (const Tensor<2,dim,Number> &ta)
+operator = (const Tensor<order+1,dim,Number> &ta)
 {
-  Assert( (dim == spacedim) && (order==1),
-          ExcMessage("Only allowed for square tensors."));
+  Assert( (dim == spacedim),
+          ExcMessage("Only allowed when dim==spacedim."));
 
-  if ((dim == spacedim) && (order==1))
+  if (dim == spacedim)
     for (unsigned int j=0; j<dim; ++j)
       (*this)[j] = ta[j];
   return *this;
@@ -276,14 +277,14 @@ DerivativeForm<order, dim, spacedim, Number>::operator Tensor<1,dim,Number>() co
 
 template <int order, int dim, int spacedim, typename Number>
 inline
-DerivativeForm<order, dim, spacedim, Number>::operator Tensor<2,dim,Number>() const
+DerivativeForm<order, dim, spacedim, Number>::operator Tensor<order+1,dim,Number>() const
 {
-  Assert( (dim == spacedim) && (order==1),
-          ExcMessage("Only allowed for square tensors."));
+  Assert( (dim == spacedim),
+          ExcMessage("Only allowed when dim==spacedim."));
 
-  Tensor<2,dim,Number> t;
+  Tensor<order+1,dim,Number> t;
 
-  if ((dim == spacedim) && (order==1))
+  if (dim == spacedim)
     for (unsigned int j=0; j<dim; ++j)
       t[j] = (*this)[j];