From 3388cdd1b654242b391fe25fb7dcc56f3cb20d8e Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 31 Aug 2010 12:59:00 +0000
Subject: [PATCH] bugfix: can't assign a row vector to a column vector
[[Imported from SVN: r6294]]
---
dune/gfe/unitvector.hh | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh
index 03e9b36a..3283955f 100644
--- a/dune/gfe/unitvector.hh
+++ b/dune/gfe/unitvector.hh
@@ -164,7 +164,10 @@ public:
row[0] = b.globalCoordinates();
row *= secondDerivativeOfArcCosSquared(sp);
- Dune::FieldMatrix<double,dim,1> column = b.projectOntoTangentSpace(a.globalCoordinates());
+ Dune::FieldVector<double,dim> tmp = b.projectOntoTangentSpace(a.globalCoordinates());
+ Dune::FieldMatrix<double,dim,1> column;
+ for (int i=0; i<dim; i++) // turn row vector into column vector
+ column[i] = tmp[i];
Dune::FieldMatrix<double,dim,dim> A;
// A = row * column
--
GitLab