[Dune] [#1099] evaluateJacobian seems to be wrong for Q1LocalBasis (Attachment added)

[Dune]

THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

A new Flyspray task has been opened.  Details are below. 

User who did this - Andreas Lauser (andlaus) 

Attached to Project - Dune
Summary - evaluateJacobian seems to be wrong for Q1LocalBasis
Task Type - Bug Report
Category - Localfunctions
Status - New
Assigned To - 
Operating System - Unspecified / All
Severity - Low
Priority - Normal
Reported Version - 2.1
Due in Version - Undecided
Due Date - Undecided
Details - I've just noticed that evaluateJacobian() for Q1 elements seem to produce jacobian vectors which depend on the position. if I'm not completely wrong (which is quite possible), this finite element is supposed to be linear in respect to all coordinates and the local gradients should be feature either -1 or 1 for all coordinates. Taking the partial derivatives on in[j] in q1localbasis.hh is also consistent with this view. I've attached a patch for dune-localfunctions (against the current dune 2.1 SVN release branch) to fix this
issue. I suppose that since this issue is not present at the corners of the element, the gradients of this finite element have probably never been tested in the interior. On the other hand I might just miss something; in this case I would like to apologize for the noise.

One or more files have been attached.

More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=1099

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.