[Dune] [#1054] use maxima instead of maple to generate quadrature rules (Attachment added)

Dune flyspray at dune-project.org
Mon Mar 19 19:30:23 CET 2012


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

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

User who did this - Ansgar Burchardt (ansgar) 

Attached to Project - Dune
Summary - use maxima instead of maple to generate quadrature rules
Task Type - Bug Report
Category - Quadrature
Status - New
Assigned To - 
Operating System - Unspecified / All
Severity - Low
Priority - Normal
Reported Version - SVN (pre2.2)
Due in Version - Undecided
Due Date - Undecided
Details - I rewrote jacobian.mw in maxima (now jacobian.mac) so that dune-geometry can be built using only free software tools.  There are small differences in the coefficients, but only at the last five decimals (of 100); also the coefficients now are written with leading zeroes instead of scientific notation (ie. 0.000123 instead of 1.23e-4).

I did run the tests with `make check' and there was an error, but the same error was also there with the old implementation, but please double-check this before committing the changes:

Error: Quadrature for (cube, 4) and order=16 does not sum to volume of RefElem
	Sums to 1( RefElem.volume() = 1)(difference -5.71765e-14)
Error: Quadrature for (cube, 4) and order=24 does not sum to volume of RefElem
	Sums to 1( RefElem.volume() = 1)(difference 9.68114e-14)
Error: Quadrature for (cube, 4) and order=25 does not sum to volume of RefElem
	Sums to 1( RefElem.volume() = 1)(difference 9.68114e-14)

Regards,
Ansgar


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=1054

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.



More information about the Dune mailing list