#1482 Problem assigning a JacobianTranspose to a FieldMatrix
Metadata
Property | Value |
---|---|
Reported by | Timo Betcke (timo.betcke@gmail.com) |
Reported at | Jul 25, 2014 11:23 |
Type | Bug Report |
Version | 2.3 |
Operating System | Unspecified / All |
Last edited by | Christoph Grüninger (gruenich@iws.uni-stuttgart.de) |
Last edited at | Jun 26, 2015 09:55 |
Closed by | Christoph Grüninger (gruenich@iws.uni-stuttgart.de) |
Closed at | Jun 26, 2015 09:55 |
Closed in version | Unknown |
Resolution | Won't fix |
Comment | Not a bug for Dune 2.4, we will not fix it for Dune 2.3. |
Description
Hi,
I have the following problem when assigning a JacobianTransposed to a FieldMatrix. The following code snippet works fine:
Dune::FieldMatrix<typename DuneGeometry::ctype,DuneGeometry::mydimension,DuneGeometry::coorddimension> duneJacobianT = duneGeo.jacobianTransposed(duneLocal);
But when I first define the variable 'duneJacobianT' and then assign the result of a call to 'jacobianTransposed' to it the resulting FieldMatrix is always filled with all ones entries and not with the correct values, that is the following code produces incorrect values in duneJacobianT
Dune::FieldMatrix<typename DuneGeometry::ctype,DuneGeometry::mydimension,DuneGeometry::coorddimension> duneJacobianT; duneJacobianT = duneGeo.jacobianTransposed(duneLocal);
It seems that the assignment operator between a Jacobian Transposed and a FieldMatrix has an issue in Dune 2.3.1 as copy constructing the object as done in the first code snippet works fine.
Best wishes
Timo