[Dune] Geometry.volume() paradigm for CurvilinearGeometry

[Aleksejs Fomins]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hi Andreas,

This "setGeometryTolerance" seems simple enough. I like it.

Best,
Aleksejs



On 23/02/15 12:13, Andreas Dedner wrote:
> Hi. I would provide the grid with a tolerance method (either a 
> "setGeometryTolerance( tolerance ) method or in the constructor). 
> Changing the volume method I don't consider an option. Best 
> Andreas
> 
> 
> On 23/02/15 11:09, Aleksejs Fomins wrote: Dear Dune,
> 
> I don't remember if I asked this question before, but it bothers
> me.
> 
> In CurvilinearGeometry, all entities are non-linear, and
> therefore, finding their volume requires an iterative procedure.
> Naturally, such computation will have finite precision. That is,
> the method geometry.volume() has to be provided with a constant
> specifying the relative integration error. Currently this is done
> internally.
> 
> The problem is that the integration time depends strongly on the 
> expected error, and, therefore, the user should have access to
> this parameter, to specify precision satisfactory for their needs.
> 
> Proposed solutions: 1) geometry.volume() becomes
> geometry.volume(tolerance)
> 
> 2) geometry.volume() returns linear volume as in other grids.
> Geometry class gains an extra method
> geometry.curvVolume(tolerance).
> 
> 3) extra command line constant (like config.h)
> 
> What are your suggestions?
> 
> Thanks, Aleksejs
>> 
>> _______________________________________________ Dune mailing
>> list Dune at dune-project.org 
>> http://lists.dune-project.org/mailman/listinfo/dune
> 
> 
> _______________________________________________ Dune mailing list 
> Dune at dune-project.org 
> http://lists.dune-project.org/mailman/listinfo/dune
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.12 (GNU/Linux)

iQIcBAEBAgAGBQJU6xDEAAoJEDkNM7UOEwMZSMkQANbdKixxzapFQetoCEw7FpiZ
FcQMa/XWXGfcNliJKNg6ZyT1gQ1UwFTpbJCryyUGtHjcHVTcLkBBIIPp5xtDMuTi
G0WMc3LbpJSbaUmmHXcVM2kUSEX2DK5x7GI/yOWcEYmQsq+Yusf/9qbU1TnQZ5fy
0Ue0yIaTSnWrCe8RdPkeirw4TnO/u+2JL/1BhBiJQG7Zqwu2tOmB4EHMUEWI+cql
xICv8o5dk70gQAdPKMkpyYlZhoHF4ymvoAZvU3gxvd9+SiySDre8jh3hwfbBN1T9
zE9aS/n+W/XV14pzY8oOqM6VrbBMgbxCnBZqUq8YauGojlmw91cMcc6kmW75nT9z
CVTz9hRfdMLuFuu6QCYUPeU9lLpliG1MZEvboL8du2wa3wtYEBsJtMcyibvJ82g6
RVIsjrgiA6+iNIlaxAO4MtZUXA9LJWgf9OGG58Pp17fNm9XtNOtSJl7js3hY+rCG
s7ZLwnvA11zS0ksCCtdSqWeD9/5cwArxKteWo6DnteCamunUH153umjyFv3DjEG6
d7CSwNE/AqRiL0k81qz8gkRKVmQ9kLrqce+y/mK8il2741xF6juXsOYs51A861qQ
c+YoQTApSW5M+4DS75VmIDZ1obVZf3gaVXB0r7cVj/eyQVtoghIc/kFTmqY4KY8s
x6obLat/F+LkuXjd2BbV
=I2F6
-----END PGP SIGNATURE-----