]> https://gitweb.dealii.org/ - dealii.git/commit
Implement a 2D-specific chart point estimate.
authorDavid Wells <wellsd2@rpi.edu>
Thu, 18 Jan 2018 18:37:43 +0000 (13:37 -0500)
committerDavid Wells <wellsd2@rpi.edu>
Fri, 19 Jan 2018 16:44:31 +0000 (11:44 -0500)
commit7a2c1d182e9d4e4dc20237acb96d5449cbecc0ce
tree95435d23dad1fcc14ad69fd9ff0b038ae4797a13
parent12aa7a60b1d4de34a8c7a6a5eb35f4c9ce1b0b9c
Implement a 2D-specific chart point estimate.

The current algorithm calculates a reasonable guess for the chart point
in 3D. Using this 3D algorithm in 2D results in relatively bad guesses:

for i = 5:
affine approx chart point: 0.5 3.64977e-14
pulled back chart point:   1 0.5
for i = 6:
affine approx chart point: 0.5 0.5
pulled back chart point:   0.5 8.678e-17
for i = 7:
affine approx chart point: 0.5 1
pulled back chart point:   0.5 1

This patch adds a 2D-specific algorithm that simply averages the two
adjacent vertices (e.g., chart point 5 is assumed to be an average of
chart points 1 and 3).

This issue was reported by Juan Carlos Araujo Cabarcas, who found a test
case where the affine guess was sufficiently inaccurate that
TransfiniteInterpolationManifold::pull_back failed to find the right
point. This change fixes that test case.
source/grid/manifold_lib.cc

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.