THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.
A new Flyspray task has been opened. Details are below.
User who did this - Bernd Flemisch (bernd)
Attached to Project - DuMuX
Summary - Geometries and Jacobians
Task Type - Bug Report
Category - General
Status - New
Assigned To -
Operating System - Unspecified
Severity - Low
Priority - Normal
Reported Version - SVN
Due in Version - Undecided
Due Date - Undecided
Details - The issue emerged when compiling test/decoupled/1p/test_1p with
optimization options >=O1. In this case the call
const JacobianTransposed& jacobianT =
eIt->geometry().jacobianTransposed(localPos);
in fvvelocity1p.hh resulted in an all-zero matrix. We can find a partial
explanation by Jö in Dune's FS#1174:
"If geometry() returns a copy and jacobianInverseTransposed() returns a
reference, then saying
const Jacobian &jac = (*it).geometry().jacobianInverseTransposed(quad(i));
will actually most likely be wrong. Your either have to do
const Geometry &geo = (*it).geometry();
const Jacobian &jac = geo.jacobianInverseTransposed(quad(i));
or
Jacobian jac = (*it).geometry().jacobianInverseTransposed(quad(i));"
The discussion resulted in returning the Jacobians by copy. Unfortunately, this
did not find its way into Dune 2.3.
We should look for other occurrences and fix them as soon as possible.
More information can be found at the following URL:
http://www.dumux.org/flyspray/index.php?do=details&task_id=240
You are receiving this message because you have requested it from the Flyspray
bugtracking system. If you did not expect this message or don't want to
receive mails in future, you can change your notification settings at the URL
shown above.
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
