#1054 use maxima instead of maple to generate quadrature rules
Metadata
| Property | Value |
|---|---|
| Reported by | Ansgar Burchardt (burchardt@igpm.rwth-aachen.de) |
| Reported at | Mar 19, 2012 18:30 |
| Type | Bug Report |
| Version | Git (pre2.4) [autotools] |
| Operating System | Unspecified / All |
| Last edited by | Christian Engwer (christi@conan.iwr.uni-heidelberg.de) |
| Last edited at | Mar 29, 2012 10:10 |
| Closed by | Christian Engwer (christi@conan.iwr.uni-heidelberg.de) |
| Closed at | Mar 29, 2012 10:10 |
| Closed in version | 2.2 |
| Resolution | Implemented |
| Comment |
Description
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